2022-02-19 11:36:08 +00:00
|
|
|
//
|
|
|
|
// Check input samplesheet and get read channels
|
|
|
|
//
|
|
|
|
|
|
|
|
|
2022-03-25 13:58:06 +00:00
|
|
|
include { SHORTREAD_FASTP } from './shortread_fastp'
|
|
|
|
include { FASTQC as FASTQC_PROCESSED } from '../../modules/nf-core/modules/fastqc/main'
|
2022-02-19 11:36:08 +00:00
|
|
|
|
2022-03-21 17:29:47 +00:00
|
|
|
workflow SHORTREAD_PREPROCESSING {
|
2022-02-19 11:36:08 +00:00
|
|
|
take:
|
|
|
|
reads // file: /path/to/samplesheet.csv
|
|
|
|
|
|
|
|
main:
|
2022-03-25 13:58:06 +00:00
|
|
|
ch_versions = Channel.empty()
|
|
|
|
ch_multiqc_files = Channel.empty()
|
2022-02-19 11:36:08 +00:00
|
|
|
|
|
|
|
//
|
|
|
|
// STEP: Read clipping and merging
|
|
|
|
//
|
|
|
|
// TODO give option to clip only and retain pairs
|
|
|
|
// TODO give option to retain singletons (probably fastp option likely)
|
|
|
|
// TODO move to subworkflow
|
|
|
|
|
2022-03-25 13:58:06 +00:00
|
|
|
if ( params.shortread_clipmerge_tool == "fastp" ) {
|
|
|
|
ch_processed_reads = SHORTREAD_FASTP ( reads ).reads
|
|
|
|
ch_versions = ch_versions.mix( SHORTREAD_FASTP.out.versions )
|
|
|
|
ch_multiqc_files = ch_multiqc_files.mix( SHORTREAD_FASTP.out.mqc )
|
2022-02-19 11:36:08 +00:00
|
|
|
} else {
|
|
|
|
ch_processed_reads = reads
|
|
|
|
}
|
|
|
|
|
2022-03-25 13:58:06 +00:00
|
|
|
//FASTQC_PROCESSED ( ch_processed_reads )
|
|
|
|
//ch_versions = ch_versions.mix( FASTQC_PROCESSED.out.versions )
|
|
|
|
//ch_multiqc_files = ch_multiqc_files.mix( FASTQC_PROCESSED.out.zip.collect{it[1]} )
|
2022-02-19 11:36:08 +00:00
|
|
|
|
|
|
|
emit:
|
2022-03-25 13:58:06 +00:00
|
|
|
// TODO: problem, this is being exported as a multi-channel output? This is why FASTQC is broken
|
2022-02-19 11:36:08 +00:00
|
|
|
reads = ch_processed_reads // channel: [ val(meta), [ reads ] ]
|
|
|
|
versions = ch_versions // channel: [ versions.yml ]
|
|
|
|
mqc = ch_multiqc_files
|
|
|
|
}
|
|
|
|
|