2022-04-02 15:02:05 +00:00
|
|
|
//
|
|
|
|
// Process short raw reads with FastP
|
|
|
|
//
|
2022-03-25 13:58:06 +00:00
|
|
|
|
|
|
|
include { FASTP as FASTP_SINGLE } from '../../modules/nf-core/modules/fastp/main'
|
|
|
|
include { FASTP as FASTP_PAIRED } from '../../modules/nf-core/modules/fastp/main'
|
|
|
|
|
|
|
|
workflow SHORTREAD_FASTP {
|
|
|
|
take:
|
2022-03-27 07:30:23 +00:00
|
|
|
reads // [[meta], [reads]]
|
2022-03-25 13:58:06 +00:00
|
|
|
|
|
|
|
main:
|
|
|
|
ch_versions = Channel.empty()
|
|
|
|
ch_multiqc_files = Channel.empty()
|
|
|
|
|
|
|
|
ch_input_for_fastp = reads
|
|
|
|
.branch{
|
|
|
|
single: it[0]['single_end'] == true
|
|
|
|
paired: it[0]['single_end'] == false
|
|
|
|
}
|
|
|
|
|
|
|
|
FASTP_SINGLE ( ch_input_for_fastp.single, false, false )
|
2022-03-27 07:30:23 +00:00
|
|
|
// Last parameter here turns on merging of PE data
|
2022-03-25 13:58:06 +00:00
|
|
|
FASTP_PAIRED ( ch_input_for_fastp.paired, false, params.shortread_clipmerge_mergepairs )
|
|
|
|
|
|
|
|
if ( params.shortread_clipmerge_mergepairs ) {
|
2022-03-28 14:38:10 +00:00
|
|
|
ch_fastp_reads_prepped_pe = FASTP_PAIRED.out.reads_merged
|
|
|
|
.map {
|
|
|
|
meta, reads ->
|
|
|
|
def meta_new = meta.clone()
|
|
|
|
meta_new['single_end'] = 1
|
|
|
|
[ meta_new, reads ]
|
2022-03-25 13:58:06 +00:00
|
|
|
}
|
2022-03-28 14:38:10 +00:00
|
|
|
|
|
|
|
ch_fastp_reads_prepped = ch_fastp_reads_prepped_pe.mix( FASTP_SINGLE.out.reads )
|
|
|
|
|
2022-03-25 13:58:06 +00:00
|
|
|
} else {
|
|
|
|
ch_fastp_reads_prepped = FASTP_PAIRED.out.reads
|
|
|
|
.mix( FASTP_SINGLE.out.reads )
|
|
|
|
}
|
|
|
|
|
|
|
|
ch_versions = ch_versions.mix(FASTP_SINGLE.out.versions.first())
|
|
|
|
ch_versions = ch_versions.mix(FASTP_PAIRED.out.versions.first())
|
|
|
|
|
2022-03-31 13:31:45 +00:00
|
|
|
ch_processed_reads = ch_fastp_reads_prepped
|
2022-03-25 13:58:06 +00:00
|
|
|
|
2022-04-02 15:02:05 +00:00
|
|
|
ch_multiqc_files = ch_multiqc_files.mix( FASTP_SINGLE.out.json )
|
|
|
|
ch_multiqc_files = ch_multiqc_files.mix( FASTP_PAIRED.out.json )
|
2022-03-25 13:58:06 +00:00
|
|
|
|
|
|
|
emit:
|
|
|
|
reads = ch_processed_reads // channel: [ val(meta), [ reads ] ]
|
|
|
|
versions = ch_versions // channel: [ versions.yml ]
|
|
|
|
mqc = ch_multiqc_files
|
|
|
|
}
|
|
|
|
|