This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is there an automatic way to search transcripts that align with intergenic sequences in a genome?

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,602 views • link
written by albersnz96 • albersnz96

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

• 0 views • link
written by PoGibas PoGibas

0 answers

No answers yet.

Log in to answer this question.

More posts like this