Do we need a step to move non-bialellic SNPs from the VCF?

[janaobsteter]

In the ancestral allele inference, we filter the variants for SNPs, but we do allow multi-allelic SNPs since that is not a problem for the ancestral inference (it's just an allele count). However, some rule might fall if there are multiallelic variants in the VCF and the INFO file (for example, the get_major rule in prepare_files_for_tsinfer.smk.

[gregorgorjanc]

We discussed in one of the meetings how we could convert multi-allelic sites to multiple bi-allelic sites. My understanding was that bcftools enables this conversion. Would that help?

[janaobsteter]

For now, let's filter the VCF to only keep biallelic SNPs to see how it behaves and then we can progress with adding other types of variants in the future.