2020-03-06 13:21:17 +00:00
|
|
|
nextflow.preview.dsl=2
|
|
|
|
|
|
|
|
process FASTQ_SCREEN {
|
2020-07-11 11:42:13 +00:00
|
|
|
|
2020-07-27 10:11:42 +00:00
|
|
|
publishDir "$outputdir",
|
|
|
|
mode: "link", overwrite: true
|
|
|
|
|
2020-07-11 11:42:13 +00:00
|
|
|
// depending on the number of genomes and the type of genome (e.g. plants!), memory needs to be ample!
|
|
|
|
// label 'bigMem'
|
|
|
|
// label 'multiCore'
|
2020-03-06 13:21:17 +00:00
|
|
|
|
|
|
|
input:
|
2020-07-11 11:42:13 +00:00
|
|
|
tuple val(name), path(reads)
|
2020-07-27 10:11:42 +00:00
|
|
|
val outputdir
|
2020-07-11 11:42:13 +00:00
|
|
|
// fastq_screen_args are best passed in to the workflow in the following manner:
|
|
|
|
// --fastq_screen_args="--subset 200000 --force"
|
2020-07-27 10:11:42 +00:00
|
|
|
val fastq_screen_args
|
|
|
|
val verbose
|
2020-03-06 13:21:17 +00:00
|
|
|
|
2020-07-11 11:42:13 +00:00
|
|
|
output:
|
|
|
|
path "*png", emit: png
|
|
|
|
path "*html", emit: html
|
|
|
|
path "*txt", emit: report
|
2020-03-06 13:21:17 +00:00
|
|
|
|
|
|
|
script:
|
2020-07-11 11:42:13 +00:00
|
|
|
println(name)
|
|
|
|
println(reads)
|
|
|
|
println(outputdir)
|
|
|
|
if (verbose){
|
|
|
|
println ("[MODULE] FASTQ SCREEN ARGS: "+ fastq_screen_args)
|
|
|
|
}
|
|
|
|
|
|
|
|
"""
|
|
|
|
module load fastq_screen
|
|
|
|
fastq_screen $fastq_screen_args $reads
|
|
|
|
"""
|
|
|
|
|
|
|
|
}
|