-
Notifications
You must be signed in to change notification settings - Fork 250
[query/combiner] Fix import_gvcf_interval's default PGT handling #13364
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
[query/combiner] Fix import_gvcf_interval's default PGT handling #13364
Conversation
PGT is a call field by default for import_vcf, the now former import_gvcfs (old combiner start point). Make import_gvcf_interval also treat PGT as a call field.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we have a GVCF with PGT? Can you test for this?
We do! Added a quick assertion to an extant test.
@@ -39,6 +39,11 @@ def test_combiner_works(): | |||
out = os.path.join(tmpdir, 'out.vds') | |||
hl.vds.new_combiner(temp_path=tmpdir, output_path=out, gvcf_paths=paths, intervals=parts, reference_genome='GRCh38').run() | |||
comb = hl.vds.read_vds(out) | |||
|
|||
# see https://github.com/hail-is/hail/issues/13367 for why these assertions are here | |||
assert 'LPGT' in comb.variant_data |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think you need variant_data.entry
or something
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
🤦 thanks!
PGT is a call field by default for import_vcf, the now former import_gvcfs (old combiner start point). Make import_gvcf_interval also treat PGT as a call field.