Skip to content

Commit

Permalink
fix: get-rename-contig function
Browse files Browse the repository at this point in the history
  • Loading branch information
famosab committed Jan 24, 2025
1 parent bb10504 commit 77f354f
Showing 1 changed file with 16 additions and 6 deletions.
22 changes: 16 additions & 6 deletions workflow/rules/common.smk
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ used_callsets = {callset for callset in callsets.keys()}
used_genomes = {benchmarks[benchmark]["genome"] for benchmark in used_benchmarks}


#TODO: can this be removed?
# TODO: can this be removed?
if any(
callset["benchmark"] == "giab-NA12878-exome" for callset in callsets.values()
) and config.get("grch37"):
Expand Down Expand Up @@ -337,12 +337,22 @@ def get_test_regions(wildcards):


def get_rename_contig_file(wildcards):
if config["variant-calls"][wildcards.callset]["genome-build"] == "grch37" and config["variant-calls"][wildcards.callset]["rename-contigs"]:
return workflow.source_path("../resources/rename-contigs/grch37_ucsc2ensembl.txt")
if config["variant-calls"][wildcards.callset]["genome-build"] == "grch38" and config["variant-calls"][wildcards.callset]["rename-contigs"]:
return workflow.source_path("../resources/rename-contigs/grch38_ucsc2ensembl.txt")
if (
config["variant-calls"][wildcards.callset]["genome-build"] == "grch37"
and config["variant-calls"][wildcards.callset]["rename-contigs"]
):
return workflow.source_path(
"../resources/rename-contigs/grch37_ucsc2ensembl.txt"
)
if (
config["variant-calls"][wildcards.callset]["genome-build"] == "grch38"
and config["variant-calls"][wildcards.callset]["rename-contigs"]
):
return workflow.source_path(
"../resources/rename-contigs/grch38_ucsc2ensembl.txt"
)
else:
config["variant-calls"][wildcards.callset]["rename-contigs"]
return config["variant-calls"][wildcards.callset]["rename-contigs"]


def get_callset_subcategory(wildcards):
Expand Down

0 comments on commit 77f354f

Please # to comment.