From 6affeb7aca71cd112803503e85b5218ad4295fbc Mon Sep 17 00:00:00 2001 From: Silas Kieser Date: Thu, 20 Jul 2023 11:16:11 +0200 Subject: [PATCH] fix bug for contig stats final --- workflow/rules/assemble.smk | 39 ++++++++++++++++++++----------------- 1 file changed, 21 insertions(+), 18 deletions(-) diff --git a/workflow/rules/assemble.smk b/workflow/rules/assemble.smk index 8dac96f5..f494bdcd 100644 --- a/workflow/rules/assemble.smk +++ b/workflow/rules/assemble.smk @@ -492,24 +492,6 @@ rule rename_contigs: " minscaf={params.minlength} &> {log} " -rule calculate_contigs_stats: - input: - "{sample}/assembly/{sample}_contigs.fasta", - output: - "{sample}/assembly/contig_stats/final_contig_stats.txt", - conda: - "../envs/required_packages.yaml" - log: - "{sample}/logs/assembly/post_process/contig_stats_final.log", - threads: 1 - resources: - mem=1, - time=config["runtime"]["simplejob"], - shell: - "stats.sh in={input} format=3 out={output} &> {log}" - - - if config["filter_contigs"]: @@ -623,6 +605,27 @@ rule finalize_contigs: os.symlink(os.path.relpath(input[0], os.path.dirname(output[0])), output[0]) + + +rule calculate_contigs_stats: + input: + "{sample}/{sample}_contigs.fasta", + output: + "{sample}/assembly/contig_stats/final_contig_stats.txt", + conda: + "../envs/required_packages.yaml" + log: + "{sample}/logs/assembly/post_process/contig_stats_final.log", + threads: 1 + resources: + mem=1, + time=config["runtime"]["simplejob"], + shell: + "stats.sh in={input} format=3 out={output} &> {log}" + + + + # generalized rule so that reads from any "sample" can be aligned to contigs from "sample_contigs" rule align_reads_to_final_contigs: input: