Commit 1f1528ee authored by Martin Larralde's avatar Martin Larralde
Browse files

Fix `OrfFinder.find_genes` not rebuilding node scores of the maximum phase in meta mode

parent c883d994
......@@ -4385,7 +4385,12 @@ cdef class OrfFinder:
)
nodes._sort()
nodes._score(sequence, tinf, closed=self.closed, is_meta=True)
# rescore nodes (this is not done in the Prodigal code, but if
# we want to provide accurate scores in the `Node` objects we
# have to do it)
nodes._record_overlapping_starts(tinf, True, self.max_overlap)
nodes._dynamic_programming(tinf, scorer, final=True)
# return the max phase on success
return max_phase
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment