@@ -443,7 +443,7 @@ ${summary.collect { k,v -> " <dt>$k</dt><dd><samp>${v ?: '<span style
443443
444444
445445/*
446- * Create BWA indices if they are not present
446+ * PREPROCESSING - Create BWA indices if they are not present
447447*/
448448
449449if (! params. bwa_index && ! params. fasta. isEmpty() && (params. aligner == ' bwa' || params. bwamem)){
@@ -528,7 +528,7 @@ process makeSeqDict {
528528}
529529
530530/*
531- * Convert BAM to FastQ if BAM input is specified instead of FastQ file(s)
531+ * PREPROCESSING - Convert BAM to FastQ if BAM input is specified instead of FastQ file(s)
532532*
533533*/
534534
@@ -554,7 +554,7 @@ process convertBam {
554554
555555
556556/*
557- * STEP 1 - FastQC
557+ * STEP 1a - FastQC
558558 */
559559process fastqc {
560560 tag " $name "
@@ -578,7 +578,7 @@ process fastqc {
578578}
579579
580580
581- /* STEP 2.0 - FastP
581+ /* STEP 1b - FastP
582582* Optional poly-G complexity filtering step before read merging/adapter clipping etc
583583* Note: Clipping, Merging, Quality Trimning are turned off here - we leave this to adapter removal itself!
584584*/
@@ -676,7 +676,7 @@ process adapter_removal {
676676
677677
678678/*
679- * STEP 2.1 - FastQC after clipping/merging (if applied!)
679+ * STEP 2b - FastQC after clipping/merging (if applied!)
680680*/
681681process fastqc_after_clipping {
682682 tag " ${ name} "
@@ -698,7 +698,7 @@ process fastqc_after_clipping {
698698}
699699
700700/*
701- Step 3: Mapping with BWA, SAM to BAM, Sort BAM
701+ Step 3a - Mapping with BWA, SAM to BAM, Sort BAM
702702*/
703703
704704process bwa {
@@ -848,7 +848,7 @@ process bwamem {
848848}
849849
850850/*
851- * Step 4 - flagstat
851+ * Step 3b - flagstat
852852*/
853853
854854process samtools_flagstat {
@@ -870,7 +870,7 @@ process samtools_flagstat {
870870
871871
872872/*
873- * Step 5: Keep unmapped/remove unmapped reads
873+ * Step 4a - Keep unmapped/remove unmapped reads
874874*/
875875
876876process samtools_filter {
@@ -960,7 +960,7 @@ process strip_input_fastq {
960960}
961961
962962/*
963- * Step 5b : Keep unmapped/remove unmapped reads flagstat
963+ * Step 4b : Keep unmapped/remove unmapped reads flagstat
964964*/
965965
966966
@@ -983,7 +983,7 @@ process samtools_flagstat_after_filter {
983983
984984
985985/*
986- Step 6 : DeDup / MarkDuplicates
986+ Step 5a : DeDup / MarkDuplicates
987987*/
988988
989989process dedup{
@@ -1026,7 +1026,7 @@ process dedup{
10261026}
10271027
10281028/*
1029- Step 5.1 : Preseq
1029+ Step 6 : Preseq
10301030*/
10311031
10321032process preseq {
@@ -1056,7 +1056,7 @@ process preseq {
10561056}
10571057
10581058/*
1059- Step 5.2 : DMG Assessment
1059+ Step 7a : DMG Assessment
10601060*/
10611061
10621062process damageprofiler {
@@ -1084,7 +1084,7 @@ process damageprofiler {
10841084}
10851085
10861086/*
1087- Step 5.3 : Qualimap
1087+ Step 8 : Qualimap
10881088*/
10891089
10901090process qualimap {
@@ -1112,7 +1112,7 @@ process qualimap {
11121112
11131113
11141114/*
1115- Step 6 : MarkDuplicates
1115+ Step 5b : MarkDuplicates
11161116 */
11171117
11181118process markDup{
@@ -1152,6 +1152,10 @@ if(!params.run_pmdtools){
11521152 ch_dedup_for_pmdtools. close()
11531153}
11541154
1155+ /*
1156+ Step 9: PMDtools
1157+ */
1158+
11551159process pmdtools {
11561160 tag " ${ bam.baseName} "
11571161 publishDir " ${ params.outdir} /pmdtools" , mode: ' copy'
@@ -1184,7 +1188,7 @@ process pmdtools {
11841188}
11851189
11861190/*
1187- * Optional BAM Trimming step using bamUtils
1191+ * Step 10 - BAM Trimming step using bamUtils
11881192* Can be used for UDGhalf protocols to clip off -n bases of each read
11891193*/
11901194
@@ -1239,7 +1243,7 @@ Downstream VCF tools:
12391243
12401244
12411245/*
1242- * STEP 3 - Output Description HTML
1246+ * Step 11a - Output Description HTML
12431247 */
12441248process output_documentation {
12451249 publishDir " ${ params.outdir} /Documentation" , mode: ' copy'
@@ -1258,7 +1262,7 @@ process output_documentation {
12581262
12591263
12601264/*
1261- * Parse software version numbers
1265+ * Step 11b - Parse software version numbers
12621266 */
12631267process get_software_versions {
12641268
@@ -1293,7 +1297,7 @@ process get_software_versions {
12931297
12941298
12951299/*
1296- * STEP 2 - MultiQC
1300+ * Step 11c - MultiQC
12971301 */
12981302process multiqc {
12991303 publishDir " ${ params.outdir} /MultiQC" , mode: ' copy'
@@ -1331,7 +1335,7 @@ process multiqc {
13311335
13321336
13331337/*
1334- * Completion e-mail notification
1338+ * Step 11d - Completion e-mail notification
13351339 */
13361340workflow. onComplete {
13371341
0 commit comments