samtools view中关于flag过滤的几个特点

以下这段翻译自SAMv1.pdf

 

__________________________

 

以下几个选项计算结果一样

-c -f1 =-c -F12

-c -f2 =-c -f2 -F4 = -c -f2 -F12

-c -F12=-c -F524  # 0x4 + 0x8 + 0x200

-c -F4 = -c -F516

-c -F260=-c -F772 # 0x4 + 0x100 + 0x200

-c -F268=-c -F780 # 0x4 + 0x8 + 0x100 + 0x200

__________________________

samtools flagstat中几个指标的计算方法

$ samtools flagstat H1650RNA_S12.hg19.bam
28400342 + 0 in total (QC-passed reads + QC-failed reads)　　#  -c
2924468 + 0 secondary 　　　　　　　　　　　　　　#  -c -f0x100
0 + 0 supplementary
0 + 0 duplicates
28188350 + 0 mapped (99.25%:-nan%)　　　　　　# -c -F4
25475874 + 0 paired in sequencing
12737937 + 0 read1
12737937 + 0 rea