PlantandFoodResearch / MCHap

Polyploid micro-haplotype assembly using Markov chain Monte Carlo simulation.
MIT License
18 stars 3 forks source link

Consider option to calculate log-likelihoods ignoring dosage for VCF output #97

Open timothymillar opened 3 years ago

timothymillar commented 3 years ago

This can be achieved by treating an allelic phenotype as a genotype (i.e. ploidy = number of unique alleles). Calculating this over all possible genotypes would give genotype likelihoods that are invariant to dosage bias. This would not be used in the haplotype assembly, only as a VCF output for downstream analysis and potential dosage correction.