samtools flagstat 统计三代测序比对结果的问题

运行:

samtools flagstat in.bam > out.txt

得到结果(自行添加部分注释):

12051191 + 0 in total (QC-passed reads + QC-failed reads)    #后三行primary+secondary+supplementary的总和,即多轮比对输入序列次数总和
4970675 + 0 primary    #原始比对时输入的序列总数
6887457 + 0 secondary
193059 + 0 supplementary
0 + 0 duplicates
0 + 0 primary duplicates
11997333 + 0 mapped (99.55% : N/A)    #多轮比对后总共比对到的序列次数
4916817 + 0 primary mapped (98.92% : N/A)    #成功比对到参考序列的序列数
0 + 0 paired in sequencing
0 + 0 read1
0 + 0 read2
0 + 0 properly paired (N/A : N/A)
0 + 0 with itself and mate mapped
0 + 0 singletons (N/A : N/A)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

总体比对比例 98.92% = 4916817(primary mapped) / 4970675(primary)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值