nf-core_modules/software/bowtie/align/main.nf

43 lines
1.3 KiB
Text
Raw Normal View History

2020-12-08 11:57:00 +00:00
// Import generic module functions
include { initOptions; saveFiles; getSoftwareName } from './functions'
params.options = [:]
def options = initOptions(params.options)
process BOWTIE_ALIGN {
tag "$meta.id"
label 'process_high'
publishDir "${params.outdir}",
mode: params.publish_dir_mode,
saveAs: { filename -> saveFiles(filename:filename, options:params.options, publish_dir:getSoftwareName(task.process), publish_id:meta.id) }
conda (params.enable_conda ? "bioconda::bowtie=1.3.0" : null)
container "quay.io/biocontainers/bowtie:1.3.0--py38hed8969a_1"
input:
tuple val(meta), path(reads)
path index
output:
tuple val(meta), path("*.sam"), emit: sam
tuple val(meta), path("*.out"), emit: log
2020-12-09 14:58:56 +00:00
path "bowtie.version.txt", emit: version
2020-12-08 11:57:00 +00:00
script:
def software = getSoftwareName(task.process)
def prefix = options.suffix ? "${meta.id}${options.suffix}" : "${meta.id}"
def unaligned = params.save_unaligned ? "--un ${prefix}.unmapped" : ''
2020-12-08 11:57:00 +00:00
"""
INDEX=`find -L ./ -name "*.1.ebwt" | sed 's/.1.ebwt//'`
bowtie \\
--threads $task.cpus \\
$options.args \\
$INDEX \\
-q ${reads} \\
$unaligned \\
> ${prefix}.sam 2> ${prefix}.out
2020-12-08 11:57:00 +00:00
bowtie --version | head -n 1 | cut -d" " -f3 > ${software}.version.txt
"""
}