Skip to content

Commit b37401d

Browse files
authored
fix MITO clinvar crdq (#735)
1 parent 6537a44 commit b37401d

File tree

2 files changed

+11
-6
lines changed

2 files changed

+11
-6
lines changed

v03_pipeline/lib/model/cached_reference_dataset_query.py

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -26,10 +26,11 @@ def access_control(self) -> AccessControl:
2626
return AccessControl.PRIVATE
2727
return AccessControl.PUBLIC
2828

29-
@property
30-
def dataset(self) -> str | None:
29+
def dataset(self, dataset_type: DatasetType) -> str | None:
3130
return {
32-
CachedReferenceDatasetQuery.CLINVAR_PATH_VARIANTS: 'clinvar',
31+
CachedReferenceDatasetQuery.CLINVAR_PATH_VARIANTS: 'clinvar_mito'
32+
if dataset_type == DatasetType.MITO
33+
else 'clinvar',
3334
CachedReferenceDatasetQuery.GNOMAD_CODING_AND_NONCODING_VARIANTS: 'gnomad_genomes',
3435
CachedReferenceDatasetQuery.GNOMAD_QC: 'gnomad_qc',
3536
CachedReferenceDatasetQuery.HIGH_AF_VARIANTS: 'gnomad_genomes',

v03_pipeline/lib/tasks/reference_data/updated_cached_reference_dataset_query.py

Lines changed: 7 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -33,7 +33,7 @@ def complete(self) -> bool:
3333
if not super().complete():
3434
return False
3535

36-
datasets_to_check = [self.crdq.dataset]
36+
datasets_to_check = [self.crdq.dataset(self.dataset_type)]
3737
crdq_globals = Globals.from_ht(
3838
hl.read_table(self.output().path),
3939
datasets_to_check,
@@ -67,7 +67,11 @@ def requires(self) -> luigi.Task:
6767
)
6868
if self.crdq.query_raw_dataset:
6969
return HailTableTask(
70-
get_ht_path(CONFIG[self.crdq.dataset][self.reference_genome.v02_value]),
70+
get_ht_path(
71+
CONFIG[self.crdq.dataset(self.dataset_type)][
72+
self.reference_genome.v02_value
73+
],
74+
),
7175
)
7276
return HailTableTask(
7377
valid_reference_dataset_collection_path(
@@ -78,7 +82,7 @@ def requires(self) -> luigi.Task:
7882
)
7983

8084
def create_table(self) -> hl.Table:
81-
dataset: str = self.crdq.dataset
85+
dataset: str = self.crdq.dataset(self.dataset_type)
8286
if self.crdq.query_raw_dataset:
8387
query_ht = import_ht_from_config_path(
8488
CONFIG[dataset][self.reference_genome.v02_value],

0 commit comments

Comments
 (0)