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

Fix error in genotype given alleles mode when input alleles have genotypes #5341

Merged
merged 2 commits into from
Oct 23, 2018

Conversation

cwhelan
Copy link
Member

@cwhelan cwhelan commented Oct 22, 2018

Fixes #5336

HaplotypeCallerGenotypingEngine.replaceWithSpanDelVC replaces the given alleles VariantContext objects with new VCs in which the alleles have been replaced with star alleles to represent spanning deletions. However, if the input variant records had genotypes, the original (non-star) deletion alleles were still being retained, causing a downstream error. This fixes that issue.

@cwhelan cwhelan requested a review from ldgauthier October 22, 2018 17:21
Copy link
Contributor

@ldgauthier ldgauthier left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for the quick turnaround. This fix should make it into today's release.

@ldgauthier
Copy link
Contributor

...I mean tomorrow's release.

@droazen droazen self-assigned this Oct 22, 2018
Copy link
Contributor

@droazen droazen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@cwhelan While we're waiting on travis to come back to life, could you add a direct unit test for HaplotypeCallerGenotypingEngine.replaceWithSpanDelVC() to check the genotypes in the returned VC?

@cwhelan
Copy link
Member Author

cwhelan commented Oct 23, 2018

@droazen added a quick unit test that failed without my change but passes now. look ok?

@codecov-io
Copy link

Codecov Report

Merging #5341 into master will increase coverage by 0.004%.
The diff coverage is 100%.

@@               Coverage Diff               @@
##              master     #5341       +/-   ##
===============================================
+ Coverage     86.821%   86.825%   +0.004%     
- Complexity     30197     30199        +2     
===============================================
  Files           1845      1845               
  Lines         139825    139852       +27     
  Branches       15412     15412               
===============================================
+ Hits          121398    121427       +29     
- Misses         12818     12819        +1     
+ Partials        5609      5606        -3
Impacted Files Coverage Δ Complexity Δ
...aller/HaplotypeCallerGenotypingEngineUnitTest.java 89.501% <100%> (+0.769%) 30 <1> (+1) ⬆️
...plotypecaller/HaplotypeCallerGenotypingEngine.java 80% <100%> (+0.126%) 39 <0> (ø) ⬇️
...e/hellbender/engine/spark/SparkContextFactory.java 71.233% <0%> (-2.74%) 11% <0%> (ø)
...oadinstitute/hellbender/utils/gcs/BucketUtils.java 81.098% <0%> (+0.61%) 42% <0%> (ø) ⬇️
...ypecaller/AssemblyBasedCallerGenotypingEngine.java 89.401% <0%> (+1.382%) 88% <0%> (+1%) ⬆️

Copy link
Contributor

@droazen droazen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

👍 New test looks good, merging

@droazen droazen merged commit 5d95d60 into master Oct 23, 2018
@droazen droazen deleted the cw_gga_allele_not_in_variant_context branch October 23, 2018 15:15
EdwardDixon pushed a commit to EdwardDixon/gatk that referenced this pull request Nov 9, 2018
…types (broadinstitute#5341)

* fix issue of alleles not being replaced with star alleles when the given allele variants have genotypes

* add unit test to replaceWithSpanDelVC
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.

HaplotypeCaller GGA mode crashes with certain spanning deletions
4 participants