Skip to content

Commit 286ec9f

Browse files
committed
fix bug in last commit
1 parent 3093eff commit 286ec9f

File tree

1 file changed

+4
-1
lines changed

1 file changed

+4
-1
lines changed

src/rgfa.cpp

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -178,7 +178,10 @@ void RGFACover::compute(const PathHandleGraph* graph,
178178
for (int64_t j = 0; j < rgfa_intervals_vector[t].size(); ++j) {
179179
// the true flag at the end disables the overlap check. since they were computed
180180
// in separate threads, snarls can overlap by a single node
181-
add_interval(this->rgfa_intervals, this->node_to_interval, rgfa_intervals_vector[t][j], true);
181+
const pair<step_handle_t, step_handle_t>& interval = rgfa_intervals_vector[t][j];
182+
if (interval.first != graph->path_end(graph->get_path_handle_of_step(interval.first))) {
183+
add_interval(this->rgfa_intervals, this->node_to_interval, rgfa_intervals_vector[t][j], true);
184+
}
182185
}
183186
rgfa_intervals_vector[t].clear();
184187
node_to_interval_vector[t].clear();

0 commit comments

Comments
 (0)