基因组学

bam 提取与过滤

2021-01-13  本文已影响0人  Silver_42ac

`

过滤掉 secondary_alignment or supplementary

sambamba view -q -f bam --filter='not (secondary_alignment or supplementary)' xxxx.minimap2.hg19.xxx.bam > xxx.minimap2.hg19.chr10.primary_aln.nosupp.bam
`

`

过滤质量只低于20的,指定区域

sambamba view -q -f bam --filter='mapping_quality >20' x.minimap2.hg19.bam 10:1-1234567 > xx.bam
or
samtools view -b -m 20 x.minimap2.hg19.bam 10:1-1234567 > xx.bam
`

参考
https://github.com/biod/sambamba/wiki/[sambamba-view]-Filter-expression-syntax
https://broadinstitute.github.io/picard/explain-flags.html

上一篇下一篇

猜你喜欢

热点阅读