From 1260bf9f5c42c5b863d68fd1077f188f4df9d9e4 Mon Sep 17 00:00:00 2001
From: Tim Dunn <me.timd1@gmail.com>
Date: Tue, 2 Jul 2024 19:18:21 +0000
Subject: [PATCH] fix: min_intron_length_ initialization in JunctionsExtractor

 - fixes #188, min_intron_length_ was accidentally set from min_anchor_len_
---
 src/junctions/junctions_extractor.h | 19 +++++++++++++++++--
 1 file changed, 17 insertions(+), 2 deletions(-)

diff --git a/src/junctions/junctions_extractor.h b/src/junctions/junctions_extractor.h
index a95d199..0b14404 100644
--- a/src/junctions/junctions_extractor.h
+++ b/src/junctions/junctions_extractor.h
@@ -196,8 +196,23 @@ class JunctionsExtractor {
             region_ = ".";
             ref_ = "NA";
         }
-        JunctionsExtractor(string bam1, string region1, int strandness1, string strand_tag1, uint32_t min_anchor_length1, uint32_t min_intron_length1, uint32_t max_intron_length1, string ref1) : 
-            bam_(bam1), region_(region1), strandness_(strandness1), strand_tag_(strand_tag1), min_anchor_length_(min_anchor_length1), min_intron_length_(min_anchor_length1), max_intron_length_(max_intron_length1), ref_(ref1){
+        JunctionsExtractor(
+                string bam1, 
+                string region1, 
+                int strandness1, 
+                string strand_tag1, 
+                uint32_t min_anchor_length1, 
+                uint32_t min_intron_length1, 
+                uint32_t max_intron_length1, 
+                string ref1) : 
+            bam_(bam1), 
+            region_(region1), 
+            strandness_(strandness1), 
+            strand_tag_(strand_tag1), 
+            min_anchor_length_(min_anchor_length1), 
+            min_intron_length_(min_intron_length1), 
+            max_intron_length_(max_intron_length1), 
+            ref_(ref1) {
             junctions_sorted_ = false;
             output_file_ = "NA";
             output_barcodes_file_ = "NA";