Skip to content
Snippets Groups Projects
Find file Normal view History Permalink
excluded_regions_processing.sh 703 B
Newer Older
Venkat Malladi's avatar
First pass at calling enhancers.
Venkat Malladi committed
1 
# Make a universe of excluded regions of the genome
Venkat Malladi's avatar
Update code for protein coding genes tss.  
Venkat Malladi committed
2 3 4 
cd MakeGencodeTSS
./make_TSS_file_from_annotation_with_confidence_better.sh ../gencode.v19.annotation.gtf ../biotypes.txt
gtf2bed  < MakeGencodeTSS/gencode.v19.annotation_capped_sites.gff > gencode.v19.annotation_capped_sites.bed
Venkat Malladi's avatar
First pass at calling enhancers.
Venkat Malladi committed
5 6 7 8 9 10 11 12 13 

## Gencode Protein Coding genes
bedtools slop -i gencode.v19.annotation_capped_sites.bed -g hg19.chrom.sizes -b 3000 | cut -f1,2,3 > gencode_tss_3kb_flanking.bed

## H3K4me3
bedtools slop -i H3K4me3_filtered_peaks.bed -g hg19.chrom.sizes -b 3000 > H3K4me3_3kb_flanking.bed

## Merge TSS and H3K4me3 peaks
bedops --everything gencode_tss_3kb_flanking.bed H3K4me3_3kb_flanking.bed | sort-bed - > excluded_3kb_flanking.bed