美文网首页生物信息学
[bcftools] 过滤vcf

[bcftools] 过滤vcf

作者: Silver_42ac | 来源:发表于2021-03-02 10:38 被阅读0次
#-i/e, --include/--exclude <expr>    
#select/exclude sites for which the expression is true (see man page for details)

bcftools view -i "DP>10" -o  ${vcf}.gz  -Oz ${vcf}

expressions

报错:

Both INFO/DP and FORMAT/DP exist, which one do you want?

这里做下修改

ls *.raw.snv_indel.vcf > all.vcf.lst

while read vcf ;do
bcftools view -i "INFO/DP>10" -o  ${vcf}.gz  -Oz ${vcf}
bcftools index ${vcf}.gz

done < all.vcf.lst

bcftools merge  -Ov  -o  all_sample.merge.vcf  *.gz 

参考:
http://www.htslib.org/doc/1.0/bcftools.html#expressions

相关文章

网友评论

    本文标题:[bcftools] 过滤vcf

    本文链接:https://www.haomeiwen.com/subject/lbetqltx.html