This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNA-seq: indexing without sorting

hello there,

In RNA-seq tuxedo pipeline I was taught to directly samtools index the bam files, obtained from the tophat alignment, without indexing. and in DNA-seq I was taught that first, we sort the bam files then only index in samtools. what's the difference? and reason to do this?

thank you.

hello rna-seq samtools index and sort
• 1,154 views • link
updated by WouterDeCoster WouterDeCoster • written by Rakesh Tiwari • Rakesh Tiwari

1 answer

It depends on the aligner. In this case, Tophat (as far as I know) creates a sorted bam by default. As such, there is no additional sorting necessary anymore. I believe that samtools index only works on coordinate sorted bams. By the way, TopHat is no longer a recommended aligner, and you should use something like STAR (for alignment) or salmon (for rapid pseudo-alignment).

• 0 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this