diff --git a/subworkflows/local/hicstuff_sub.nf b/subworkflows/local/hicstuff_sub.nf
index 09cd683e3ee376045ce662d68ae4628721b681df..60f1f7903d757471577ae13fa582cdd0182fc135 100644
--- a/subworkflows/local/hicstuff_sub.nf
+++ b/subworkflows/local/hicstuff_sub.nf
@@ -3,6 +3,7 @@ include { BOWTIE2_ALIGNMENT } from '../../modules/local/hicstuff/align_bowtie2'
 include { FRAGMENT_ENZYME } from '../../modules/local/hicstuff/fragment_enzyme'
 include { BAM2PAIRS } from '../../modules/local/hicstuff/bam2pairs'
 include { BUILD_MATRIX } from '../../modules/local/hicstuff/build_matrix'
+include { BUILD_MATRIX_COOL } from '../../modules/local/hicstuff/build_matrix_cool'
 
 // Paired-end to Single-end
 def pairToSingle(row, mates) {
@@ -65,4 +66,9 @@ workflow HICSTUFF_SUB {
         BAM2PAIRS.out.idx_pairs,
         FRAGMENT_ENZYME.out.fragments_list
     )
+
+    BUILD_MATRIX_COOL(
+        BAM2PAIRS.out.idx_pairs,
+        FRAGMENT_ENZYME.out.fragments_list
+    )
 }