美文网首页funny生物信息plink
2020-01-14 如何合并VCF文件

2020-01-14 如何合并VCF文件

作者: YIYI_91fd | 来源:发表于2020-01-14 20:33 被阅读0次

因为一直在练手GWAS分析,所以刚开始Call SNP的时候是单个个体call的,后面就需要把这些不同vcf合并。

琢磨两天后,发现bcftools可以搞定。

首先对每个个体分别压缩。再建立索引。

bgzip XZ0201.call.vcf

tabix -p vcf XZ0201.call..vcf.gz

最后,bcftools合并即可。

bcftools -merge XZ0201.call.vcf.gz XZ0209.call.vcf.gz XZ0215.call.vcf.gz XZ0224.call.vcf.gz XZ0230.call.vcf.gz XZ0410.call.vcf.gz XZ0416.call.vcf.gz XZ0423.call.vcf.gz > merge.vcf

相关文章

网友评论

    本文标题:2020-01-14 如何合并VCF文件

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