I'm using Tophat to align RNA-Seq reads to a genome. I want to know how many reads aligned, but I didn't know about the -g option when I ran Tophat, so the normal commands like "samtools -c accepted_hits.bam" isn't giving me the number of unique reads in the file, it's giving me the total number of alignments. Is there any way I can get this information without re-running Tophat?
Question: Get # of unique mapped reads from Tophat output
5.3 years ago by
colin.kern • 920
colin.kern • 920 wrote:
ADD COMMENT • link •
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 1643 users visited in the last hour