Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Improve handling haploid calls #230

Open
wants to merge 4 commits into
base: main
Choose a base branch
from

Conversation

VorontsovIE
Copy link
Contributor

@VorontsovIE VorontsovIE commented Jul 29, 2020

Here are four fixes which should improve handling haploid genotypes:

  • Function revise_genotypes now skips haploid calls. Earlier it was failing due to incorrect number of genotypes (n_allele haploid genotypes instead of n_allele*(n_allele+1)/2 diploid genotypes).
  • Introduce new value NoCallReason::HaploidCall to mark alleles missing because they were missing in original haploid call. This reason is denoted by letter 'H' in RNC field.
  • FormatHelpers for Number=G[enotype] fields just skip and censor out samples which has haploid records for a site.
  • Haploid calls get genotypes A/. instead of ./A.

This patch should fix #162.

@coveralls
Copy link

coveralls commented Jul 29, 2020

Coverage Status

Coverage decreased (-0.09%) to 93.315% when pulling 47e9465 on VorontsovIE:fix_haploid_calls into ca0e9b7 on dnanexus-rnd:master.

…st allele position and to pad the second allele position with bcf_gt_missing.
@mlin
Copy link
Contributor

mlin commented Aug 13, 2020

@VorontsovIE my apologies for the slow uptake on this. It looks just right to me, but I'd like to add a new gVCF test case to lock it down against future changes. I do still have those DRAGEN gVCFs for CEU trio, but I need to get to it.........has this diff unblocked you for the time being?

Also, I really appreciate how you were able to dig into the different parts of the genotyper code. You might have seen we're considering important changes to how it works in the medium-term future, specifically around multiallelic sites. Really eager for community feedback on that from any interested parties!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Dragen JgVCF merging
3 participants