Sorted terminal telomeres

vgl-hub · Jan 13, 2025 · e652845 · e652845
1 parent 9129d2e
commit e652845
Showing 1 changed file with 8 additions and 3 deletions.
diff --git a/src/teloscope.cpp b/src/teloscope.cpp
@@ -192,11 +192,11 @@ std::vector<TelomereBlock> Teloscope::filterTerminalBlocks(const std::vector<Tel
     }
 
     // Assign best 'p' and 'q' blocks per path
-    if (has_p1 && has_q1) {
+    if (has_p1 && has_q1) { 
         filteredBlocks.push_back(best_p);
         filteredBlocks.push_back(best_q);
     }
-    else {
+    else { 
         if (has_p1) {
             filteredBlocks.push_back(best_p);
             if (has_p2) {
@@ -211,10 +211,15 @@ std::vector<TelomereBlock> Teloscope::filterTerminalBlocks(const std::vector<Tel
         }
     }
 
+    // Sort by coords
+    std::sort(filteredBlocks.begin(), filteredBlocks.end(),
+            [](const TelomereBlock &a, const TelomereBlock &b) {
+                return a.start < b.start;
+            });
+
     return filteredBlocks;
 }
 
-
 std::vector<TelomereBlock> Teloscope::filterInterstitialBlocks(
             const std::vector<TelomereBlock>& interstitialBlocks,
             const std::vector<TelomereBlock>& terminalBlocks) {