使用bcftools查看vcf文件所含位点数量信息


shell命令:

bcftools +counts $filename

其中counts是bcftools的一个插件(plugin)
后面接上的$filename即为需要查看的vcf文件

示例:
这里使用千牛基因组Run9数据做例子
shell命令和运行结果:

> bcftools +counts ChrMT-Run9-TAUIND-raw-toDistribute.vcf.gz

Number of samples: 6191
Number of SNPs:    5739
Number of INDELs:  930
Number of MNPs:    0
Number of others:  0
Number of sites:   6242

REF

https://samtools.github.io/bcftools/bcftools.html


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM