mirror of
https://github.com/MillironX/nf-core_modules.git
synced 2024-11-14 05:43:08 +00:00
16d20a7cc4
* New last/lastdb module to index sequences before alignment. The `lastdb` command creates a sequence index for the LAST aligner (https://gitlab.com/mcfrith/last). Input can be in FASTA or FASTQ format, and compression is handled automagically. DNA or protein sequences can be indexed. The sequence index is a collection of files sharing the same basename. This module sets the basename to the sample identifier (`$meta.id`) and creates the index in a directory always called `lastdb`. The module's output channel then conveys a copy of the metadata and the path to the `lastdb` directory. Other modules will follow (see Issue #464). The LAST aligner can align proteins to proteins, DNA to DNA and can translate DNA align to proteins. * Remove trailing whitespace. * Apply suggestions from code review Co-authored-by: Harshil Patel <drpatelh@users.noreply.github.com> * Update as suggested in PR. * Attempt to pass linting. Co-authored-by: Harshil Patel <drpatelh@users.noreply.github.com> Co-authored-by: Harshil Patel <drpatelh@users.noreply.github.com>
41 lines
1.3 KiB
Text
41 lines
1.3 KiB
Text
// Import generic module functions
|
|
include { initOptions; saveFiles; getSoftwareName } from './functions'
|
|
|
|
params.options = [:]
|
|
options = initOptions(params.options)
|
|
|
|
process LAST_LASTDB {
|
|
tag "$meta.id"
|
|
label 'process_medium'
|
|
publishDir "${params.outdir}",
|
|
mode: params.publish_dir_mode,
|
|
saveAs: { filename -> saveFiles(filename:filename, options:params.options, publish_dir:getSoftwareName(task.process), meta:meta, publish_by_meta:['id']) }
|
|
|
|
conda (params.enable_conda ? "bioconda::last=1219" : null)
|
|
if (workflow.containerEngine == 'singularity' && !params.singularity_pull_docker_container) {
|
|
container "https://depot.galaxyproject.org/singularity/last:1219--h2e03b76_0"
|
|
} else {
|
|
container "quay.io/biocontainers/last:1219--h2e03b76_0"
|
|
}
|
|
|
|
input:
|
|
tuple val(meta), path(fastx)
|
|
|
|
output:
|
|
tuple val(meta), path("lastdb"), emit: index
|
|
path "*.version.txt" , emit: version
|
|
|
|
script:
|
|
def software = getSoftwareName(task.process)
|
|
def prefix = options.suffix ? "${meta.id}${options.suffix}" : "${meta.id}"
|
|
"""
|
|
mkdir lastdb
|
|
lastdb \\
|
|
$options.args \\
|
|
-P $task.cpus \\
|
|
lastdb/${prefix} \\
|
|
$fastx
|
|
|
|
echo \$(lastdb --version 2>&1) | sed 's/lastdb //' > ${software}.version.txt
|
|
"""
|
|
}
|