Relax template alignment positions to account for optimize positions

wolongac · Sep 5, 2020 · 2c12c35 · 2c12c35
1 parent c1329f7
commit 2c12c35
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/utgcns/unitigConsensus.C b/src/utgcns/unitigConsensus.C
@@ -480,7 +480,7 @@ unitigConsensus::generateTemplateStitch(void) {
     bool   moreToExtend  = (readEnd < readLen);
 
 
-    int32 maxDifference = min(1000, (int32)ceil(0.1*olapLen));
+    int32 maxDifference = min(2500, (int32)ceil(0.30*olapLen));
     //  Reset if the edit distance is waay more than our error rate allows or it's very short and we haven't topped out on error.  This seems to be a quirk with
     //  edlib when aligning to N's - I got startLocation = endLocation = 0 and editDistance = alignmentLength.
     if ((double)result.editDistance / result.alignmentLength > bandErrRate ||