From 66f59323b52eb3551b51d6fe7fd39a2e2587abae Mon Sep 17 00:00:00 2001 From: szhan Date: Fri, 8 May 2026 15:03:38 +0100 Subject: [PATCH] Skip multi-parent recombinants --- sc2ts/debug.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/sc2ts/debug.py b/sc2ts/debug.py index 68545225..468f20d0 100644 --- a/sc2ts/debug.py +++ b/sc2ts/debug.py @@ -913,7 +913,9 @@ def node_info(node, label): breakpoint_intervals.append(node_md["breakpoint_intervals"]) # Only deal with 2 parents recombs for now. - assert self.nodes_num_parents[u] == 2 + if self.nodes_num_parents[u] > 2: + logger.warning(f"RE: {u} has more than two parents. Skipping.") + continue # assert len(set(hmm_matches)) == 1 # assert len(set(breakpoint_intervals)) == 1 hmm_match = hmm_matches[0]