kkonganti@11: process QUAST { kkonganti@11: tag "$meta.id" kkonganti@11: label "process_micro" kkonganti@11: kkonganti@11: module (params.enable_module ? "${params.swmodulepath}${params.fs}quast${params.fs}5.2.0" : null) kkonganti@11: conda (params.enable_conda ? "bioconda::quast=5.2.0 conda-forge::libgcc-ng" : null) kkonganti@11: container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ? kkonganti@11: 'https://depot.galaxyproject.org/singularity/quast:5.2.0--py39pl5321h2add14b_1' : kkonganti@16: 'quay.io/biocontainers/quast:5.2.0--py39pl5321h2add14b_1' }" kkonganti@11: kkonganti@11: input: kkonganti@11: tuple val(meta) , path(consensus) kkonganti@11: tuple val(meta2), path(fasta) kkonganti@11: tuple val(meta3), path(gff) kkonganti@11: kkonganti@11: output: kkonganti@11: tuple val(meta), path("${prefix}") , emit: results kkonganti@11: tuple val(meta), path("${prefix}.quastreport.tsv") , emit: tsv kkonganti@11: tuple val(meta), path("${prefix}_transcriptome.tsv"), emit: transcriptome, optional: true kkonganti@11: tuple val(meta), path("${prefix}_misassemblies.tsv"), emit: misassemblies, optional: true kkonganti@11: tuple val(meta), path("${prefix}_unaligned.tsv") , emit: unaligned , optional: true kkonganti@11: path "versions.yml" , emit: versions kkonganti@11: kkonganti@11: when: kkonganti@11: task.ext.when == null || task.ext.when kkonganti@11: kkonganti@11: script: kkonganti@11: def args = task.ext.args ?: '' kkonganti@11: prefix = task.ext.prefix ?: "${meta.id}" kkonganti@11: def reference = fasta ? "-r $fasta" : '' kkonganti@11: def features = gff ? "--features $gff" : '' kkonganti@11: """ kkonganti@11: quast.py \\ kkonganti@11: -l $prefix \\ kkonganti@11: --output-dir $prefix \\ kkonganti@11: $reference \\ kkonganti@11: $features \\ kkonganti@11: --threads $task.cpus \\ kkonganti@11: $args \\ kkonganti@11: ${consensus.join(' ')} kkonganti@11: kkonganti@11: ln -s ${prefix}/report.tsv ${prefix}.quastreport.tsv kkonganti@11: [ -f ${prefix}/contigs_reports/all_alignments_transcriptome.tsv ] && ln -s ${prefix}/contigs_reports/all_alignments_transcriptome.tsv ${prefix}_transcriptome.tsv kkonganti@11: [ -f ${prefix}/contigs_reports/misassemblies_report.tsv ] && ln -s ${prefix}/contigs_reports/misassemblies_report.tsv ${prefix}_misassemblies.tsv kkonganti@11: [ -f ${prefix}/contigs_reports/unaligned_report.tsv ] && ln -s ${prefix}/contigs_reports/unaligned_report.tsv ${prefix}_unaligned.tsv kkonganti@11: kkonganti@11: cat <<-END_VERSIONS > versions.yml kkonganti@11: "${task.process}": kkonganti@11: quast: \$(quast.py --version 2>&1 | sed 's/^.*QUAST v//; s/ .*\$//') kkonganti@11: bash: \$( bash --version 2>&1 | sed '1!d; s/^.*version //; s/ (.*\$//' ) kkonganti@11: END_VERSIONS kkonganti@11: kkonganti@11: zcmd="" kkonganti@11: zver="" kkonganti@11: kkonganti@11: if type pigz > /dev/null 2>&1; then kkonganti@11: zcmd="pigz" kkonganti@11: zver=\$( echo \$( \$zcmd --version 2>&1 ) | sed -e '1!d' | sed "s/\$zcmd //" ) kkonganti@11: elif type gzip > /dev/null 2>&1; then kkonganti@11: zcmd="gzip" kkonganti@11: kkonganti@11: if [ "${workflow.containerEngine}" != "null" ]; then kkonganti@11: zver=\$( echo \$( \$zcmd --help 2>&1 ) | sed -e '1!d; s/ (.*\$//' ) kkonganti@11: else kkonganti@11: zver=\$( echo \$( \$zcmd --version 2>&1 ) | sed "s/^.*(\$zcmd) //; s/\$zcmd //; s/ Copyright.*\$//" ) kkonganti@11: fi kkonganti@11: fi kkonganti@11: kkonganti@11: cat <<-END_VERSIONS >> versions.yml kkonganti@11: \$zcmd: \$zver kkonganti@11: END_VERSIONS kkonganti@11: """ kkonganti@11: }