vcf格式转fasta格式

参考此网站改善的:https://www.jianshu.com/p/a23d6f1226a1

## 这个比较简单好用。

## 第一步

cd /datapool/gesong/mengqinglin/sh/

bgzip -c 你的.vcf > 你的.vcf.gz

tabix -p vcf 你的.vcf.gz

##第二步

cat 你的参考基因组.fasta \

| bcftools consensus \

你的.vcf.gz > 你的.fa

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。