variant analysis config order in common.smk

sanjaynagi · Oct 17, 2024 · 38fad19 · 38fad19
1 parent 27587d0
commit 38fad19
Showing 1 changed file with 40 additions and 38 deletions.
diff --git a/workflow/rules/common.smk b/workflow/rules/common.smk
@@ -163,15 +163,6 @@ def GetDesiredOutputs(wildcards):
             )
         )
 
-        if config['QualityControl']['coverage']['activate']:
-            wanted_input.extend(
-                expand(
-                    [
-                        "results/qc/coverage/{sample}.mosdepth.summary.txt",
-                    ],
-                sample=samples,
-                )
-            )
 
         if config["VariantAnalysis"]["selection"]["activate"]:
             wanted_input.extend(
@@ -189,19 +180,53 @@ def GetDesiredOutputs(wildcards):
                 )
             )
 
-    if config["VariantAnalysis"]["ancestry"]["activate"]:
+        if config["VariantAnalysis"]["ancestry"]["activate"]:
+            wanted_input.extend(
+                expand(
+                    [
+                        "results/variantAnalysis/ancestry/AIMs_summary.tsv",
+                        "results/variantAnalysis/ancestry/AIM_fraction_whole_genome.svg",
+                        "results/variantAnalysis/ancestry/n_AIMS_per_chrom.tsv",
+                        "results/variantAnalysis/ancestry/AIM_fraction_{contig}.tsv",
+                    ],
+                    contig=config["contigs"],
+                )
+            )
+
+        if config["VariantAnalysis"]["karyotype"]["activate"]:
+            wanted_input.extend(
+                expand(
+                    [
+                        "results/karyotype/{karyo}.{dataset}.karyo.txt",
+                        "results/karyotype/karyoFreqs.svg",
+                    ],
+                    karyo=config["VariantAnalysis"]["karyotype"]["inversions"],
+                    dataset=config["dataset"],
+                )
+            )
+
+        if (
+        config["VariantAnalysis"]['selection']["activate"]
+        and config["DifferentialExpression"]["GSEA"]["activate"]
+        ):
+            wanted_input.extend(
+                expand(
+                    ["results/gsea/fst/{comp}.fst.tsv"],
+                    comp=config["contrasts"],
+                )
+            )
+
+    if config['QualityControl']['coverage']['activate']:
         wanted_input.extend(
             expand(
                 [
-                    "results/variantAnalysis/ancestry/AIMs_summary.tsv",
-                    "results/variantAnalysis/ancestry/AIM_fraction_whole_genome.svg",
-                    "results/variantAnalysis/ancestry/n_AIMS_per_chrom.tsv",
-                    "results/variantAnalysis/ancestry/AIM_fraction_{contig}.tsv",
+                    "results/qc/coverage/{sample}.mosdepth.summary.txt",
                 ],
-                contig=config["contigs"],
+            sample=samples,
             )
         )
 
+
     if config["miscellaneous"]["VariantsOfInterest"]["activate"]:
         wanted_input.extend(
             [
@@ -221,29 +246,6 @@ def GetDesiredOutputs(wildcards):
             )
         )
 
-    if (
-        config["VariantAnalysis"]['selection']["activate"]
-        and config["DifferentialExpression"]["GSEA"]["activate"]
-    ):
-        wanted_input.extend(
-            expand(
-                ["results/gsea/fst/{comp}.fst.tsv"],
-                comp=config["contrasts"],
-            )
-        )
-
-    if config["VariantAnalysis"]["karyotype"]["activate"]:
-        wanted_input.extend(
-            expand(
-                [
-                    "results/karyotype/{karyo}.{dataset}.karyo.txt",
-                    "results/karyotype/karyoFreqs.svg",
-                ],
-                karyo=config["VariantAnalysis"]["karyotype"]["inversions"],
-                dataset=config["dataset"],
-            )
-        )
-
     if config["miscellaneous"]["sweeps"]["activate"]:
         if config["DifferentialExpression"]['gene-level']["activate"]:
             wanted_input.extend(