emit last range, if any

pangenome · Mar 14, 2024 · cdc1733 · cdc1733
1 parent ff88bfc
commit cdc1733
Showing 1 changed file with 16 additions and 0 deletions.
diff --git a/src/subcommand/paths_main.cpp b/src/subcommand/paths_main.cpp
@@ -505,6 +505,22 @@ int main_paths(int argc, char** argv) {
                     }
                 });
             }
+
+            // Emit last non reference range, if any
+            if (end > start && (end - start) >= min_size_in_bp) {
+                if (_show_step_ranges) {
+                    std::string step_range_str = "";
+                    for (auto& step : step_range) {
+                        const handle_t handle = graph.get_handle_of_step(step);
+                        step_range_str += std::to_string(graph.get_id(handle)) + (graph.get_is_reverse(handle) ? "-" : "+") + ",";
+                    }
+                    #pragma omp critical (cout)
+                        std::cout << graph.get_path_name(path) << "\t" << start << "\t" << end << "\t" << step_range_str.substr(0, step_range_str.size() - 1) << std::endl; // trim the trailing comma from step_range
+                } else {
+                    #pragma omp critical (cout)
+                        std::cout << graph.get_path_name(path) << "\t" << start << "\t" << end << std::endl;
+                }
+            }
         }
     }