Question: Splitting .bed file by genomic location
0
gravatar for rbronste
9 months ago by
rbronste170
rbronste170 wrote:

Hi, pretty simple question - wondering about the easiest way to take a single .bed file and split the ranges into those at defined mm10 promoters and those that are intergenic/intronic. With the promoter being +/- 2kb from the TSS? Thanks!

bed chip-seq promoter R • 328 views
ADD COMMENTlink modified 9 months ago • written 9 months ago by rbronste170
1

Do you already have the co-ordinates of the defined promoters?; or are you just saying that, generally, you are assuming that they are +/- 2kb from the TSS?

ADD REPLYlink written 9 months ago by Kevin Blighe30k

I am just generally trying to cover a large enough area to exclude most promoter elements, maybe +/-2.5kb is even more appropriate.

ADD REPLYlink written 9 months ago by rbronste170

Any other ideas? Thanks!

ADD REPLYlink written 9 months ago by rbronste170

You could download a reference template of genes and promoter regions from the UCSC Table Browser in BED format, and then use bedtools intersect, or something like that, in order to overlap your current regions to the actual defined promoters.

How can I fetch promoter sequence upstream of a gene?

Kevin

ADD REPLYlink written 9 months ago by Kevin Blighe30k
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: 1792 users visited in the last hour