diff --git a/assembly.smk b/assembly.smk index d290c20d4fd38110a139d99c864b0dd2852b09c3..bb03a030997eeeee404c842c08573e8eea2ea912 100644 --- a/assembly.smk +++ b/assembly.smk @@ -93,3 +93,9 @@ rule metaspades: "-i work/metaSPADES/{wildcards.reads}/contigs.fasta " "-o {output.contigs} " "--min-length {params.min_len} " + +def contigs_input(wildcards): + if config["ASSEMBLER"] == "megahit": + return ["work/megahit/{wildcards.sample}/{wildcards.sample}.contigs.fa".format(wildcards=wildcards)] + elif config["ASSEMBLER"] == "metaspades": + return ["work/metaSPADES/{wildcards.sample}/{wildcards.sample}.contigs.fasta".format(wildcards=wildcards)] diff --git a/global.smk b/global.smk index dc8d4ba750556baa1201e96fbf616ca55ccb3a20..460dcdc3597e7f0d1f2cb5e3aac2f9c2099c22e0 100644 --- a/global.smk +++ b/global.smk @@ -8,6 +8,7 @@ rule all: input: "report/multiqc_report.html", expand("report/reads_{sample}-krona.html", sample=config["SAMPLES"]), + "report/contigs_coassembly-taxNames.tsv", include: "../workflow_metagenomics/quality.smk" include: "../workflow_metagenomics/preprocess.smk"