构建一致性序列方法如下:
1.将vcf文件压缩
bgzip -c raw.filtered.snp.vcf > raw.filtered.snp.vcf.gz
2.建立索引
tabix raw.filtered.snp.vcf.gz
3.构建一致性序列
bcftools consensus -f Ghirsutum_527_v2.0.fa raw.filtered.snp.vcf.gz > new.fa
351 篇文章
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!