Converting a VCF into a FASTA given some reference !
Samtools/BCFtools (Heng Li) provides a Perl script vcfutils.pl which does this, the function vcf2fq (lines 469-528) This script has been modified by others to convert InDels as well, e.g. this by David Eccles ./vcf2fq.pl -f <input.fasta> <all-site.vcf> ...Tags: Convert, VCF, FASTA, reference, genome, assembly, phase
2110 days ago