$ samtools flagstat
Usage: samtools flagstat [options] <in.bam>
--input-fmt-option OPT[=VAL]
Specify a single input file format option in the form
of OPTION or OPTION=VALUE
-@, --threads INT
Number of additional threads to use [0]
$ samtools flagstat -@ 8 EV-1_HISAT2.bam
59306432 + 0 in total (QC-passed reads + QC-failed reads)
4546096 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
55162761 + 0 mapped (93.01% : N/A) #匹配到reference的reads百分比（55162761/59306432）
54760336 + 0 paired in sequencing #pair數(shù)目
27380168 + 0 read1
27380168 + 0 read2
48652166 + 0 properly paired (88.85% : N/A)
49197764 + 0 with itself and mate mapped #pair reads都成功匹配到reference sequence的對(duì)數(shù)
1418901 + 0 singletons (2.59% : N/A) #pair reads只有一條匹配到reference sequence的條數(shù)
106592 + 0 with mate mapped to a different chr #paired reads分別匹配到不同染色體
82634 + 0 with mate mapped to a different chr (mapQ>=5)
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡(jiǎn)書系信息發(fā)布平臺(tái)，僅提供信息存儲(chǔ)服務(wù)。
