HipHap: Haplotype Assignment and Confidence Scoring for Diploid Reference Genomes
Heinz, J. M.; Marin, M. G.; Meyerson, M.; Li, H.
Show abstract
Diploid genome assemblies are now routinely available, but most read aligners were designed for haploid references, which have long been the gold standard. When reads are aligned to a diploid assembly, the aligner sees two nearly identical alignments to either haplotype, thus reducing the mapping quality (MapQ) score to reflect this ambiguity. This can cause downstream tools to discard reads from easily mappable regions. Here, we present HipHap (HIgh- Performance HAPlotype assigner) to resolve this issue by aligning reads to each haplotype assembly separately and assigning each read to its best-supported haplotype. HipHap increased the fraction of reads mapping to the diploid genome with high MapQ, outperforming alignment to either parental haplotype. Additionally, we introduce a haplotype assignment quality score (HapQ) in HipHap to quantify confidence in the haplotype of origin of a read. HipHap is implemented in Rust, supports SAM, BAM, CRAM, and PAF formats, and is freely available at: https://github.com/jheinz27/hiphap.
Matching journals
The top 5 journals account for 50% of the predicted probability mass.
Similar papers in this journal
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.