Sep 10, 2024: Correct gapLen

vgl-hub · Sep 11, 2024 · 08216bc · 08216bc
1 parent 9301b2e
commit 08216bc
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/input.cpp b/src/input.cpp
@@ -99,7 +99,7 @@ bool Teloscope::walkPath(InPath* path, std::vector<InSegment*> &inSegments, std:
         }else if (component->componentType == GAP){
 
             auto inGap = find_if(inGaps.begin(), inGaps.end(), [cUId](InGap& obj) {return obj.getuId() == cUId;}); // given a node Uid, find it
-            gapLen += inGap->getDist(component->start - component->end);
+            gapLen = inGap->getDist(component->start - component->end);
 
             absPos += gapLen;
 

diff --git a/src/teloscope.cpp b/src/teloscope.cpp
@@ -175,7 +175,7 @@ std::vector<WindowData> Teloscope::analyzeSegment(std::string &sequence, UserInp
         }
 
         // Recycle the overlapping string sequence
-        currentWindowSize = std::min(userInput.windowSize, static_cast<uint32_t>(sequence.size() - windowStart));
+        currentWindowSize = std::min(userInput.windowSize, static_cast<uint32_t>(sequence.size() - windowStart)); // CHECK
 
         if (currentWindowSize == userInput.windowSize) {
             window = window.substr(userInput.step) + sequence.substr(windowStart + userInput.windowSize - userInput.step, userInput.step);