MCHap icon indicating copy to clipboard operation
MCHap copied to clipboard

Option to drop uncalled haplotypes

Open timothymillar opened this issue 1 year ago • 1 comments

MCHap can report haplotypes that have not been called in the posterior mode of any individual (i.e., and allele count of 0). This is necessary when reporting statistics over the posterior (such as the posterior mean allele frequencies) but is not needed if we just want genotype calls. It would be nice to have an option to exclude uncalled alleles from the output VCF file. E.g. a simple flag like --drop-uncalled-haplotypes. This could potentially raise an error if used combination with an option to report posterior summary stats.

timothymillar avatar Feb 23 '23 00:02 timothymillar