File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -601,24 +601,26 @@ process adapter_removal {
601601
602602 script:
603603 base = reads[0 ]. baseName
604+ // This checks whether we skip trimming and defines a variable respectively
605+ trim_me = params. skip_trim ? ' ' : ' '
604606
605607 // PE, collapse & trim reads
606- if (! params. singleEnd && ! params. skip_collapse && ! params. skip_trim && ! params . skip_adapterremoval ){
608+ if (! params. singleEnd && ! params. skip_collapse && ! params. skip_trim){
607609 """
608610 mkdir -p output
609611 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
610612 #Combine files
611613 zcat *.collapsed.gz *.collapsed.truncated.gz *.singleton.truncated.gz *.pair1.truncated.gz *.pair2.truncated.gz | gzip > output/${ base} .combined.fq.gz
612614 """
613615 // PE, don't collapse, but trim reads
614- } else if (! params. singleEnd && params. skip_collapse && ! params. skip_trim && ! params . skip_adapterremoval ) {
616+ } else if (! params. singleEnd && params. skip_collapse && ! params. skip_trim) {
615617 """
616618 mkdir -p output
617619 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}
618620 mv ${ base} .pair*.truncated.gz output/
619621 """
620622 // PE, collapse, but don't trim reads
621- } else if (! params. singleEnd && ! params. skip_collapse && params. skip_trim && ! params . skip_adapterremoval ) {
623+ } else if (! params. singleEnd && ! params. skip_collapse && params. skip_trim) {
622624 bogus_adaptor = " NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN"
623625 """
624626 mkdir -p output
You can’t perform that action at this time.
0 commit comments