bcftools concensus -H for non-diploid/polyploids?

[peterdfields]

Hi,

I was just wondering if there is any way to apply bcftools concensus in order to output distinct haplotypes with IUPAC codes for for unphased sites when dealing with non-diploid data? Specifically, when I set -H 3 or -H 4 the approach works fine, but if I try -H 3pIu or -H 4pIu I get an error like the following:
Error: Could not parse --haplotype 4pIu, expected numeric argument
Maybe there's a simple fix to allow the same behavior as -H 1pIu or 2pIu?

[pd3]

This is now supported. Thank you for the issue

[peterdfields]

@pd3 thank you for the modification!