For past two hours I tried to browse around for samtools flagstat output information, unfortunately could not find much related to my problem - this is my samtool flagstat output for single end reads - does anyone think that the tophat alignment did not work well. Please let me know. Any possible replies would be greatly appreciated. Thanks in advance.
$ samtools flagstat accepted_hits.bam
18109074 + 0 in total (QC-passed reads + QC-failed reads) 0 + 0 duplicates 18109074 + 0 mapped (100.00%:nan%) 0 + 0 paired in sequencing 0 + 0 read1 0 + 0 read2 0 + 0 properly paired (nan%:nan%) 0 + 0 with itself and mate mapped 0 + 0 singletons (nan%:nan%) 0 + 0 with mate mapped to a different chr 0 + 0 with mate mapped to a different chr (mapQ>=5)
1 answer
It has to be multi-reads.
try
samtools view accepted_hits.bam | cut -f1 | sort | uniq | wc -l
this should give you the count of unique reads
•
written
by
Wen.Huang
Log in to answer this question.
-
Samtools: How can I extract properly-paired QC-passed reads instead of extracting only properly-pai…
written by bioinfo8
-
BWA mem paired end vs single end shows unusual flagstat summary
written by rrdavis •
-
RNA seq samtools flagstat
written by pshubhamoy •
-
Tophat mapping less reads (and less properly paired) when more real -r is specified?
written by manekineko
-
Samtools Flagstat Output Format
written by thsay •
-
samtool flagstat question
written by alex
-
How to extract properly paired alignments from a bam file using samtools?
written by pedrodcb
-
How To Check That Tophat Works Well
written by M K
-
Thorough Documentation Of Samtools Flagstat
written by Biomonika (Noolean)
-
Does The New Tophat2 - Bowtie2 Pipeline Really Map 100% Of The Rna-Seq Reads ? Is This Real ?
written by biorepine
