Skip to content

Commit

Permalink
fix merge
Browse files Browse the repository at this point in the history
  • Loading branch information
skoren committed Jul 25, 2024
1 parent 0cc4322 commit e225af6
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion src/scripts/resolve_triplets_kmerify.py
Original file line number Diff line number Diff line change
Expand Up @@ -218,14 +218,16 @@ def get_valid_triplets(node, edges, paths_crossing, min_edge_support, min_covera

def resolve_hairpins(nodelength, nodes, paths_crossing, node_seqs, node_lens, edges, maybe_resolvable, min_edge_support, min_coverage, removable_nodes):
hairpins = set()
unresolvable_hairpins = set()

for node in gf.iterate_deterministic(nodes):
if ">" + node not in edges: continue
if "<" + node not in edges: continue
if len(edges[">" + node]) >= 2 and "<" + node in edges[">" + node]:
unresolvable_hairpins.add(node)
if len(edges["<" + node]) >= 2 and ">" + node in edges["<" + node]:
unresolvable_hairpins.add(node)
if len(edges[">" + node]) == 1 and gf.getone(edges[">" + node]) == "<" + node:
if len(edges[">" + node]) == 1 and gf.getone(edges[">" + node]) == "<" + node:
hairpins.add(">" + node)
if len(edges["<" + node]) == 1 and gf.getone(edges["<" + node]) == ">" + node:
hairpins.add("<" + node)
Expand Down

0 comments on commit e225af6

Please sign in to comment.