Build BED files for genomic regions enrichment analysis
Specific analysis of H3K27ac status in super-enhancer regions
It is time-consuming for analysis large datasets.
It needs about 3 seconds per genome region and per tissue/cell type. What should be mentioned is that there maybe more than one dataset for each tissue/cell type.
For test, we submited the "Example BED" files with 8 human (hg19) genome regions, and selected all the 99 tissue/cell types. Then, it takes more than 40 minutes to get the result of overlap table and heeatmap for average status of H3K27ac.
Example for Regions:
#bed file description
#<chrom> <start> <end> <region discription>
chr3 181403691 181476186 region_1
chr3 181403893 181455295 region_2
chr3 181404243 181455400 region_3