Skip to content

Commit b87bd2e

Browse files
committed
do alleles field validation only if it exists on ht
1 parent bfde429 commit b87bd2e

File tree

20 files changed

+69
-0
lines changed

20 files changed

+69
-0
lines changed

v03_pipeline/lib/misc/validation.py

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -31,6 +31,9 @@ def validate_allele_type(
3131
**_: Any,
3232
) -> None:
3333
ht = t.rows() if isinstance(t, hl.MatrixTable) else t
34+
if not hasattr(ht, 'alleles'):
35+
return
36+
3437
ht = ht.filter(
3538
dataset_type.invalid_allele_types.contains(
3639
hl.numeric_allele_type(ht.alleles[0], ht.alleles[1]),
Lines changed: 63 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,63 @@
1+
import unittest
2+
from unittest.mock import patch
3+
4+
import hail as hl
5+
6+
from v03_pipeline.lib.model import ReferenceGenome
7+
from v03_pipeline.lib.reference_datasets.reference_dataset import ReferenceDataset
8+
9+
TEST_GNOMAD_SVS_RAW_HT = (
10+
'v03_pipeline/var/test/reference_datasets/raw/gnomad_svs_from_vcf.ht'
11+
)
12+
13+
14+
class GnomadSVsTest(unittest.TestCase):
15+
@patch('v03_pipeline.lib.reference_datasets.gnomad_svs.vcf_to_ht')
16+
def test_gnomad_svs(self, mock_vcf_to_ht):
17+
mock_vcf_to_ht.return_value = hl.read_table(TEST_GNOMAD_SVS_RAW_HT)
18+
ht = ReferenceDataset.gnomad_svs.get_ht(ReferenceGenome.GRCh38)
19+
self.assertEqual(
20+
ht.collect(),
21+
[
22+
hl.Struct(
23+
KEY='gnomAD-SV_v3_BND_chr1_1a45f73a',
24+
AF=0.11413399875164032,
25+
AC=8474,
26+
AN=74246,
27+
N_HET=8426,
28+
N_HOMREF=28673,
29+
),
30+
hl.Struct(
31+
KEY='gnomAD-SV_v3_BND_chr1_3fa36917',
32+
AF=0.004201000090688467,
33+
AC=466,
34+
AN=110936,
35+
N_HET=466,
36+
N_HOMREF=55002,
37+
),
38+
hl.Struct(
39+
KEY='gnomAD-SV_v3_BND_chr1_7bbf34b5',
40+
AF=0.03698499873280525,
41+
AC=3119,
42+
AN=84332,
43+
N_HET=3115,
44+
N_HOMREF=39049,
45+
),
46+
hl.Struct(
47+
KEY='gnomAD-SV_v3_BND_chr1_933a2971',
48+
AF=0.3238990008831024,
49+
AC=21766,
50+
AN=67200,
51+
N_HET=21616,
52+
N_HOMREF=11909,
53+
),
54+
hl.Struct(
55+
KEY='gnomAD-SV_v3_DUP_chr1_01c2781c',
56+
AF=0.0019970000721514225,
57+
AC=139,
58+
AN=69594,
59+
N_HET=139,
60+
N_HOMREF=34658,
61+
),
62+
],
63+
)
Binary file not shown.
Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,3 @@
1+
This folder comprises a Hail (www.hail.is) native Table or MatrixTable.
2+
Written with version 0.2.133-4c60fddb171a
3+
Created at 2025/03/04 14:19:46

v03_pipeline/var/test/reference_datasets/raw/gnomad_svs_from_vcf.ht/_SUCCESS

Whitespace-only changes.

0 commit comments

Comments
 (0)