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

change 0/0 GQ0 sites to no-calls #7190

Merged
merged 1 commit into from
Apr 8, 2021
Merged

change 0/0 GQ0 sites to no-calls #7190

merged 1 commit into from
Apr 8, 2021

Conversation

kcibul
Copy link
Contributor

@kcibul kcibul commented Apr 8, 2021

@@ -365,7 +365,8 @@ private void processSampleRecordsForLocation(final long location, final Iterable
currentPositionHasVariant = true;
break;
case "0": // Non Variant Block with GQ < 10
unmergedCalls.add(createRefSiteVariantContextWithGQ(sampleName, contig, currentPosition, refAllele, 0));
// Reference calls with GQ 0 should be rendered as no-call (#271)
Copy link
Member

Choose a reason for hiding this comment

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

this will functionally include reference calls with GQ >0 but <10, right? that's intended?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Right -- post reblocking there are only GQs in increments of 10 (ie 0, 10, 20...). So we are changing the GQ 0s to no calls here, even if that was a GQ2 pre reblocking

Copy link
Contributor Author

Choose a reason for hiding this comment

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

There is no separate GQ here, for reference calls that is encoded in the state (pet)

Copy link
Contributor

@RoriCremer RoriCremer left a comment

Choose a reason for hiding this comment

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

context in comments is helpful

note to self: dig into what we are doing in createRefSiteVariantContextWithGQ

@kcibul kcibul merged commit 5387de5 into ah_var_store Apr 8, 2021
@kcibul kcibul deleted the kc_no_calls branch April 8, 2021 18:55
This was referenced Mar 17, 2023
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.

3 participants