Wondering about a straightforward way to take BED files with peak locations for ChIP-seq data and see where they fall against a list of highly expressed genes (in terms of proximity) from RNA-seq data? Thanks
Firstly, you need to define the distance, such as promoter, enhancer; then you can judge which peaks are located in the gene feature aboved-mentioned with custome script.
Probably a more elegant solution available, but:
Grab gene locations in BED format. Generally you can do this fairly easily using the UCSC Table Browser and their identifiers option.
Use bedops closest-feature with the --dist flag to calculate distances. Make sure to sort files correctly before running this step using bedops sort-bed function or you'll get weird results.
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy