modules.config

/*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    Config file for defining DSL2 per module options and publishing paths
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    Available keys to override module options:
        ext.args   = Additional arguments appended to command in module.
        ext.args2  = Second set of arguments appended to command in module (multi-tool modules).
        ext.args3  = Third set of arguments appended to command in module (multi-tool modules).
        ext.prefix = File name prefix for output files.
----------------------------------------------------------------------------------------
*/

process {

    publishDir = [
        path: { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" },
        mode: params.publish_dir_mode,
        saveAs: { filename -> filename.equals('versions.yml') ? null : filename }
    ]

    withName: 'MULTIQC' {
        ext.args   = { params.multiqc_title ? "--title \"$params.multiqc_title\"" : '' }
        publishDir = [
            path: { "${params.outdir}/multiqc" },
            mode: params.publish_dir_mode,
            saveAs: { filename -> filename.equals('versions.yml') ? null : filename }
        ]
    }

    //
    // CONVERT INPUT BAM
    //
    withName: SAMTOOLS_CONVERT_BAM_INPUT {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: CAT_FASTQ_CONVERTED_BAM {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            enabled: false
        ]
    }

    //
    // READ PREPROCESSING
    //
    withName: FASTQC {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = '--quiet'
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/preprocessing/fastqc_raw/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: FASTQC_PROCESSED {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = '--quiet'
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/preprocessing/fastqc_processed/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: FALCO {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = '--quiet'
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/preprocessing/falco_raw/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: FALCO_PROCESSED {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = '--quiet'
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/preprocessing/falco_processed/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: FASTP_SINGLE {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args   = [
            // TRIMMING OPTIONS
            params.preprocessing_trim5p != 0 ? "--trim_front1 ${params.preprocessing_trim5p}" : "",
            params.preprocessing_trim3p != 0 ? "--trim_tail1 ${params.preprocessing_trim3p}" : "",
            params.preprocessing_skipadaptertrim ? "--disable_adapter_trimming" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter1 ? "--adapter_sequence ${params.preprocessing_adapter1}" : "",  // adding adapter list happens at module input channel level
            // FILTERING OPTIONS
            "--length_required ${params.preprocessing_minlength}",
            params.preprocessing_fastp_complexityfilter ? "--low_complexity_filter --complexity_threshold ${params.preprocessing_fastp_complexityfilter_threshold}" : ""
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/preprocessing/fastp/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.preprocessing_savepreprocessedreads
            ],
            [
                //stats
                path: { "${params.outdir}/preprocessing/fastp/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.{log,html,json}'
            ]
        ]
    }

    withName: FASTP_PAIRED {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args   = [
            // COLLAPSING OPTIONS - option to retain singletons
            params.preprocessing_excludeunmerged ? "" : "--include_unmerged",
            // TRIMMING OPTIONS
            params.preprocessing_trim5p != 0 ? "--trim_front1 ${params.preprocessing_trim5p} --trim_front2 ${params.preprocessing_trim5p}" : "",
            params.preprocessing_trim3p != 0 ? "--trim_tail1 ${params.preprocessing_trim3p} --trim_tail2 ${params.preprocessing_trim3p}" : "",
            params.preprocessing_skipadaptertrim ? "--disable_adapter_trimming" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter1 ? "--adapter_sequence ${params.preprocessing_adapter1}" : "",  // adding adapter list happens at module input channel level
            params.preprocessing_skipadaptertrim ? "" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter2 ? "--adapter_sequence_r2 ${params.preprocessing_adapter2}" : "",  // adding adapter list happens at module input channel level
            // FILTERING OPTIONS--disable_adapter_trimming
            "--length_required ${params.preprocessing_minlength}",
            params.preprocessing_fastp_complexityfilter ? "--low_complexity_filter --complexity_threshold ${params.preprocessing_fastp_complexityfilter_threshold}" : ""
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/preprocessing/fastp/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.preprocessing_savepreprocessedreads
            ],
            [
                // stats
                path: { "${params.outdir}/preprocessing/fastp/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.{log,html,json}'
            ]
        ]
    }

    withName: ADAPTERREMOVAL_SINGLE {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args   = [
            // trimming options - note: adding adapter list happens at module input channel level
            params.preprocessing_skipadaptertrim ? "--adapter1 ''" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter1 ? "--adapter1 ${params.preprocessing_adapter1}" : "", // adding adapter list happens at module input channel level
            params.preprocessing_skipadaptertrim ? "--adapter2 ''" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter2 ? "--adapter2 ${params.preprocessing_adapter2}" : "",
            "--minadapteroverlap ${params.preprocessing_adapterremoval_adapteroverlap}",
            params.preprocessing_adapterremoval_preserve5p ? "--preserve5p" : "",
            params.preprocessing_trim5p != 0 ? "--trim5p ${params.preprocessing_trim5p}" : "",
            params.preprocessing_trim3p != 0 ? "--trim3p ${params.preprocessing_trim3p}" : "",
            !params.preprocessing_adapterremoval_skipqualitytrimming ? "--trimqualities --minquality ${params.preprocessing_adapterremoval_trimbasequalitymin}" : "",
            !params.preprocessing_adapterremoval_skipntrimming ? "--trimns" : "",
            "--qualitymax ${params.preprocessing_adapterremoval_qualitymax}",
            // filtering options
            "--minlength ${params.preprocessing_minlength}"
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/preprocessing/adapterremoval/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.preprocessing_savepreprocessedreads
            ],
            [
                // stats
                path: { "${params.outdir}/preprocessing/adapterremoval/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.settings'
            ]
        ]
    }

    withName: ADAPTERREMOVAL_PAIRED {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args   = [
            // COLLAPSING OPTIONS
            params.preprocessing_skippairmerging ? "" : "--collapse",
            // TRIMMING OPTIONS  - note: adding adapter list happens at module input channel level
            params.preprocessing_skipadaptertrim ? "--adapter1 ''" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter1 ? "--adapter1 ${params.preprocessing_adapter1}" : "",
            params.preprocessing_skipadaptertrim ? "--adapter2 ''" : params.preprocessing_adapterlist ? "" : params.preprocessing_adapter2 ? "--adapter2 ${params.preprocessing_adapter2}" : "",
            "--minadapteroverlap ${params.preprocessing_adapterremoval_adapteroverlap}",
            params.preprocessing_adapterremoval_preserve5p ? "--preserve5p" : "",
            params.preprocessing_trim5p != 0 ? "--trim5p ${params.preprocessing_trim5p}" : "",
            params.preprocessing_trim3p != 0 ? "--trim3p ${params.preprocessing_trim3p}" : "",
            !params.preprocessing_adapterremoval_skipqualitytrimming ? "--trimqualities --minquality ${params.preprocessing_adapterremoval_trimbasequalitymin}" : "",
            !params.preprocessing_adapterremoval_skipntrimming ? "--trimns" : "",
            "--qualitymax ${params.preprocessing_adapterremoval_qualitymax}",
            // FILTERING OPTIONS
            "--minlength ${params.preprocessing_minlength}"
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/preprocessing/adapterremoval/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.preprocessing_savepreprocessedreads
            ],
            [
                // stats
                path: { "${params.outdir}/preprocessing/adapterremoval/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.settings'
            ]
        ]
    }

    withName: CAT_FASTQ_ADAPTERREMOVAL {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
                [
                    // data
                    path: { "${params.outdir}/preprocessing/adapterremoval/data/" },
                    mode: params.publish_dir_mode,
                    pattern: '*.fastq.gz',
                    enabled: params.preprocessing_savepreprocessedreads
                ]
            ]
    }


    //
    // REFERENCE PREPARATION
    //
    withName: GUNZIP_FASTA {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/fasta/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{fasta,fna,fas,fa}',
                enabled: params.save_reference
            ]
        ]
    }

    withName: 'GUNZIP_PMDFASTA|GUNZIP_PMDBED|GUNZIP_SNPBED' {
        publishDir = [
            enabled: false
        ]
    }

    withName: SAMTOOLS_FAIDX {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/fasta/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fai',
                enabled: params.save_reference
            ]
        ]
    }

    withName: PICARD_CREATESEQUENCEDICTIONARY {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/fasta/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.dict',
                enabled: params.save_reference
            ]
        ]
    }

    withName: BOWTIE2_BUILD {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/bowtie2/data/" },
                mode: params.publish_dir_mode,
                pattern: 'bowtie2',
                enabled: params.save_reference
            ]
        ]
    }

    withName: BWA_INDEX {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/bwa/data/" },
                mode: params.publish_dir_mode,
                pattern: 'bwa',
                enabled: params.save_reference
            ]
        ]
    }

    // Reference elongation and indexing for circular mapping
    withName: GUNZIP_ELONGATED_FASTA {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/fasta/data/" },
                mode: params.publish_dir_mode,
                pattern: '*_*[0-9].{fasta,fna,fas,fa}',
                enabled: params.save_reference
            ]
        ]
    }

    withName: CIRCULARMAPPER_CIRCULARGENERATOR {
        tag = { "${meta.id}_${params.fasta_circularmapper_elongationfactor}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/fasta/data/" },
                mode: params.publish_dir_mode,
                pattern: '*_*[0-9].fasta',
                enabled: params.save_reference
            ]
        ]
    }

    withName: BWA_INDEX_CIRCULARISED {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/bwa/data/" },
                mode: params.publish_dir_mode,
                pattern: 'bwa',
                enabled: params.save_reference
            ]
        ]
    }

    //
    // BAM INPUT
    //
    withName: 'SAMTOOLS_FLAGSTATS_BAM_INPUT' {
        // NOTE This step becomes obsolete once a lane-merging step is added for input BAMs.
        // TODO Once a lane-merging step is added for input BAMs, the lane should be dropped from this tag.
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.lane}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/mapping/bam_input/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

    withName: SAMTOOLS_INDEX_BAM_INPUT {
        // NOTE This step becomes obsolete once a lane-merging step is added for input BAMs.
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        publishDir = [
            enabled: false
        ]
    }

    //
    // BAM FILTERING
    //
    withName: FILTER_BAM_FRAGMENT_LENGTH {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = "-l ${params.bamfiltering_minreadlength}"
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/filtered_bams/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.filtered.bam',
                enabled: params.bamfiltering_savefilteredbams
            ]
        ]
    }

    withName: SAMTOOLS_VIEW_BAM_FILTERING {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = [
            "-q ${params.bamfiltering_mappingquality}",
            params.bamfiltering_retainunmappedgenomicbam ? '' : "-F ${params.bamfilter_genomicbamfilterflag}",
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_filtered" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/filtered_bams/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.bam',
                enabled: params.bamfiltering_savefilteredbams
            ]
        ]
    }

    withName: 'SAMTOOLS_LENGTH_FILTER_INDEX|SAMTOOLS_FILTER_INDEX' {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_filtered" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/filtered_bams/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    withName: SAMTOOLS_FLAGSTAT_FILTERED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_filtered" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/read_filtering/filtered_bams/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

    withName: SAMTOOLS_FASTQ_MAPPED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = [
            params.metagenomics_input == 'all' ? '' : '-F 4',
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_mapped" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/fastq/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.bamfiltering_generatemappedfastq
            ]
        ]
    }

    withName: SAMTOOLS_FASTQ_UNMAPPED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = [
            '-f 4',
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_unmapped" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/fastq/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.bamfiltering_generateunmappedfastq
            ]
        ]
    }

    withName: 'CAT_FASTQ_UNMAPPED|CAT_FASTQ_MAPPED' {
        tag = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/fastq/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.preprocessing_savepreprocessedreads
            ]
        ]
    }

    //
    // SHARDING FASTQS
    //
    withName: SEQKIT_SPLIT2 {
        tag = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.prefix = "out"
        ext.args = "-s ${params.fastq_shard_size}"
        publishDir = [
            enabled: false
        ]
    }

    //
    // READ MAPPING
    //
    withName: BWA_ALN {
        tag = { "${meta.id_index}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = { "-n ${params.mapping_bwaaln_n} -k ${params.mapping_bwaaln_k} -l ${params.mapping_bwaaln_l} -o ${params.mapping_bwaaln_o}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.id_index}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: 'BWA_SAMSE|BWA_SAMPE' {
        tag = { "${meta.id_index}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = {
            se_pe_string = meta.single_end ? "SE" : "PE"
            [
                "-r '@RG\\tID:ILLUMINA-${meta.sample_id}_${meta.library_id}\\tSM:${meta.sample_id}\\tLB:${meta.library_id}\\tPL:illumina\\tPU:ILLUMINA-${meta.library_id}-${meta.strandedness}_stranded-${se_pe_string}'"
            ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.id_index}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: ".*MAP:FASTQ_ALIGN_BWAALN:SAMTOOLS_INDEX" {
        tag = { "${meta.id_index}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: BWA_MEM {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = {
            se_pe_string = meta.single_end ? "SE" : "PE"
            [
                "-k ${params.mapping_bwamem_k}",
                "-r ${params.mapping_bwamem_r}",
                "-R '@RG\\tID:ILLUMINA-${meta.sample_id}_${meta.library_id}\\tSM:${meta.sample_id}\\tLB:${meta.library_id}\\tPL:illumina\\tPU:ILLUMINA-${meta.library_id}-${meta.strandedness}_stranded-${se_pe_string}'"
            ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: BOWTIE2_ALIGN {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = {
            se_pe_string = meta.single_end ? "SE" : "PE"
            [
                "-N ${params.mapping_bowtie2_n}",
                "-L ${params.mapping_bowtie2_l}",
                "--trim5 ${params.mapping_bowtie2_trim5}",
                "--trim3 ${params.mapping_bowtie2_trim3}",
                "--maxins ${params.mapping_bowtie2_maxins}",
                "--rg-id ILLUMINA-${meta.library_id}",
                "--rg SM:${meta.sample_id}",
                "--rg LB:${meta.library_id}",
                "--rg PL:illumina",
                "--rg PU:ILLUMINA-${meta.sample_id}_${meta.library_id}-${meta.strandedness}_stranded-${se_pe_string}",
                "${params.mapping_bowtie2_alignmode}" == 'local' ? "--local --${params.mapping_bowtie2_sensitivity}-local" : ( "${params.mapping_bowtie2_alignmode}" == 'end-to-end' ? "--end-to-end --${params.mapping_bowtie2_sensitivity}" : "" )
        ].join(' ').trim()}
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: 'SAMTOOLS_INDEX_MEM|SAMTOOLS_INDEX_BT2' {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: SAMTOOLS_MERGE_LANES {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
        ext.args = { params.run_fastq_sharding ? "-c -p" : "" }
    }

    withName: SAMTOOLS_SORT_MERGED_LANES {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_sorted" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/mapping/${params.mapping_tool}/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bam}'
            ]
        ]
    }

    withName: SAMTOOLS_INDEX_MERGED_LANES {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/mapping/${params.mapping_tool}/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    withName: SAMTOOLS_FLAGSTAT_MERGED_LANES {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_sorted" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/mapping/${params.mapping_tool}/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

    // Circular mapping
    // Configuration for BWA_ALN and BWA_SAMSE/SAMPE is the same as for the non-circular mapping
    withName: CIRCULARMAPPER_REALIGNSAMFILE {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = { params.mapping_circularmapper_circularfilter ? "-f true -x true" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: ".*MAP:CIRCULARMAPPER:FASTQ_ALIGN_BWAALN_ELONGATED:SAMTOOLS_INDEX" {
        tag = { "${meta.id_index}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: ".*MAP:CIRCULARMAPPER:SAMTOOLS_INDEX_REALIGNED" {
        tag = { "${meta.id_index}|${meta.sample_id}_${meta.library_id}_L${meta.lane}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_L${meta.lane}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    //
    // DEDUPLICATION
    //
    withName: PICARD_MARKDUPLICATES {
        tag = { "${meta.reference}:${meta.genomic_region}|${meta.sample_id}_${meta.library_id}" }
        ext.args = [
            "--REMOVE_DUPLICATES",
            "--VALIDATION_STRINGENCY SILENT",
            "--ASSUME_SORTED"
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_${meta.genomic_region}_MarkDuplicates" }
        publishDir = [
            enabled: false
        ]
    }

    withName: DEDUP {
        tag = { "${meta.reference}:${meta.genomic_region}|${meta.sample_id}_${meta.library_id}" }
        ext.args = "--merged"
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_${meta.genomic_region}_DeDup" }
        publishDir = [
            enabled: false
        ]
    }

    withName: SAMTOOLS_MERGE_DEDUPPED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

        withName: ".*DEDUPLICATE:BUILD_INTERVALS" {
        publishDir = [
            enabled: false
        ]
    }

    // Overwrite default SWF tag and prefix
    withName: ".*BAM_SPLIT_BY_REGION:SAMTOOLS_VIEW" {
        tag = { "${meta.reference}:${meta.genomic_region}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_${meta.genomic_region}_dedupped" }
        publishDir = [
            enabled: false
        ]
    }

    withName: ".*BAM_SPLIT_BY_REGION:SAMTOOLS_INDEX" {
        // The BAM_SPLIT_BY_REGION SWF only works with bais, so `params.fasta_largeref` should not be passed to it.
        tag = { "${meta.reference}:${meta.genomic_region}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_${meta.genomic_region}_dedupped" }
        publishDir = [
            enabled: false
        ]
    }

    withName: SAMTOOLS_SORT_DEDUPPED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_dedupped" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/deduplication/${params.deduplication_tool}/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.bam'
            ]
        ]
    }

    withName: SAMTOOLS_INDEX_DEDUPPED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_dedupped" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/deduplication/${params.deduplication_tool}/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    withName: SAMTOOLS_FLAGSTAT_DEDUPPED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_dedupped" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/deduplication/${params.deduplication_tool}/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

    //
    // HOST REMOVAL
    //
    withName: HOST_REMOVAL {
        tag = { "${meta.reference}|${meta_fastqs.sample_id}_${meta_fastqs.library_id}_L${meta_fastqs.lane}" }
        ext.args = {[
            "-m ${params.host_removal_mode}",
            "${meta_fastqs.single_end}" == false && "${params.preprocessing_skippairmerging}" == false ? "-merged" : ""
        ].join(' ').trim()}
        publishDir = [
            [
                // data
                path: { "${params.outdir}/read_filtering/host_removal/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fq.gz'
            ]
        ]
    }

    //
    // AUTHENTICATION
    //
    withName: ENDORSPY {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/endorspy/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.json'
            ]
        ]
    }

    withName: PRESEQ_CCURVE {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = "-B -s ${params.mapstats_preseq_stepsize}"
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/preseq/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: PRESEQ_LCEXTRAP {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args   = [
            "-B",
            "-s ${params.mapstats_preseq_stepsize}",
            "-e ${params.mapstats_preseq_maxextrap}",
            "-x ${params.mapstats_preseq_terms}",
            "-n ${params.mapstats_preseq_bootstrap}",
            "-c ${params.mapstats_preseq_cval}",
            params.mapstats_preseq_defects_mode ? '-D' : ''
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/preseq/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: SAMTOOLS_VIEW_GENOME {
        tag = { "${meta.reference}|${meta.sample_id}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: BEDTOOLS_COVERAGE_DEPTH {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = '-mean -nonamecheck'
        ext.prefix = { "${meta.sample_id}_${meta.reference}_depth" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/bedtools/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: BEDTOOLS_COVERAGE_BREADTH {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = '-nonamecheck'
        ext.prefix = { "${meta.sample_id}_${meta.reference}_breadth" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/bedtools/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: ANGSD_DOCOUNTS {
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = [
            "-iCounts 1",
            "-r ${params.contamination_estimation_angsd_chrom_name}:${params.contamination_estimation_angsd_range_from}-${params.contamination_estimation_angsd_range_to}",
            "-minMapQ ${params.contamination_estimation_angsd_mapq}",
            "-minQ ${params.contamination_estimation_angsd_minq}"
        ].join(' ').trim()
        publishDir = [
            enabled: false
        ]
    }

    withName: ANGSD_CONTAMINATION {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/angsd_nuclear_contamination/stats/"},
                mode: params.publish_dir_mode,
                pattern: '*.txt'
            ]
        ]
    }

    withName: PRINT_CONTAMINATION_ANGSD {
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/angsd_nuclear_contamination/stats/"},
                mode: params.publish_dir_mode,
                pattern: 'nuclear_contamination*'
            ]
        ]
    }

    withName: MTNUCRATIO {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: 'QUALIMAP_BAMQC_WITHBED|QUALIMAP_BAMQC_NOBED' {
        tag = { "${meta.reference}|${meta.sample_id}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/qualimap/stats/${meta.reference}/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: DAMAGEPROFILER {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args   = [
            "-l ${params.damagecalculation_damageprofiler_length}",
            "-t ${params.damagecalculation_xaxis}",
            "-yaxis_dp_max ${params.damagecalculation_yaxis}"
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/damageprofiler/stats/" },
                mode: params.publish_dir_mode,
                saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
            ]
        ]
    }

    withName: CALCULATE_MAPDAMAGE2 {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { [
            "--no-stats",
            "-y ${params.damagecalculation_yaxis}",
            params.damagecalculation_mapdamage_downsample != 0 ? "-n ${params.damagecalculation_mapdamage_downsample} --downsample-seed=1" : "",
            { meta.strandedness } == "single" ? '--single-stranded' : '',
            "-m ${params.damagecalculation_xaxis}"
        ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/mapdamage2/stats/" },
                mode: params.publish_dir_mode,
                pattern: 'results_*/*'
            ]
        ]
    }

    withName: SAMTOOLS_DEPTH_SEXDETERRMINE {
        tag = { "${meta1.reference}|${meta1.sample_id}" }
        ext.prefix = { "${meta2.id}_samtoolsdepth" }
        ext.args = '-aa -q30 -Q30 -H'
        publishDir = [
            enabled: false
        ]
    }

    withName: SEXDETERRMINE {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.reference}_sexdeterrmine" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/authentication/sexdeterrmine/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*{_sexdeterrmine}*',
            ]
        ]
    }

    //
    // DAMAGE MANIPULATION
    //
    withName: BEDTOOLS_MASKFASTA {
        ext.prefix = { "${meta.id}.masked" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/references/masked_reference/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.masked.fa'
            ]
        ]
    }

    withName: MAPDAMAGE2 {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { [
            "--rescale",
            // "--rescale-out=${task.ext.prefix}.bam", // This doesn't work, because the output is expected to be in the mapdamage results dir.
            "--seq-length=${params.damage_manipulation_rescale_seqlength}",
            params.damage_manipulation_rescale_length_3p != 0 ? "--rescale-length-3p=${params.damage_manipulation_rescale_length_3p}" : "",
            params.damage_manipulation_rescale_length_5p != 0 ? "--rescale-length-5p=${params.damage_manipulation_rescale_length_5p}" : "",
            { meta.strandedness } == "single" ? '--single-stranded' : ''
        ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_rescaled" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/damage_manipulation/mapdamage2/data/" },
                mode: params.publish_dir_mode,
                pattern: 'result*/*.rescaled.bam',
                saveAs: { "${meta.sample_id}_${meta.library_id}_${meta.reference}_rescaled.bam" } // Needed to save the bam directly in the output dir instead of within the mapdamage results dir.
            ],
            [
                // stats
                path: { "${params.outdir}/damage_manipulation/mapdamage2/stats/" },
                mode: params.publish_dir_mode,
                pattern: 'results_*/Stats_out_MCMC_*'
            ]
        ]
    }

    withName: SAMTOOLS_INDEX_DAMAGE_RESCALED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/damage_manipulation/mapdamage2/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}',
                saveAs: { "${meta.sample_id}_${meta.library_id}_${meta.reference}_rescaled.bam.bai" } // ext.prefix cannot be used here, so rename for publishing instead.
            ]
        ]
    }

    withName: PMDTOOLS_FILTER {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args2 = {[
            "${meta.damage_treatment}" == 'none' ? '--UDGminus' : "${meta.damage_treatment}" == 'half' ? '--UDGhalf' : '--UDGplus' ,
            "${meta.strandedness}" == 'single' ? '--ss' : ''
        ].join(' ').trim()}
        ext.args3 = { "-h" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_pmdfiltered" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/damage_manipulation/pmdtools/data/" },
                mode: params.publish_dir_mode,
                pattern: '*_pmdfiltered.bam'
            ]
        ]
    }

    withName: SAMTOOLS_INDEX_DAMAGE_FILTERED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/damage_manipulation/pmdtools/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    withName: SAMTOOLS_FLAGSTAT_DAMAGE_FILTERED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_pmdfiltered" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/damage_manipulation/pmdtools/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

    withName: BAMUTIL_TRIMBAM {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = [
            params.damage_manipulation_bamutils_softclip ? '-c' : ''
        ].join(' ').trim()
        ext.prefix = { params.run_pmd_filtering ? "${meta.sample_id}_${meta.library_id}_${meta.reference}_pmdfiltered_trimmed" : "${meta.sample_id}_${meta.library_id}_${meta.reference}_trimmed" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/damage_manipulation/bamutils_trimbam/data/" },
                mode: params.publish_dir_mode,
                pattern: '*_trimmed.bam'
            ]
        ]
    }

    withName: SAMTOOLS_INDEX_DAMAGE_TRIMMED {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/damage_manipulation/bamutils_trimbam/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    //
    // METAGENOMIC SCREENING
    //
    withName: PRINSEQPLUSPLUS {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args =  [
                params.metagenomics_prinseq_mode == 'dust' ? "-lc_dust=${params.metagenomics_prinseq_dustscore}" : "-lc_entropy=${params.metagenomics_complexity_entropy}",
                "-trim_qual_left=0 -trim_qual_left=0 -trim_qual_window=0 -trim_qual_step=0",
            ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_complexity" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/metagenomics/prinseq/data/" },
                mode: params.publish_dir_mode,
                pattern: '*{_good_out.fastq.gz,_good_out_R1.fastq.gz,_good_out_R2.fastq.gz}',
                enabled: params.metagenomics_complexity_savefastq
            ],
            [
                // stats
                path: { "${params.outdir}/metagenomics/prinseq/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*log',
                enabled: params.metagenomics_complexity_savefastq
            ]
        ]
    }

    withName: BBMAP_BBDUK {
        tag = { "${meta.reference}|${meta.sample_id}_${meta.library_id}" }
        ext.args = { "entropymask=f entropy=${params.metagenomics_complexity_entropy}" }
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}_complexity" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/metagenomics/bbduk/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz',
                enabled: params.metagenomics_complexity_savefastq
            ],
            [
                // stats
                path: { "${params.outdir}/metagenomics/bbduk/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.log',
                enabled: params.metagenomics_complexity_savefastq
            ]
        ]
    }

    withName: MALT_RUN {
        // TODO needs a tag
        ext.args = [
            "-m ${params.metagenomics_malt_mode}",
            "-at ${params.metagenomics_malt_alignmentmode}",
            "-top ${params.metagenomics_malt_toppercent}",
            "-id ${params.metagenomics_malt_minpercentidentity}",
            "-mq ${params.metagenomics_malt_maxqueries}",
            "--memoryMode ${params.metagenomics_malt_memorymode}",
            params.metagenomics_malt_minsupportmode == "percent" ? "-supp ${params.metagenomics_malt_minsupportpercent}" : "-sup ${params.metagenomics_malt_minsupportreads}",
            params.metagenomics_malt_savereads ? "--alignments ./" : ""
        ].join(' ').trim()
        ext.prefix = { "${meta.label}_${meta.id}-run" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/metagenomics/malt/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{rma6,sam.gz}'
            ],
            [
                // stats
                path: { "${params.outdir}/metagenomics/malt/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.log'
            ]
        ]
    }

    withName: CAT_CAT_MALT {
        // TODO needs a tag
        ext.prefix = { "${meta.id}_runtime_log_concatenated.log" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/metagenomics/malt/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.log'
            ]
        ]
    }

    withName: UNTAR_METAGENOMICS {
        publishDir = [
            enabled: false
        ]
    }

    withName: KRAKEN2_KRAKEN2 {
        ext.args = [
            params.metagenomics_kraken2_saveminimizers ? "--report-minimizer-data" : ""
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/metagenomics/kraken2/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz'
            ],
            [
                // stats
                path: { "${params.outdir}/metagenomics/kraken2/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.txt'
            ]
        ]
    }

    withName: KRAKENUNIQ_PRELOADEDKRAKENUNIQ {
        publishDir = [
            [
                // data
                path: { "${params.outdir}/metagenomics/krakenuniq/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.fastq.gz'
            ],
            [
                // stats
                path: { "${params.outdir}/metagenomics/krakenuniq/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.txt'
            ]
        ]
    }

    withName: METAPHLAN_METAPHLAN {
        ext.prefix = { "${meta.sample_id}_${meta.library_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/metagenomics/metaphlan/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.{biom,txt}'
            ]
        ]
    }

    withName: MALTEXTRACT {
        // TODO needs a tag
        ext.args = [
            "-f ${params.metagenomics_maltextract_filter}",
            "-a ${params.metagenomics_maltextract_toppercent}",
            "--minPI ${params.metagenomics_maltextract_minpercentidentity}",
            params.metagenomics_maltextract_destackingoff ? "--destackingOff" : "",
            params.metagenomics_maltextract_downsamplingoff ? "--downSampOff" : "",
            params.metagenomics_maltextract_duplicateremovaloff ? "--dupRemOff" : "",
            params.metagenomics_maltextract_matches ? "--matches" : "",
            params.metagenomics_maltextract_megansummary ? "--meganSummary" : "",
            params.metagenomics_maltextract_usetopalignment ?  "--useTopAlignment" : "",
            { meta.strandedness } == "single" ? '--singleStranded' : '',
        ].join(' ').trim()
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/metagenomics/maltextract/stats/" },
                mode: params.publish_dir_mode,
                pattern: 'results',
                saveAs: { "${meta.id}" }
            ]
        ]
    }

    withName: MEGAN_RMA2INFO {
        tag = {"${meta.id}"}
        ext.args = "-c2c Taxonomy"
        ext.prefix = { "${meta.id}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/metagenomics/megan_summaries/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.{txt.gz,megan}'
            ]
        ]
    }

    withName: AMPS {
        publishDir = [
            path: { "${params.outdir}/metagenomics/amps/stats/" },
            mode: params.publish_dir_mode,
            pattern: 'results'
        ]
        errorStrategy = 'ignore' // required as it fails the run for low reads: https://github.com/rhuebler/HOPS/issues/9
    }

    withName: TAXPASTA_MERGE {
        ext.args = { "--profiler ${meta.profiler} --output ${meta.profiler}_taxpasta_table.tsv" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/metagenomics/taxpasta/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.{csv,tsv,ods,xlsx,arrow,parquet,biom}'
            ]
        ]
    }

    withName: TAXPASTA_STANDARDISE {
        ext.args = { "--profiler ${meta.profiler} --output ${meta.profiler}taxpasta_table.tsv" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/metagenomics/taxpasta/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.{csv,tsv,ods,xlsx,arrow,parquet,biom}'
            ]
        ]
    }

    //
    // LIBRARY MERGE
    //
    withName: ".*MERGE_LIBRARIES:SAMTOOLS_MERGE_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}_unsorted" }
        publishDir = [
            enabled: false
        ]
    }

    withName: ".*MERGE_LIBRARIES:SAMTOOLS_SORT_MERGED_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/final_bams/raw/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.bam'
            ]
        ]
    }

    withName: ".*MERGE_LIBRARIES:SAMTOOLS_INDEX_MERGED_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/final_bams/raw/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    withName: ".*MERGE_LIBRARIES:SAMTOOLS_FLAGSTAT_MERGED_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/final_bams/raw/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

        withName: ".*MERGE_LIBRARIES:SAMTOOLS_MERGE_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}_unsorted" }
        publishDir = [
            enabled: false
        ]
    }

    withName: ".*MERGE_LIBRARIES_GENOTYPING:SAMTOOLS_SORT_MERGED_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/final_bams/${params.genotyping_source}/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.bam'
            ]
        ]
    }

    withName: ".*MERGE_LIBRARIES_GENOTYPING:SAMTOOLS_INDEX_MERGED_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = { params.fasta_largeref ? "-c" : "" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/final_bams/${params.genotyping_source}/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bai,csi}'
            ]
        ]
    }

    withName: ".*MERGE_LIBRARIES_GENOTYPING:SAMTOOLS_FLAGSTAT_MERGED_LIBRARIES" {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/final_bams/${params.genotyping_source}/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.flagstat'
            ]
        ]
    }

    //
    // GENOTYPING
    //
    withName: SAMTOOLS_MPILEUP_PILEUPCALLER {
        tag = { "${meta.reference}|${meta.strandedness}" }
        ext.args = [
            "-B",
            "-q ${params.genotyping_pileupcaller_min_base_quality}",
            "-Q ${params.genotyping_pileupcaller_min_map_quality}",
        ].join(' ').trim()
        ext.prefix = { "${meta.strandedness}_${meta.reference}" }
        publishDir = [
            enabled: false
        ]
    }

    withName: SEQUENCETOOLS_PILEUPCALLER {
        tag = { "${meta.reference}|${meta.strandedness}" }
        ext.args = {[
            "--${params.genotyping_pileupcaller_method}",
            params.genotyping_pileupcaller_transitions_mode == "SkipTransitions" ? "--skipTransitions" : params.genotyping_pileupcaller_transitions_mode == "TransitionsMissing" ? "--transitionsMissing" : "",
            "${meta.strandedness}" == 'single' ? "--singleStrandMode" : "" ,
            "--sampleNames", meta.sample_id.join(","),
            "-e pileupcaller.${meta.strandedness}.${meta.reference}"
        ].join(' ').trim() }
        ext.prefix = { "${meta.strandedness}_${meta.reference}" }
        publishDir = [
            enabled: false // Not published because the output goes through COLLECT_GENOTYPES
        ]
    }

    withName: COLLECT_GENOTYPES {
        tag = { "${meta.reference}" }
        ext.prefix = { "pileupcaller_genotypes_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/pileupcaller/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{geno,snp,ind}',
            ]
        ]
    }

    withName: EIGENSTRATDATABASETOOLS_EIGENSTRATSNPCOVERAGE {
        tag = { "${meta.reference}" }
        ext.args = { "-j ${prefix}.json" }
        ext.prefix = { "pileupcaller_genotypes_${meta.reference}_coverage" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/genotyping/pileupcaller/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.tsv',
            ]
        ]
    }

    withName: GATK_REALIGNERTARGETCREATOR {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = [
            params.genotyping_gatk_ug_defaultbasequalities > 0 ? "--defaultBaseQualities ${params.genotyping_gatk_ug_defaultbasequalities}" : "", // Empty string since GATK complains if its default of -1 is provided.
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.reference}_realigntarget" }
        publishDir = [
            enabled: false
        ]
    }

    withName: GATK_INDELREALIGNER {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = [
            params.genotyping_gatk_ug_defaultbasequalities > 0 ? "--defaultBaseQualities ${params.genotyping_gatk_ug_defaultbasequalities}" : "", // Empty string since GATK complains if its default of -1 is provided.
        ].join(' ').trim()
        ext.prefix = { "${meta.sample_id}_${meta.reference}_realigned" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/gatk_ug/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{bam,bai}',
                enabled: params.genotyping_gatk_ug_keeprealignbam,
            ]
        ]
    }

    withName: GATK_UNIFIEDGENOTYPER {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = {[
            "--sample_ploidy ${meta2.ploidy}",
            "-stand_call_conf ${params.genotyping_gatk_call_conf}",
            "-dcov ${params.genotyping_gatk_ug_downsample}",
            "--output_mode ${params.genotyping_gatk_ug_out_mode}",
            "--genotype_likelihoods_model ${params.genotyping_gatk_ug_genotype_mode}",
            params.genotyping_gatk_ug_defaultbasequalities > 0 ? "--defaultBaseQualities ${params.genotyping_gatk_ug_defaultbasequalities}" : "", // Empty string since GATK complains if its default of -1 is provided.
        ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/gatk_ug/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.vcf.gz',
            ]
        ]
    }

    withName: BCFTOOLS_INDEX_UG {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = "--tbi" //tbi indices for consistency with GATK HC
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/gatk_ug/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.vcf.gz.tbi',
            ]
        ]
    }

    withName: GATK4_HAPLOTYPECALLER {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = {[
            // Option names have changed from underscore_separated to hyphen-separated in GATK4
            "--sample-ploidy ${meta2.ploidy}",
            "-stand-call-conf ${params.genotyping_gatk_call_conf}",
            "--output-mode ${params.genotyping_gatk_hc_out_mode}",
            "--emit-ref-confidence ${params.genotyping_gatk_hc_emitrefconf}",
        ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/gatk_hc/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{vcf.gz,vcf.gz.tbi}',
            ]
        ]
    }

    withName: FREEBAYES {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = {[
            "-p ${ref_meta.ploidy}",
            "-C ${params.genotyping_freebayes_min_alternate_count}",
            params.genotyping_freebayes_skip_coverage != 0 ? "-g ${params.genotyping_freebayes_skip_coverage}" : "",
        ].join(' ').trim() }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/freebayes/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.vcf.gz',
            ]
        ]
    }

    withName: BCFTOOLS_INDEX_FREEBAYES {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.args = "--tbi" //tbi indices for consistency with GATK HC
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/freebayes/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.vcf.gz.tbi',
            ]
        ]
    }

    withName: BCFTOOLS_STATS_GENOTYPING {
        tag = { "${meta.reference}|${meta.sample_id}" }
        ext.prefix = { "${meta.sample_id}_${meta.reference}" }
        publishDir = [
            [
                // stats
                path: { "${params.outdir}/genotyping/bcftools/stats/" },
                mode: params.publish_dir_mode,
                pattern: '*.txt',
            ]
        ]
    }

    withName: ANGSD_GL {
        tag = { "${meta.reference}" }
        ext.args = {
            gl_model = params.genotyping_angsd_glmodel == 'samtools' ? 1 : params.genotyping_angsd_glmodel == 'gatk' ? 2 : params.genotyping_angsd_glmodel == 'soapsnp' ? 3 : 4
            gl_format = params.genotyping_angsd_glformat == 'binary' ? 1 : params.genotyping_angsd_glformat == 'beagle_binary' ? 2 : params.genotyping_angsd_glformat == 'binary_three' ? 3 : 4
            [
                ( gl_format == 2 || gl_format == 3 ) ? '-doMajorMinor 1': '',
                "-GL ${gl_model}",
                "-doGlf ${gl_format}",
            ].join(' ').trim()
        }
        ext.prefix = { "angsd_${meta.reference}" }
        publishDir = [
            [
                // data
                path: { "${params.outdir}/genotyping/angsd/data/" },
                mode: params.publish_dir_mode,
                pattern: '*.{glf,beagle}.gz',
            ]
        ]
    }

    withName: MULTIVCFANALYZER {
        tag = { "${meta.reference}" }
        ext.prefix = { "multivcfanalyzer_${meta.reference}" }
        publishDir = [
            path: { "${params.outdir}/consensus_sequence/" },
            mode: params.publish_dir_mode,
            enabled: true,
            pattern: '*.{fasta.gz,tsv,txt}'
        ]
    }
}