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
4.2 years ago by
colin.kern • 420
colin.kern • 420 wrote:
ADD COMMENT • link •
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 1096 users visited in the last hour