From a5103d7112b7216b3ba7338185b32637d527fcd6 Mon Sep 17 00:00:00 2001 From: Dmitry Antipov Date: Tue, 13 Aug 2024 17:34:14 -0400 Subject: [PATCH] experimentally allowing hap-dip for multimappers only --- src/scripts/scaffolding/scaffold_graph.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/scripts/scaffolding/scaffold_graph.py b/src/scripts/scaffolding/scaffold_graph.py index 146c99b7..818141e4 100755 --- a/src/scripts/scaffolding/scaffold_graph.py +++ b/src/scripts/scaffolding/scaffold_graph.py @@ -424,11 +424,11 @@ def findExtension(self, cur_path_id, type): best_ind = -1 second_best_ind = -1 for i in range (0, len(local_scores)): - if not self.forbiddenPair(cur_path_id, local_scores[i][0]): + if (type != "unique_weight") or not self.forbiddenPair(cur_path_id, local_scores[i][0]): best_ind = i break for j in range (best_ind+1, len(local_scores)): - if not self.forbiddenPair(cur_path_id, local_scores[j][0]): + if (type != "unique_weight") or not self.forbiddenPair(cur_path_id, local_scores[j][0]): second_best_ind = j break if len(local_scores) == 0: