2021-05-07 15:37:35 +00:00
|
|
|
process ARTIC_MINION {
|
|
|
|
tag "$meta.id"
|
|
|
|
label 'process_high'
|
|
|
|
|
|
|
|
conda (params.enable_conda ? "bioconda::artic=1.2.1" : null)
|
2021-11-26 07:58:40 +00:00
|
|
|
container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
|
|
|
|
'https://depot.galaxyproject.org/singularity/artic:1.2.1--py_0' :
|
|
|
|
'quay.io/biocontainers/artic:1.2.1--py_0' }"
|
2021-05-07 15:37:35 +00:00
|
|
|
|
|
|
|
input:
|
|
|
|
tuple val(meta), path(fastq)
|
|
|
|
path fast5_dir
|
|
|
|
path sequencing_summary
|
|
|
|
path ("primer-schemes/${scheme}/V${scheme_version}/${scheme}.reference.fasta")
|
|
|
|
path ("primer-schemes/${scheme}/V${scheme_version}/${scheme}.scheme.bed")
|
2022-02-18 17:36:09 +00:00
|
|
|
path medaka_model_file
|
|
|
|
val medaka_model_string
|
2021-05-07 15:37:35 +00:00
|
|
|
val scheme
|
|
|
|
val scheme_version
|
|
|
|
|
|
|
|
output:
|
|
|
|
tuple val(meta), path("${prefix}.*") , emit: results
|
|
|
|
tuple val(meta), path("${prefix}.sorted.bam") , emit: bam
|
|
|
|
tuple val(meta), path("${prefix}.sorted.bam.bai") , emit: bai
|
|
|
|
tuple val(meta), path("${prefix}.trimmed.rg.sorted.bam") , emit: bam_trimmed
|
|
|
|
tuple val(meta), path("${prefix}.trimmed.rg.sorted.bam.bai") , emit: bai_trimmed
|
|
|
|
tuple val(meta), path("${prefix}.primertrimmed.rg.sorted.bam") , emit: bam_primertrimmed
|
|
|
|
tuple val(meta), path("${prefix}.primertrimmed.rg.sorted.bam.bai"), emit: bai_primertrimmed
|
|
|
|
tuple val(meta), path("${prefix}.consensus.fasta") , emit: fasta
|
|
|
|
tuple val(meta), path("${prefix}.pass.vcf.gz") , emit: vcf
|
|
|
|
tuple val(meta), path("${prefix}.pass.vcf.gz.tbi") , emit: tbi
|
|
|
|
tuple val(meta), path("*.json"), optional:true , emit: json
|
2021-10-01 13:04:56 +00:00
|
|
|
path "versions.yml" , emit: versions
|
2021-05-07 15:37:35 +00:00
|
|
|
|
2022-02-04 08:53:32 +00:00
|
|
|
when:
|
|
|
|
task.ext.when == null || task.ext.when
|
|
|
|
|
2021-05-07 15:37:35 +00:00
|
|
|
script:
|
2021-12-07 14:00:43 +00:00
|
|
|
def args = task.ext.args ?: ''
|
|
|
|
prefix = task.ext.prefix ?: "${meta.id}"
|
2021-05-07 15:37:35 +00:00
|
|
|
def version = scheme_version.toString().toLowerCase().replaceAll('v','')
|
2021-11-26 07:58:40 +00:00
|
|
|
def fast5 = fast5_dir ? "--fast5-directory $fast5_dir" : ""
|
|
|
|
def summary = sequencing_summary ? "--sequencing-summary $sequencing_summary" : ""
|
2021-05-07 15:37:35 +00:00
|
|
|
def model = ""
|
2021-11-26 07:58:40 +00:00
|
|
|
if (args.tokenize().contains('--medaka')) {
|
2021-05-07 15:37:35 +00:00
|
|
|
fast5 = ""
|
|
|
|
summary = ""
|
2022-02-18 17:36:09 +00:00
|
|
|
model = medaka_model_file ? "--medaka-model ./$medaka_model_file" : "--medaka-model $medaka_model_string"
|
2021-05-07 15:37:35 +00:00
|
|
|
}
|
2022-01-13 01:05:01 +00:00
|
|
|
def hd5_plugin_path = task.ext.hd5_plugin_path ? "export HDF5_PLUGIN_PATH=" + task.ext.hd5_plugin_path : "export HDF5_PLUGIN_PATH=/usr/local/lib/python3.6/site-packages/ont_fast5_api/vbz_plugin"
|
2021-05-07 15:37:35 +00:00
|
|
|
"""
|
2022-01-11 21:10:26 +00:00
|
|
|
$hd5_plugin_path
|
|
|
|
|
2021-05-07 15:37:35 +00:00
|
|
|
artic \\
|
|
|
|
minion \\
|
2021-11-26 07:58:40 +00:00
|
|
|
$args \\
|
2021-05-07 15:37:35 +00:00
|
|
|
--threads $task.cpus \\
|
|
|
|
--read-file $fastq \\
|
|
|
|
--scheme-directory ./primer-schemes \\
|
|
|
|
--scheme-version $version \\
|
|
|
|
$model \\
|
|
|
|
$fast5 \\
|
|
|
|
$summary \\
|
|
|
|
$scheme \\
|
|
|
|
$prefix
|
|
|
|
|
2021-09-27 08:41:24 +00:00
|
|
|
cat <<-END_VERSIONS > versions.yml
|
2021-11-26 07:58:40 +00:00
|
|
|
"${task.process}":
|
|
|
|
artic: \$(artic --version 2>&1 | sed 's/^.*artic //; s/ .*\$//')
|
2021-09-27 08:41:24 +00:00
|
|
|
END_VERSIONS
|
2021-05-07 15:37:35 +00:00
|
|
|
"""
|
|
|
|
}
|