-
Laurent Modolo authoredLaurent Modolo authored
multiqc_paired.nf 800 B
params.fastq = "$baseDir/data/fastq/*_{1,2}.fastq"
log.info "fastq files : ${params.fastq}"
Channel
.fromFilePairs( params.fastq )
.ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" }
.set { fastq_files }
process fastqc_fastq {
tag "$pair_id"
publishDir "results/fastq/fastqc/", mode: 'copy'
input:
set pair_id, file(reads) from fastq_files
output:
file "*.{zip,html}" into fastqc_report
script:
"""
fastqc --quiet --threads ${task.cpus} --format fastq --outdir ./ \
${reads[0]} ${reads[1]}
"""
}
process multiqc {
tag "$report[0].baseName"
publishDir "results/fastq/multiqc/", mode: 'copy'
cpus = 1
input:
file report from fastqc_report.collect()
output:
file "*multiqc_*" into multiqc_report
script:
"""
multiqc -f .
"""
}