Peak ident near highly expressed genes
1
0
Entering edit mode
4.2 years ago
rbronste ▴ 400

Hi,

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

Rob

ChIP-Seq RNA-Seq bedops • 769 views
ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
1
Entering edit mode
4.2 years ago
Sinji ★ 3.1k

Probably a more elegant solution available, but:

  1. Grab gene locations in BED format. Generally you can do this fairly easily using the UCSC Table Browser and their identifiers option.

  2. 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.

ADD COMMENT

Login before adding your answer.

Traffic: 2264 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6