swissprot

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"}
 

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