Question: Splitting .bed file by genomic location
0
gravatar for rbronste
3 months ago by
rbronste160
rbronste160 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 • 239 views
ADD COMMENTlink modified 3 months ago • written 3 months ago by rbronste160
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 3 months ago by Kevin Blighe17k

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 3 months ago by rbronste160

Any other ideas? Thanks!

ADD REPLYlink written 3 months ago by rbronste160

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 3 months ago by Kevin Blighe17k
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: 1931 users visited in the last hour