Skip to content
Snippets Groups Projects
Commit 06c9e00f authored by seehagec01's avatar seehagec01
Browse files

kraken kinda working

parent 8047f854
No related branches found
No related tags found
No related merge requests found
...@@ -2,30 +2,18 @@ rule build: ...@@ -2,30 +2,18 @@ rule build:
input: input:
"results/fastp/trimmed/ERR4082860.1.fastq" "results/fastp/trimmed/ERR4082860.1.fastq"
output: output:
"resources/db/yourdb" "results/krakendb/db"
conda:
"../envs/yourenv.yaml"
shell:
"kraken-build --standard --db resources/krakendb/db"
rule classification:
input:
"results/fastp/trimmed/ERR4082860.1.fastq",
"results/fastp/trimmed/ERR4082860.2.fastq"
output:
"results/kraken/ERR4082860_kraken_cseqs#.fq"
threads:4 threads:4
conda: conda:
"../envs/yourenv.yaml" "../envs/yourenv.yaml"
shell: shell:
"kraken2 --gzip-compressed --paired --classified-out {output} {input[0]} {input[1]} --db resources/db/yourdb --report ERR4082860_report.txt --threads {threads}" "kraken2-build --standard --threads {threads} --db db"
rule report: rule report:
input: input:
db="resources/db/yourdb", db="results/krakendb/db",
sample=["results/fastp/trimmed/{sample}_1.fastq", "results/fastp/trimmed/{sample}_2.fastq"] sample=["results/fastp/trimmed/{sample}.1.fastq", "results/fastp/trimmed/{sample}.2.fastq"]
output: output:
"results/kraken/{sample}_report.txt" "results/kraken/{sample}_report.txt"
conda: conda:
...@@ -33,4 +21,4 @@ rule report: ...@@ -33,4 +21,4 @@ rule report:
log: log:
"workflow/report/kraken/{sample}.log" "workflow/report/kraken/{sample}.log"
shell: shell:
"kraken-report --db {input.db} {input.sample} > {output} 2> {log}" "kraken2 --report --db {input.db} {input.sample} > {output} 2> {log}"
\ No newline at end of file \ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment