Question: Peak ident near highly expressed genes
0
gravatar for rbronste
24 months ago by
rbronste240
rbronste240 wrote:

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

rna-seq chip-seq bedops • 550 views
ADD COMMENTlink modified 24 months ago • written 24 months ago by rbronste240

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 REPLYlink written 24 months ago by Ben50
1
gravatar for Sinji
24 months ago by
Sinji2.8k
UT Southwestern Medical Center
Sinji2.8k wrote:

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 COMMENTlink modified 24 months ago • written 24 months ago by Sinji2.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1707 users visited in the last hour