diff --git a/tools/eggnogmapper2/diamond/Snakefile b/tools/eggnogmapper2/diamond/Snakefile new file mode 100644 index 0000000000000000000000000000000000000000..a387a6d30e40adf344b90fa5ede33427707e0c6c --- /dev/null +++ b/tools/eggnogmapper2/diamond/Snakefile @@ -0,0 +1,22 @@ +__eggnogmapper2_exec_command = config.get('eggnogmapper2', {}).get('exec_command', 'emapper.py') +__eggnogmapper2_diamond_options = config.get('eggnogmapper2', {}).get('diamond', {}).get('options', '') +__eggnogmapper2_diamond_threads = config.get('eggnogmapper2', {}).get('diamond', {}).get('threads', 4) + +rule eggnogmapper2_diamond: + """ + Run diamond mode of eggnogmapper v2 + """ + input: + __eggnogmapper2_diamond_input + output: + __eggnogmapper2_diamond_output + params: + exec_command = __eggnogmapper2_exec_command + ' -m diamond', + output_prefix = __eggnogmapper2_diamond_output_prefix, + options = __eggnogmapper2_diamond_options + threads: + __eggnogmapper2_diamond_threads + shell: + """ + {params.exec_command} {params.options} --cpu {threads} -i {input} -o {params.output_prefix} + """ diff --git a/tools/eggnogmapper2/diamond/example_usage/Snakefile b/tools/eggnogmapper2/diamond/example_usage/Snakefile new file mode 100644 index 0000000000000000000000000000000000000000..ef94ddfdd73ee745c20fa9a4405198280bdaf119 --- /dev/null +++ b/tools/eggnogmapper2/diamond/example_usage/Snakefile @@ -0,0 +1,22 @@ +configfile: "config.yaml" + +# ==== Snakefile path ==== +__eggnogmapper2_diamond_rules = config.get("snakefiles", {}).get("eggnogmapper2_diamond") + +__main_output_dir = config.get('output_dir', 'output') + +# ==== Main config ==== +SAMPLES = config.get('samples') +__input_dir = config.get('input_dir', 'data') + +# ==== EggNOGmapper2 Diamond ==== +__eggnogmapper2_output_dir = __main_output_dir + "/eggnogmapper2" +__eggnogmapper2_diamond_output_dir = __eggnogmapper2_output_dir + "/diamond" + +__eggnogmapper2_diamond_input = "{dir}/{{sample}}.fa".format(dir=__input_dir) +__eggnogmapper2_diamond_output_prefix = "{dir}/{{sample}}".format(dir=__eggnogmapper2_diamond_output_dir) +__eggnogmapper2_diamond_output = "{dir}/{{sample}}.emapper.seed_orthologs".format(dir=__eggnogmapper2_diamond_output_dir) +include: __eggnogmapper2_diamond_rules + +rule all: + input: expand("{dir}/{{sample}}.emapper.seed_orthologs".format(dir=__eggnogmapper2_diamond_output_dir), sample=SAMPLES) diff --git a/tools/eggnogmapper2/diamond/example_usage/config.yaml b/tools/eggnogmapper2/diamond/example_usage/config.yaml new file mode 100644 index 0000000000000000000000000000000000000000..9e007791b3a169384c84937698b223922fc9d66f --- /dev/null +++ b/tools/eggnogmapper2/diamond/example_usage/config.yaml @@ -0,0 +1,15 @@ +snakefiles: + eggnogmapper2_diamond: /pasteur/projets/policy01/Atm/snakemake/tools/eggnogmapper2/diamond/Snakefile + +input_dir: /pasteur/projets/policy01/Atm/kenzo/sandbox/20200210_test_snakemake/test_output/split_fasta +output_dir: /pasteur/projets/policy01/Atm/kenzo/sandbox/20200210_test_snakemake/test_output + +samples: +- test_00000 +- test_00001 +- test_00002 + +eggnogmapper2: + exec_command: "/pasteur/homes/kehillio/venv/eggnog-mapper-v2/bin/python /pasteur/homes/kehillio/tools/eggnog-mapper/emapper.py" + diamond: + options: "--no_annot --no_file_comments"