Is there an automatic way to search transcripts that align with intergenic sequences in a genome?
0
0
Entering edit mode
9.4 years ago
albersnz96 • 0

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.

RNA-Seq IGV Transcripts alignment Intergenic • 2.2k views
ADD COMMENT
1
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1830 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6