Question: Intersect Gene Annotation With Specific Position Or Genomic Interval
0
gravatar for Nicolas Rosewick
7.2 years ago by
Belgium, Brussels
Nicolas Rosewick9.2k wrote:

Hi,

I've several genomic interval and I want to check if they are overlapping with known gene. I've a gtf file with the coordinates of gene exons. My idea was to use intersectBed from bedtools but I've a little problem with small genomic interval that are are overlapping intron coordinates and not exons ( it do ot report me the gene where this interval is). Is it possible to specifiy to intersectBed to take into account introns ? or is there an another tool ?

Thanks

N.

bedtools intersect • 3.1k views
ADD COMMENTlink modified 7.2 years ago by Chris Miller21k • written 7.2 years ago by Nicolas Rosewick9.2k
0
gravatar for Chris Miller
7.2 years ago by
Chris Miller21k
Washington University in St. Louis, MO
Chris Miller21k wrote:

Instead of intersecting with some monolithic gene track containing everything, you're going to want to intersect with a track containing only exons.

Since I don't know exactly what your data looks like, I can't tell you exactly how to accomplish this. If you don't have coordinates for specific exons in your data, you can download such tracks from UCSC genome browser or Ensembl easily enough.

ADD COMMENTlink written 7.2 years ago by Chris Miller21k

I've ensembl gtf file

ADD REPLYlink written 7.2 years ago by Nicolas Rosewick9.2k

I'm not immediately familiar with that format, but it probably contains lines labelled "intron", "cds_exon", "rna", etc. Grep out the lines you want, do your intersection, then collapse by gene name if necessary

ADD REPLYlink written 7.2 years ago by Chris Miller21k
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: 1271 users visited in the last hour