Remove limit of AS_MAX_READLEN on utgcns

marbl · Mar 10, 2022 · 7041856 · 7041856
1 parent 001fa55
commit 7041856
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/utgcns/unitigConsensus.C b/src/utgcns/unitigConsensus.C
@@ -298,7 +298,7 @@ unitigConsensus::generateTemplateStitch(void) {
   char        *fragment = seq->getBases();
   uint32       readLen  = seq->length();
 
-  uint32       tigmax = AS_MAX_READLEN;  //  Must be at least AS_MAX_READLEN, else resizeArray() could fail
+  uint32       tigmax = std::max(readLen, AS_MAX_READLEN);  //  Must be at least AS_MAX_READLEN, else resizeArray() could fail
   uint32       tiglen = 0;
   char        *tigseq = NULL;