使用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

posted @   Minerw  阅读(4665)  评论(0编辑  收藏  举报
编辑推荐:
· 一个费力不讨好的项目,让我损失了近一半的绩效!
· .NET Core 托管堆内存泄露/CPU异常的常见思路
· PostgreSQL 和 SQL Server 在统计信息维护中的关键差异
· C++代码改造为UTF-8编码问题的总结
· DeepSeek 解答了困扰我五年的技术问题
阅读排行:
· 一个费力不讨好的项目,让我损失了近一半的绩效!
· 清华大学推出第四讲使用 DeepSeek + DeepResearch 让科研像聊天一样简单!
· 实操Deepseek接入个人知识库
· CSnakes vs Python.NET:高效嵌入与灵活互通的跨语言方案对比
· Plotly.NET 一个为 .NET 打造的强大开源交互式图表库
点击右上角即可分享
微信分享提示
主题色彩