使用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
本文来自博客园,作者:Minerw,转载请注明原文链接:https://www.cnblogs.com/minerw/p/get-vcf-info-with-bcftools.html