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