Hi Everyone!
I am writing down some stats for my RNAseq data.
Briefly, I have used Tophat2+samtools+HTSeq count for mapping and gene counting (PE, GRCh37).
When i check the samtools flagstat everything looks terrific...great mapping rates as well as pairing between reads.
I was trying to get my unique mapping rates stats:
I've been using the accepted_hits.bam file and run
samtools view accepted_hits.bam | cut -f1 | sort | uniq | wc -l
My naive question is:
the number I get, is it considering the pair or the number of reads???
Anyone has a better solution for counting unique mapping rate for PE?
Thanks!
Manu
I am writing down some stats for my RNAseq data.
Briefly, I have used Tophat2+samtools+HTSeq count for mapping and gene counting (PE, GRCh37).
When i check the samtools flagstat everything looks terrific...great mapping rates as well as pairing between reads.
I was trying to get my unique mapping rates stats:
I've been using the accepted_hits.bam file and run
samtools view accepted_hits.bam | cut -f1 | sort | uniq | wc -l
My naive question is:
the number I get, is it considering the pair or the number of reads???
Anyone has a better solution for counting unique mapping rate for PE?
Thanks!
Manu
Comment