kkonganti@105: process CAT_FASTQ {
kkonganti@105:     tag "$meta.id"
kkonganti@105:     label 'process_micro'
kkonganti@105: 
kkonganti@105:     conda (params.enable_conda ? "conda-forge::sed=4.7 conda-forge::gzip" : null)
kkonganti@105:     container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
kkonganti@105:         'https://containers.biocontainers.pro/s3/SingImgsRepo/biocontainers/v1.2.0_cv1/biocontainers_v1.2.0_cv1.img' :
kkonganti@105:         'biocontainers/biocontainers:v1.2.0_cv1' }"
kkonganti@105: 
kkonganti@105:     input:
kkonganti@105:         tuple val(meta), path(reads, stageAs: "input*/*")
kkonganti@105: 
kkonganti@105:     output:
kkonganti@105:         tuple val(meta), path("*.merged.fastq.gz"), emit: catted_reads
kkonganti@105:         path "versions.yml"                       , emit: versions
kkonganti@105: 
kkonganti@105:     when:
kkonganti@105:         task.ext.when == null || task.ext.when
kkonganti@105: 
kkonganti@105:     script:
kkonganti@105:         def args = task.ext.args ?: ''
kkonganti@105:         def prefix = task.ext.prefix ?: "${meta.id}"
kkonganti@105:         def readList = reads.collect{ it.toString() }
kkonganti@105:         def is_in_gz = readList[0].endsWith('.gz')
kkonganti@105:         def gz_or_ungz = (is_in_gz ? '' : ' | gzip')
kkonganti@105:         def pigz_or_ungz = (is_in_gz ? '' : " | pigz -p ${task.cpus}")
kkonganti@105:         if (meta.single_end) {
kkonganti@105:             if (readList.size > 1) {
kkonganti@105:                 """
kkonganti@105:                 zcmd="gzip"
kkonganti@105:                 zver=""
kkonganti@105: 
kkonganti@105:                 if type pigz > /dev/null 2>&1; then
kkonganti@105:                     cat ${readList.join(' ')} ${pigz_or_ungz} > ${prefix}.merged.fastq.gz
kkonganti@105:                     zcmd="pigz"
kkonganti@105:                     zver=\$( echo \$( \$zcmd --version 2>&1 ) | sed -e '1!d' | sed "s/\$zcmd //" )
kkonganti@105:                 else
kkonganti@105:                     cat ${readList.join(' ')} ${gz_or_ungz} > ${prefix}.merged.fastq.gz
kkonganti@105:                     zcmd="gzip"
kkonganti@105: 
kkonganti@105:                     if [ "${workflow.containerEngine}" != "null" ]; then
kkonganti@105:                         zver=\$( echo \$( \$zcmd --help 2>&1 ) | sed -e '1!d; s/ (.*\$//' )
kkonganti@105:                     else
kkonganti@105:                         zver=\$( echo \$( \$zcmd --version 2>&1 ) | sed "s/^.*(\$zcmd) //; s/\$zcmd //; s/ Copyright.*\$//" )
kkonganti@105:                     fi
kkonganti@105:                 fi
kkonganti@105: 
kkonganti@105:                 cat <<-END_VERSIONS > versions.yml
kkonganti@105:                 "${task.process}":
kkonganti@105:                     cat: \$( echo \$(cat --version 2>&1) | sed 's/^.*coreutils) //; s/ .*\$//' )
kkonganti@105:                     \$zcmd: \$zver
kkonganti@105:                 END_VERSIONS
kkonganti@105:                 """
kkonganti@105:             }
kkonganti@105:         } else {
kkonganti@105:             if (readList.size > 2) {
kkonganti@105:                 def read1 = []
kkonganti@105:                 def read2 = []
kkonganti@105:                 readList.eachWithIndex{ v, ix -> ( ix & 1 ? read2 : read1 ) << v }
kkonganti@105:                 """
kkonganti@105:                 zcmd="gzip"
kkonganti@105:                 zver=""
kkonganti@105: 
kkonganti@105:                 if type pigz > /dev/null 2>&1; then
kkonganti@105:                     cat ${read1.join(' ')} ${pigz_or_ungz} > ${prefix}_1.merged.fastq.gz
kkonganti@105:                     cat ${read2.join(' ')} ${pigz_or_ungz} > ${prefix}_2.merged.fastq.gz
kkonganti@105:                     zcmd="pigz"
kkonganti@105:                     zver=\$( echo \$( \$zcmd --version 2>&1 ) | sed -e '1!d' | sed "s/\$zcmd //" )
kkonganti@105:                 else
kkonganti@105:                     cat ${read1.join(' ')} ${gz_or_ungz} > ${prefix}_1.merged.fastq.gz
kkonganti@105:                     cat ${read2.join(' ')} ${gz_or_ungz} > ${prefix}_2.merged.fastq.gz
kkonganti@105:                     zcmd="gzip"
kkonganti@105: 
kkonganti@105:                     if [ "${workflow.containerEngine}" != "null" ]; then
kkonganti@105:                         zver=\$( echo \$( \$zcmd --help 2>&1 ) | sed -e '1!d; s/ (.*\$//' )
kkonganti@105:                     else
kkonganti@105:                         zver=\$( echo \$( \$zcmd --version 2>&1 ) | sed "s/^.*(\$zcmd) //; s/\$zcmd //; s/ Copyright.*\$//" )
kkonganti@105:                     fi
kkonganti@105:                 fi
kkonganti@105: 
kkonganti@105:                 cat <<-END_VERSIONS > versions.yml
kkonganti@105:                 "${task.process}":
kkonganti@105:                     cat: \$( echo \$(cat --version 2>&1) | sed 's/^.*coreutils) //; s/ .*\$//' )
kkonganti@105:                     \$zcmd: \$zver
kkonganti@105:                 END_VERSIONS
kkonganti@105:                 """
kkonganti@105:             }
kkonganti@105:         }
kkonganti@105: }