diff --git a/Snakefile b/Snakefile index 7ceb77189ddf8f4cb7c559f13c6cf089b8556ad3..254f1b9cb60b64107af19b3406863d0e00e9ede1 100755 --- a/Snakefile +++ b/Snakefile @@ -726,7 +726,7 @@ if config["seacr"]["do"]: seacr_output = os.path.join(analysis_dir, "06-PeakCalling/seacr/{{SAMPLE}}.{}.bed".format(config["seacr"]["threshold"])) seacr_logs_std = os.path.join(analysis_dir, "06-PeakCalling/seacr/logs/{SAMPLE}_seacr_calling.out") seacr_logs_err = os.path.join(analysis_dir, "06-PeakCalling/seacr/logs/{SAMPLE}_seacr_calling.err") - final_output.extend(expand(seacr_output, SAMPLE=ALL_IP_PC)) + final_output.extend(expand(seacr_output, SAMPLE=IP_ALL)) include: os.path.join(RULES, "seacr.rules") @@ -739,7 +739,7 @@ if config["macs2"]["do"] or config["seacr"]["do"] : if wildcards.CALLER == "macs2": return expand(os.path.join(analysis_dir, "06-PeakCalling/{{CALLER}}/%s/{IP_REP}_peaks.{{MOD}}Peak" % (model_dir)), IP_REP=ALL_IP_PC) elif wildcards.CALLER == "seacr": - return expand(os.path.join(analysis_dir, "06-PeakCalling/{{CALLER}}/{IP_REP}.{{MOD}}.bed"), IP_REP=ALL_IP_PC) + return expand(os.path.join(analysis_dir, "06-PeakCalling/{{CALLER}}/{IP_REP}.{{MOD}}.bed"), IP_REP=IP_ALL) stats_peakCalling_input = stats_pc_input stats_peakCalling_csv = os.path.join(analysis_dir, "{CALLER}_{MOD}_Peaks_metrics_mqc.out")