diff --git a/annotation.smk b/annotation.smk index 44290d94a437f300573d397a1f28721a97d870e2..8c72abb383877834c044eed9fbb2d7c9f2c200d4 100644 --- a/annotation.smk +++ b/annotation.smk @@ -63,6 +63,8 @@ def gff_input(wildcards): def db_input(wildcards): if wildcards.db == "nr": return {"db": "/db/outils/diamond-0.9.24/nr_tax.dmnd"} + elif wildcards.db == "swissprot": + return {"db": "/db/outils/diamond-0.9.26/uniprot_swissprot.dmnd"} elif wildcards.db == "phrogs": return {"db": "/projet/maiage/save/cmidoux/PhROGs/phrogs.dmnd"} diff --git a/global.smk b/global.smk index f9dc334bb0a4eed8e350b0888efa14958c8d843c..d76dd76209c27c74363ded137b102a395a3beebd 100644 --- a/global.smk +++ b/global.smk @@ -13,6 +13,7 @@ rule all: expand("report/count-contigs-coassembly-{sample}.tsv", sample=config["SAMPLES"]), expand("report/count-genes-coassembly-{sample}.tsv", sample=config["SAMPLES"]), "report/diamond_nr_coassembly.tsv", + "report/diamond_swissprot_coassembly.tsv", "report/contigs_coassembly-taxNames.tsv", "report/genes_coassembly-taxNames.tsv", #catalogue