构建consensus序列

2020-04-27 本文已影响0人曲凉不见

bcftools view sample.vcf -Oz -o sample.vcf.gz
bcftools index sample.vcf.gz

samtools faidx ref.fna refgenome:2253-3709 | bcftools consensus sample.vcf.gz > consensus.sample.fa