Skip to content

Commit

Permalink
skip nodes completely composed of a gapfill
Browse files Browse the repository at this point in the history
  • Loading branch information
skoren committed Jan 31, 2025
1 parent 9a78f76 commit 2c7336d
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions src/scripts/get_layout_from_mbg.py
Original file line number Diff line number Diff line change
Expand Up @@ -263,9 +263,11 @@ def get_exact_match_length(clusters):
pieceid = pieceid + 1
pathname = f"piece{pieceid:06d}"

contig_pieces[fullname].append(pathname)

(path, overlaps) = get_leafs(re.findall(r"[<>][^<>]+", pp), node_mapping, edge_overlaps, raw_node_lens)
# skip a path if the only thing in it is a gapfill
if len(path) == 1 and path[0][0][1:4] == "gap":
continue
contig_pieces[fullname].append(pathname)

contig_nodeseqs[pathname] = path
contig_nodeoverlaps[pathname] = overlaps
Expand Down

0 comments on commit 2c7336d

Please sign in to comment.