This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tophat vs Tophat-fusion concordant alignment rate

Hi everyone!

I am working with RNA-seq paired-end data. When I run tophat-fusion on a dataset, I get mapping percentage of 70-80% for both left and right reads and the concordant pair alignment rate is about 72%. However, if I run Tophat without fusion search, I get the same mapping percentage for left and right reads, but the concodant alignment rate decreases to 30-40%. I am running both commands with the same arguments.

Does anyone know why the concordant pair alignment rate decreases with tophat-fusion? My guess would be that it should increase the concordant alignment rate as it is going to be able to map fusion events that would be unmapped by tophat.

Thank you very much!

rna-seq tophat tophat-fusion
• 2,301 views • link
written by IP IP

1 answer

Tophat fusion removes reads that map to more than two places in the genome, thus, the tophat mapping percentage can be higher because it won't remove multimappers by default

• 0 views • link
written by IP IP

Log in to answer this question.

More posts like this