Question: Get # of unique mapped reads from Tophat output
3
gravatar for colin.kern
4.2 years ago by
colin.kern420
United States
colin.kern420 wrote:

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?

rna-seq tophat • 2.0k views
ADD COMMENTlink modified 4.2 years ago by Devon Ryan92k • written 4.2 years ago by colin.kern420
3
gravatar for Devon Ryan
4.2 years ago by
Devon Ryan92k
Freiburg, Germany
Devon Ryan92k wrote:

samtools view -c -F 256 accepted_hits.bam

ADD COMMENTlink written 4.2 years ago by Devon Ryan92k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1096 users visited in the last hour