@@ -602,37 +602,37 @@ process adapter_removal {
602602 script:
603603 base = reads[0 ]. baseName
604604 // This checks whether we skip trimming and defines a variable respectively
605- trim_me = params. skip_trim ? ' ' : ' '
605+ trim_me = params. skip_trim ? ' ' : " --trimns --trimqualities --adapter1 ${ params.clip_forward_adaptor} --adapter2 ${ params.clip_reverse_adaptor} --minlength ${ params.clip_readlength} --minquality ${ params.clip_min_read_quality} --minadapteroverlap ${ params.min_adap_overlap} "
606+ collapse_me = params. skip_collapse ? ' ' : ' --collapse'
606607
607- // PE, collapse & trim reads
608+ // PE mode, dependent on trim_me and collapse_me the respective procedure is run or not :-)
608609 if (! params. singleEnd && ! params. skip_collapse && ! params. skip_trim){
609610 """
610611 mkdir -p output
611- AdapterRemoval --file1 ${ reads[0]} --file2 ${ reads[1]} --basename ${ base} --gzip --threads ${ task.cpus} --trimns --trimqualities --adapter1 ${ params.clip_forward_adaptor } --adapter2 ${ params.clip_reverse_adaptor } --minlength ${ params.clip_readlength } --minquality ${ params.clip_min_read_quality } --minadapteroverlap ${ params.min_adap_overlap } --collapse
612+ AdapterRemoval --file1 ${ reads[0]} --file2 ${ reads[1]} --basename ${ base} ${ trim_me } --gzip --threads ${ task.cpus} ${ collapse_me }
612613 #Combine files
613614 zcat *.collapsed.gz *.collapsed.truncated.gz *.singleton.truncated.gz *.pair1.truncated.gz *.pair2.truncated.gz | gzip > output/${ base} .combined.fq.gz
614615 """
615616 // PE, don't collapse, but trim reads
616617 } else if (! params. singleEnd && params. skip_collapse && ! params. skip_trim) {
617618 """
618619 mkdir -p output
619- AdapterRemoval --file1 ${ reads[0]} --file2 ${ reads[1]} --basename ${ base} --gzip --threads ${ task.cpus} --trimns --trimqualities --adapter1 ${ params.clip_forward_adaptor } --adapter2 ${ params.clip_reverse_adaptor } --minlength ${ params.clip_readlength } --minquality ${ params.clip_min_read_quality } --minadapteroverlap ${ params.min_adap_overlap }
620+ AdapterRemoval --file1 ${ reads[0]} --file2 ${ reads[1]} --basename ${ base} --gzip --threads ${ task.cpus} ${ trim_me } ${ collapse_me }
620621 mv ${ base} .pair*.truncated.gz output/
621622 """
622623 // PE, collapse, but don't trim reads
623624 } else if (! params. singleEnd && ! params. skip_collapse && params. skip_trim) {
624- bogus_adaptor = " NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN"
625625 """
626626 mkdir -p output
627- AdapterRemoval --file1 ${ reads[0]} --file2 ${ reads[1]} --basename ${ base} --gzip --threads ${ task.cpus} --basename ${ base} --collapse --adapter1 $b ogus_adaptor --adapter2 $b ogus_adaptor
627+ AdapterRemoval --file1 ${ reads[0]} --file2 ${ reads[1]} --basename ${ base} --gzip --threads ${ task.cpus} --basename ${ base} ${ collapse_me } ${ trim_me }
628628
629629 mv ${ base} .pair*.truncated.gz output/
630630 """
631631 } else {
632632 // SE, collapse not possible, trim reads
633633 """
634634 mkdir -p output
635- AdapterRemoval --file1 ${ reads[0]} --basename ${ base} --gzip --threads ${ task.cpus} --trimns --trimqualities --adapter1 ${ params.clip_forward_adaptor } --minlength ${ params.clip_readlength } --minquality ${ params.clip_min_read_quality }
635+ AdapterRemoval --file1 ${ reads[0]} --basename ${ base} --gzip --threads ${ task.cpus} ${ trim_me }
636636
637637 mv *.truncated.gz output/
638638 """
0 commit comments