Dear bcftools team,
this was brought up in a chat with @pd3 about a month ago — it would be cool to be able to call variants at ploidy levels above 2 with the bcftools mpileup | bcftools call pipeline! I am aware of freebayes and gatk but I got to like bcftools variant calls for their balance of sensitivity and not too heavy assumptions of the input data which altogether works well for non-human population datasets. Is that something that is on the roadmap or at least feasible to implement?
Cheers, Nikita