Skip to content
Snippets Groups Projects
Commit 87d0a723 authored by Rachel  LEGENDRE's avatar Rachel LEGENDRE
Browse files

get good log output for multiQC

parent 5c425553
No related branches found
No related tags found
No related merge requests found
...@@ -29,7 +29,7 @@ rule sortmerna: ...@@ -29,7 +29,7 @@ rule sortmerna:
fasta = sortmerna_fasta fasta = sortmerna_fasta
output: output:
rRNA = temp(sortmerna_outfile_rRNA), rRNA = temp(sortmerna_outfile_rRNA),
no_rRNA = sortmerna_outfile_no_rRNA no_rRNA = temp(sortmerna_outfile_no_rRNA)
params: params:
tmpdir = config['tmpdir'] tmpdir = config['tmpdir']
singularity: singularity:
...@@ -54,7 +54,6 @@ rule sortmerna: ...@@ -54,7 +54,6 @@ rule sortmerna:
tmp="{output.no_rRNA}" tmp="{output.no_rRNA}"
outnorRNA=($tmp) outnorRNA=($tmp)
fasta="{input.fasta}" fasta="{input.fasta}"
idxdir=$(dirname ${{fasta}})"/sortmerna/" idxdir=$(dirname ${{fasta}})"/sortmerna/"
...@@ -66,6 +65,8 @@ rule sortmerna: ...@@ -66,6 +65,8 @@ rule sortmerna:
ext="" ext=""
fi fi
if [ -d {params.tmpdir}/{wildcards.SAMPLE} ]; then rm -rf {params.tmpdir}/{wildcards.SAMPLE} ; fi
if [[ ${{#infiles[@]}} -eq 2 ]]; if [[ ${{#infiles[@]}} -eq 2 ]];
then then
sortmerna --ref ${{fasta}} --idx-dir ${{idxdir}} --workdir {params.tmpdir}/{wildcards.SAMPLE} -threads {threads} --reads ${{infiles[0]}} --reads ${{infiles[1]}} --aligned out_paired_ribo --fastx --out2 --paired_out --other out_paired -v > {log.out} 2> {log.err} sortmerna --ref ${{fasta}} --idx-dir ${{idxdir}} --workdir {params.tmpdir}/{wildcards.SAMPLE} -threads {threads} --reads ${{infiles[0]}} --reads ${{infiles[1]}} --aligned out_paired_ribo --fastx --out2 --paired_out --other out_paired -v > {log.out} 2> {log.err}
...@@ -73,14 +74,17 @@ rule sortmerna: ...@@ -73,14 +74,17 @@ rule sortmerna:
mv out_paired_rev.fq${{ext}} {{outnorRNA[1]}} mv out_paired_rev.fq${{ext}} {{outnorRNA[1]}}
mv out_paired_ribo_fwd.fq${{ext}} {{outrRNA[0]}} mv out_paired_ribo_fwd.fq${{ext}} {{outrRNA[0]}}
mv out_paired_ribo_rev.fq${{ext}} {{outrRNA[1]}} mv out_paired_ribo_rev.fq${{ext}} {{outrRNA[1]}}
mv output.rRNA.log {log.out}
else else
sortmerna --ref ${{fasta}} --idx-dir ${{idxdir}} --workdir {params.tmpdir}/{wildcards.SAMPLE} -threads {threads} --reads {input.fastq} --aligned output.rRNA --fastx --other output.no_rRNA -v > {log.out} 2> {log.err} sortmerna --ref ${{fasta}} --idx-dir ${{idxdir}} --workdir {params.tmpdir}/{wildcards.SAMPLE} -threads {threads} --reads {input.fastq} --aligned output.rRNA --fastx --other output.no_rRNA -v > {log.out} 2> {log.err}
mv output.no_rRNA.fq${{ext}} {output.no_rRNA} mv output.no_rRNA.fq${{ext}} {output.no_rRNA}
mv output.rRNA.fq${{ext}} {output.rRNA} mv output.rRNA.fq${{ext}} {output.rRNA}
mv output.rRNA.log {log.out}
fi fi
rm -rf {params.tmpdir}/{wildcards.SAMPLE} rm -rf {params.tmpdir}/{wildcards.SAMPLE}
""" """
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment