Skip to content

Commit

Permalink
allows to parametrize import qc
Browse files Browse the repository at this point in the history
importqc_params: 
"iupacToN=t touppercase=t qout=33 addslash=t trimreaddescription=t"
  • Loading branch information
SilasK committed Oct 27, 2020
1 parent a744644 commit 17e0784
Show file tree
Hide file tree
Showing 3 changed files with 5 additions and 7 deletions.
2 changes: 1 addition & 1 deletion atlas/__init__.py
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
__version__ = "2.4.3"
__version__ = "2.4.4"

from .scripts import utils

Expand Down
1 change: 1 addition & 0 deletions atlas/default_values.py
Original file line number Diff line number Diff line change
Expand Up @@ -129,6 +129,7 @@ def make_default_config():

conf["simplejob_mem"]=10
conf["simplejob_threads"]=4
conf["importqc_params"]= "iupacToN=t touppercase=t qout=33 addslash=t trimreaddescription=t"

conf["mem"] = MEM
conf["large_mem"] = 100
Expand Down
9 changes: 3 additions & 6 deletions atlas/rules/qc.snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,8 @@ rule initialize_qc:
inputs = lambda wc, input: io_params_for_tadpole(input,'in'),
interleaved = lambda wc: "t" if config.get('interleaved_fastqs',False) else "f",
outputs = lambda wc, output: io_params_for_tadpole(output,'out'),
verifypaired = "t" if PAIRED_END else "f"
verifypaired = "t" if PAIRED_END else "f",
extra = config['importqc_params']
log:
"{sample}/logs/QC/init.log"
conda:
Expand All @@ -89,13 +90,9 @@ rule initialize_qc:
reformat.sh {params.inputs} \
interleaved={params.interleaved} \
{params.outputs} \
iupacToN=t \
touppercase=t \
qout=33 \
{params.extra} \
overwrite=true \
verifypaired={params.verifypaired} \
addslash=t \
trimreaddescription=t \
threads={threads} \
-Xmx{resources.java_mem}G 2> {log}
"""
Expand Down

0 comments on commit 17e0784

Please sign in to comment.