Question: Splitting .bed file by genomic location
0
gravatar for rbronste
15 months ago by
rbronste240
rbronste240 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 • 449 views
ADD COMMENTlink modified 15 months ago • written 15 months ago by rbronste240
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 15 months ago by Kevin Blighe41k

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 15 months ago by rbronste240

Any other ideas? Thanks!

ADD REPLYlink written 15 months ago by rbronste240

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 15 months ago by Kevin Blighe41k
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: 2031 users visited in the last hour