Practical Haplotype Graph v2 not finding correct paths

Question

[beantkapoor16]

Hello, I am using Practical Haplotype Graph v2.2.85.134 to build a pangenome graph using six diploid plant species (12 haplotypes). I was able to go through their Build and Load module. Then, I simulated 10 million WGS Illumina Novaseq paired end reads from two of my haplotypes (australasica_primary and fallglo_primary) and mapped it to the pangenome using their Imputation module. The idea was that the most haplotype paths in the imputed vcf file should originate from those two haplotypes however that's not the case. Below are the number of haplotype paths from each of my haplotypes.

• 5,232 australasica_alternate
• 325 australasica_primary
• 11,390 australis_alternate
• 14,105 australis_primary
• 1,439 fallglo_alternate
• 10,086 fallglo_primary
• 12,617 fortune_alternate
• 5 inodora_alternate
• 17 inodora_primary

I understand there are a lot of moving parts to this pipeline and I would be happy to provide more details if requested. Any suggestion is highly appreciated. Thank you.

Answer 1

[pjb39]

There have been a couple of bugs found that affect imputation accuracy, one in build-kmer-index and one in find-paths. The find-paths bug only affects haploid path finding. Those should be fixed in the next few days. I will post here when the fix is available.

[beantkapoor16]

Hello, just wanted to check in if the code has been fixed yet? Thank you.