Commit d086f56c by Eric Coissac

Patches two bugs related to circular contigs

- Patch a bug leading to delete a base at the end of circular contigs
- Patch a bug leading to give different ids to reverse-complement
circular contigs
parent 1d92aa97
......@@ -469,10 +469,14 @@ cdef class StemIterator:
if circle:
minabs = min(abs(i) for i in cc)
print(minabs,file=sys.stderr)
if minabs in cc:
first = minabs
print("in",file=sys.stderr)
else:
first = sons(graph,-minabs).pop()
# first = sons(graph,-minabs).pop()
first = -minabs
print("out",file=sys.stderr)
nn = first
path=[first]
lsequence=[]
......@@ -484,6 +488,9 @@ cdef class StemIterator:
son=sons(graph,nn).pop()
last = first
path.append(first)
lsequence.append(graph.getEdgeAttr(nn,first)['ext'])
sequence=b''.join(lsequence)
stem = buildstem(self._assembler,
......
Markdown is supported
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