參考此網(wǎng)站改善的:http://www.itdecent.cn/p/a23d6f1226a1
## 這個(gè)比較簡(jiǎn)單好用。
## 第一步
cd /datapool/gesong/mengqinglin/sh/
bgzip -c 你的.vcf > 你的.vcf.gz
tabix -p vcf 你的.vcf.gz
##第二步
cat 你的參考基因組.fasta \
| bcftools consensus \
你的.vcf.gz > 你的.fa
參考此網(wǎng)站改善的:http://www.itdecent.cn/p/a23d6f1226a1
## 這個(gè)比較簡(jiǎn)單好用。
## 第一步
cd /datapool/gesong/mengqinglin/sh/
bgzip -c 你的.vcf > 你的.vcf.gz
tabix -p vcf 你的.vcf.gz
##第二步
cat 你的參考基因組.fasta \
| bcftools consensus \
你的.vcf.gz > 你的.fa