Open hyanwong opened 5 months ago
At the moment, the write_vcf() function simply outputs the 0th allele from variant.alleles (i.e. the ancestral state) as the REF. If we have defined a reference sequence, then presumably we should be using the actual REF instead, if it exists.
write_vcf()
Either way, the behaviour should be documented.
At the moment, the
write_vcf()
function simply outputs the 0th allele from variant.alleles (i.e. the ancestral state) as the REF. If we have defined a reference sequence, then presumably we should be using the actual REF instead, if it exists.Either way, the behaviour should be documented.