Skip to content

Commit

Permalink
Updating the output structure for rules prep_metaquast and metaquast
Browse files Browse the repository at this point in the history
  • Loading branch information
skchronicles committed Mar 22, 2022
1 parent af75f16 commit 5044819
Show file tree
Hide file tree
Showing 2 changed files with 9 additions and 8 deletions.
4 changes: 2 additions & 2 deletions workflow/Snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -116,10 +116,10 @@ rule all:
),
# Metaquast, quality report of assembly
expand(
join(workpath,"temp","{name}","{name}.metaquast.fa"),
join(workpath,"{name}","temp","{name}.metaquast.fa"),
name=samples
),
expand(
join(workpath,"output","{name}","{name}_metaquast","report.html"),
join(workpath,"{name}","output","{name}_metaquast","report.html"),
name=samples
)
13 changes: 7 additions & 6 deletions workflow/rules/paired-end.smk
Original file line number Diff line number Diff line change
Expand Up @@ -528,12 +528,12 @@ rule prep_metaquast:
input:
report=join(workpath,"{name}","info","{name}.reads_kraken2_report.txt"),
output:
txt=join(workpath,"temp","{name}","{name}_reads_class_names.txt"),
fa=join(workpath,"temp","{name}","{name}.metaquast.fa"),
txt=join(workpath,"{name}","temp","{name}_reads_class_names.txt"),
fa=join(workpath,"{name}","temp","{name}.metaquast.fa"),
params:
rname='prepmetaq',
ncbi_viral=config['references']['ncbi_viral_fasta'],
outdir=join(workpath,"temp","{name}","metaquastref"),
outdir=join(workpath,"{name}","temp","metaquastref"),
threads: int(allocated("threads", "prep_metaquast", cluster))
envmodules:
config['tools']['ucsc']
Expand Down Expand Up @@ -570,12 +570,13 @@ rule metaquast:
input:
metaspades=join(workpath,"{name}","output","{name}.metaspades.contigs.fa"),
megahit=join(workpath,"{name}","output","{name}.megahit.contigs.fa"),
dep=join(workpath,"{name}","temp","{name}.metaquast.fa"),
output:
report=join(workpath,"output","{name}","{name}_metaquast","report.html")
report=join(workpath,"{name}","output","{name}_metaquast","report.html")
params:
rname='metaq',
ref=join(workpath,"temp","{name}","metaquastref"),
outdir=join(workpath,"output","{name}","{name}_metaquast"),
ref=join(workpath,"{name}","temp","metaquastref"),
outdir=join(workpath,"{name}","output","{name}_metaquast"),
threads: int(allocated("threads", "metaquast", cluster))
container:
config['images']['metaquast']
Expand Down

0 comments on commit 5044819

Please sign in to comment.