Doubt samtools flagstat
0
0
Entering edit mode
6 weeks ago
Etienne • 0

I'd like to see the percentage of how many sequences align with my decrementing sequence and I've come to this sample table. But wanted to know what use? The mapped (80.94% : N/A)or properly paired (0.06% : N/A) percentage?

1036193 + 0 in total (QC-passed reads + QC-failed reads)
197709 + 0 primary
0 + 0 secondary
838484 + 0 supplementary
0 + 0 duplicates
0 + 0 primary duplicates
838710 + 0 mapped (80.94% : N/A)
226 + 0 primary mapped (0.11% : N/A)
197709 + 0 paired in sequencing
124 + 0 properly paired (0.06% : N/A)
226 + 0 with itself and mate mapped
0 + 0 singletons (0.00% : N/A)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

paired Samtools flagstat read • 276 views
1
Entering edit mode

I'd start by examining why the properly paired % is so low.

1
Entering edit mode

Since 80% of reads are mapped it is likely that the reads are out of sync if the R1/R2 files were processed/trimmed independently. Fix using repair.sh from BBMap in that case.

0
Entering edit mode

I am comparing a variant of the ITS region that is not very common.