Question: Is there an automatic way to search transcripts that align with intergenic sequences in a genome?
0
gravatar for albersnz96
4.4 years ago by
albersnz960
Mexico
albersnz960 wrote:

I have a problem with the IGV program. I am comparing the information of all the annotated genes for a yeast (.gtf file) and the mapped RNA-seq reads with its genome performed with Tophat (the genome is a .fa file and the mapped sequences is a .bam file). I am looking for transcripts that align with intergenic sequences to correct the annotations in the .gtf file; therefore, I would like to check the number and localization of each transcript that align with an intergenic sequence. Could somebody tell me if there is a way of doing this analysis in an automatic form using IGV or if I need the use of another Bioinformatics tool?. Thank you for your help.

ADD COMMENTlink written 4.4 years ago by albersnz960
1

Try bedtools intersect (what you want is to extract regions with reads that doesn't intersect annotated genes).

ADD REPLYlink modified 4.4 years ago • written 4.4 years ago by PoGibas4.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: 1253 users visited in the last hour