kkonganti@17: process SALMON_QUANT {
kkonganti@17:     tag "$meta.id"
kkonganti@17:     label "process_micro"
kkonganti@17: 
kkonganti@17:     module (params.enable_module ? "${params.swmodulepath}${params.fs}salmon${params.fs}1.10.0" : null)
kkonganti@17:     conda (params.enable_conda ? 'conda-forge::libgcc-ng bioconda::salmon=1.10.1' : null)
kkonganti@17:     container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
kkonganti@17:         'https://depot.galaxyproject.org/singularity/salmon:1.10.1--h7e5ed60_1' :
kkonganti@17:         'quay.io/biocontainers/salmon:1.10.1--h7e5ed60_1' }"
kkonganti@17:     input:
kkonganti@17:         tuple val(meta), path(reads_or_bam), path(index_or_tr_fasta)
kkonganti@17: 
kkonganti@17:     output:
kkonganti@17:         tuple val(meta), path("${meta.id}_salmon_res"), emit: results
kkonganti@17:         path  "versions.yml"                          , emit: versions
kkonganti@17: 
kkonganti@17:     when:
kkonganti@17:         task.ext.when == null || task.ext.when
kkonganti@17: 
kkonganti@17:     script:
kkonganti@17:         def args = task.ext.args   ?: ''
kkonganti@17:         def prefix   = task.ext.prefix ?: "${meta.id}_salmon_res"
kkonganti@17:         def reference   = "--index $index_or_tr_fasta"
kkonganti@17:         def lib_type = (meta.salmon_lib_type ?: '')
kkonganti@17:         def alignment_mode = (meta.salmon_alignment_mode ?: '')
kkonganti@17:         def gtf = (meta.salmon_gtf ? "--geneMap ${meta.salmon_gtf}" : '')
kkonganti@17:         def input_reads =(meta.single_end ? "-r $reads_or_bam" : "-1 ${reads_or_bam[0]} -2 ${reads_or_bam[1]}")
kkonganti@17: 
kkonganti@17:         // Use path(reads_or_bam) to point to BAM and path(index_or_tr_fasta) to point to transcript fasta
kkonganti@17:         // if using salmon DSL2 module in alignment-based mode.
kkonganti@17:         // By default, this module will be run in selective-alignment-based mode of salmon.
kkonganti@17:         if (alignment_mode) {
kkonganti@17:             reference   = "-t $index_or_tr_fasta"
kkonganti@17:             input_reads = "-a $reads_or_bam"
kkonganti@17:         }
kkonganti@17: 
kkonganti@17:         def strandedness_opts = [
kkonganti@17:             'A', 'U', 'SF', 'SR',
kkonganti@17:             'IS', 'IU' , 'ISF', 'ISR',
kkonganti@17:             'OS', 'OU' , 'OSF', 'OSR',
kkonganti@17:             'MS', 'MU' , 'MSF', 'MSR'
kkonganti@17:         ]
kkonganti@17: 
kkonganti@17:         def strandedness =  'A'
kkonganti@17: 
kkonganti@17:         if (lib_type) {
kkonganti@17:             if (strandedness_opts.contains(lib_type)) {
kkonganti@17:                 strandedness = lib_type
kkonganti@17:             } else {
kkonganti@17:                 log.info "[Salmon Quant] Invalid library type specified '--libType=${lib_type}', defaulting to auto-detection with '--libType=A'."
kkonganti@17:             }
kkonganti@17:         } else {
kkonganti@17:             strandedness = meta.single_end ? 'U' : 'IU'
kkonganti@17:             if (meta.strandedness == 'forward') {
kkonganti@17:                 strandedness = meta.single_end ? 'SF' : 'ISF'
kkonganti@17:             } else if (meta.strandedness == 'reverse') {
kkonganti@17:                 strandedness = meta.single_end ? 'SR' : 'ISR'
kkonganti@17:             }
kkonganti@17:         }
kkonganti@17:         """
kkonganti@17:         salmon quant \\
kkonganti@17:             --threads $task.cpus \\
kkonganti@17:             --libType=$strandedness \\
kkonganti@17:             $gtf \\
kkonganti@17:             $args \\
kkonganti@17:             -o $prefix \\
kkonganti@17:             $reference \\
kkonganti@17:             $input_reads
kkonganti@17: 
kkonganti@17:         cat <<-END_VERSIONS > versions.yml
kkonganti@17:         "${task.process}":
kkonganti@17:             salmon: \$(echo \$(salmon --version) | sed -e "s/salmon //g")
kkonganti@17:         END_VERSIONS
kkonganti@17:         """
kkonganti@17: }