@@ -12,7 +12,7 @@ include { STARFUSION } from '../modules/local/starfusion/de
1212include { STRINGTIE_WORKFLOW } from ' ../subworkflows/local/stringtie_workflow/main'
1313include { FUSIONCATCHER_WORKFLOW } from ' ../subworkflows/local/fusioncatcher_workflow'
1414include { FUSIONINSPECTOR_WORKFLOW } from ' ../subworkflows/local/fusioninspector_workflow'
15- include { FUSIONREPORT_WORKFLOW } from ' ../subworkflows/local/fusionreport_workflow '
15+ include { FUSIONREPORT_DETECT } from ' ../modules/nf-core/fusionreport/detect/main '
1616include { FASTQC } from ' ../modules/nf-core/fastqc/main'
1717include { MULTIQC } from ' ../modules/nf-core/multiqc/main'
1818include { STAR_ALIGN } from ' ../modules/nf-core/star/align/main'
@@ -308,21 +308,26 @@ workflow RNAFUSION {
308308 def ch_fusion_list_filtered = Channel . empty()
309309 def ch_fusionreport_report = Channel . empty()
310310 def ch_fusionreport_csv = Channel . empty()
311- if (! params. skip_vis && tools. contains(" fusionreport" )) {
311+ if (! params. skip_vis && tools. contains(" fusionreport" )) {
312312 if (! fusions_created) {
313313 error(" Could not find any fusion files. Please generate some with `--tools arriba`, `--tools starfusion` and/or `--tools fusioncatcher`" )
314314 }
315- FUSIONREPORT_WORKFLOW (
315+
316+ def ch_fusions = ch_arriba_fusions
317+ .join(ch_starfusion_fusions, failOnMismatch :true , failOnDuplicate :true )
318+ .join(ch_fusioncatcher_fusions, failOnMismatch :true , failOnDuplicate :true )
319+
320+ FUSIONREPORT_DETECT (
321+ ch_fusions,
316322 BUILD_REFERENCES . out. fusionreport_ref,
317- ch_arriba_fusions,
318- ch_starfusion_fusions,
319- ch_fusioncatcher_fusions
323+ params. tools_cutoff
320324 )
321- ch_versions = ch_versions. mix(FUSIONREPORT_WORKFLOW . out. versions)
322- ch_fusion_list = FUSIONREPORT_WORKFLOW . out. fusion_list
323- ch_fusion_list_filtered = FUSIONREPORT_WORKFLOW . out. fusion_list_filtered
324- ch_fusionreport_report = FUSIONREPORT_WORKFLOW . out. report
325- ch_fusionreport_csv = FUSIONREPORT_WORKFLOW . out. csv
325+
326+ ch_versions = ch_versions. mix(FUSIONREPORT_DETECT . out. versions)
327+ ch_fusion_list = FUSIONREPORT_DETECT . out. fusion_list
328+ ch_fusion_list_filtered = FUSIONREPORT_DETECT . out. fusion_list_filtered
329+ ch_fusionreport_report = FUSIONREPORT_DETECT . out. report
330+ ch_fusionreport_csv = FUSIONREPORT_DETECT . out. csv
326331 } else if (params. fusioninspector_fusions) {
327332 def input_fusions = file(params. fusioninspector_fusions, checkIfExists :true )
328333 ch_fusion_list = ch_reads. map { it -> [ it[0 ], input_fusions ] }
0 commit comments