Skip to content

Commit

Permalink
Remove output files from all_outputs list that do not exist anymore
Browse files Browse the repository at this point in the history
  • Loading branch information
verku committed Oct 9, 2023
1 parent 92f7aa9 commit 4c49098
Show file tree
Hide file tree
Showing 2 changed files with 0 additions and 8 deletions.
4 changes: 0 additions & 4 deletions workflow/rules/12_snpEff.smk
Original file line number Diff line number Diff line change
Expand Up @@ -4,23 +4,19 @@

# Code collecting output files from this part of the pipeline
if os.path.exists(config["historical_samples"]) and os.path.exists(config["modern_samples"]):
all_outputs.append(expand("results/{dataset}/vcf/" + REF_NAME + "/stats/multiqc/multiqc_report.html",
dataset=["historical", "modern"]))
all_outputs.append(expand("results/all/snpEff/" + REF_NAME + ".all.fmissing{fmiss}.{chr}.snpEff_variant_impact_plot.pdf",
fmiss=config["f_missing"],
chr=CHR,))
all_outputs.append("results/historical/snpEff/" + REF_NAME + "/multiqc/multiqc_report.html")
all_outputs.append("results/modern/snpEff/" + REF_NAME + "/multiqc/multiqc_report.html")

elif os.path.exists(config["historical_samples"]):
all_outputs.append("results/historical/vcf/" + REF_NAME + "/stats/multiqc/multiqc_report.html")
all_outputs.append(expand("results/historical/snpEff/" + REF_NAME + ".historical.fmissing{fmiss}.{chr}.snpEff_variant_impact_plot.pdf",
fmiss=config["f_missing"],
chr=CHR,))
all_outputs.append("results/historical/snpEff/" + REF_NAME + "/multiqc/multiqc_report.html")

elif os.path.exists(config["modern_samples"]):
all_outputs.append("results/modern/vcf/" + REF_NAME + "/stats/multiqc/multiqc_report.html")
all_outputs.append(expand("results/modern/snpEff/" + REF_NAME + ".modern.fmissing{fmiss}.{chr}.snpEff_variant_impact_plot.pdf",
fmiss=config["f_missing"],
chr=CHR,))
Expand Down
4 changes: 0 additions & 4 deletions workflow/rules/13_GERP.smk
Original file line number Diff line number Diff line change
Expand Up @@ -7,24 +7,20 @@ all_outputs.append(expand("results/gerp/" + REF_NAME + ".{chr}.ancestral.rates.g
chr=CHR,))

if os.path.exists(config["historical_samples"]) and os.path.exists(config["modern_samples"]):
all_outputs.append(expand("results/{dataset}/vcf/" + REF_NAME + "/stats/multiqc/multiqc_report.html",
dataset=["historical", "modern"]))
all_outputs.append(expand("results/gerp/all/" + REF_NAME + ".all.fmissing{fmiss}.{chr}.relative_mutational_load.gerp_{minGERP}_{maxGERP}_plot.pdf",
fmiss=config["f_missing"],
chr=CHR,
minGERP=config["min_gerp"],
maxGERP=config["max_gerp"],))

elif os.path.exists(config["historical_samples"]):
all_outputs.append("results/historical/vcf/" + REF_NAME + "/stats/multiqc/multiqc_report.html")
all_outputs.append(expand("results/gerp/historical/" + REF_NAME + ".historical.fmissing{fmiss}.{chr}.relative_mutational_load.gerp_{minGERP}_{maxGERP}_plot.pdf",
fmiss=config["f_missing"],
chr=CHR,
minGERP=config["min_gerp"],
maxGERP=config["max_gerp"],))

elif os.path.exists(config["modern_samples"]):
all_outputs.append("results/modern/vcf/" + REF_NAME + "/stats/multiqc/multiqc_report.html")
all_outputs.append(expand("results/gerp/modern/" + REF_NAME + ".modern.fmissing{fmiss}.{chr}.relative_mutational_load.gerp_{minGERP}_{maxGERP}_plot.pdf",
fmiss=config["f_missing"],
chr=CHR,
Expand Down

0 comments on commit 4c49098

Please sign in to comment.