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

Maximum alt alleles guardrail [VS-1530] #9033

Open
wants to merge 1 commit into
base: ah_var_store
Choose a base branch
from

Conversation

mcovarr
Copy link
Collaborator

@mcovarr mcovarr commented Nov 5, 2024

Extract and WeightedSplitIntervals don't deal well with arbitrarily large numbers of alt alleles, so add an (overridable) guardrail of 1000 alt alleles max if the user neglects to specify the 100 alt alleles max recommended in our documentation for AoU callsets.

@mcovarr mcovarr changed the base branch from master to ah_var_store November 5, 2024 13:27
@@ -135,9 +135,9 @@ Once the VAT table is created and a tsv is exported, the AoU research workbench
- Specify the same `call_set_identifier`, `dataset_name`, `project_id`, `extract_table_prefix`, and `interval_list` that were used in the `GvsPrepareRangesCallset` run documented above.
- Specify the `interval_weights_bed` appropriate for the PGEN / VCF extraction run you are performing. `gs://gvs_quickstart_storage/weights/gvs_full_vet_weights_1kb_padded_orig.bed` is the interval weights BED used for Quickstart.
- For both `GvsExtractCallset` and `GvsExtractCallsetPgenMerged`, select the workflow option "Retry with more memory" and choose a "Memory retry factor" of 1.5
- `GvsExtractCallset` currently defaults to 1000 alt alleles maximum, which means that any sites having more than that number of alt alleles will be dropped. Make sure to specify the appropriate `maximum_alternate_alleles` value (currently 100).
Copy link
Collaborator

Choose a reason for hiding this comment

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

I'm confused. "currently 100?"

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

This is carried over from line 140. The value of max alt alleles we were told to use for Echo was 100, but in my testing we seem to be able to go to 1000 at least without WeightedSplitIntervals breaking down. I don't know if variants with > 100 alt alleles are actually useful though.

@mcovarr mcovarr changed the title Maximum alt alleles guardrail [VS-1491] Maximum alt alleles guardrail [VS-1530] Nov 5, 2024
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.

2 participants