Commit f61e648a by Eric Coissac

Fix a bug in fillgap created by the change of API

parent 53852200
......@@ -490,7 +490,7 @@ def fillGaps(self,minlink=5,
for e2 in range(nei):
connected = et[e1][1]==ei[e2][0]
if not connected:
linkedby,ml,sl,delta = pairEndedConnected(self,assgraph,etid[e1],eiid[e2],back)
linkedby,ml,sl,pdelta = pairEndedConnected(self,assgraph,etid[e1],eiid[e2],back) # @UnusedVariable
if linkedby >= minlink and abs(etid[e1]) <= abs(eiid[e2]):
extended.add(etid[e1])
......@@ -1840,7 +1840,7 @@ def fillGaps2(self,minlink=5,
for e2 in range(nei):
connected = et[e1][1]==ei[e2][0]
if not connected:
linkedby,ml,sl,delta = pairEndedConnected(self,assgraph,etid[e1],eiid[e2],back)
linkedby,ml,sl,pdelta = pairEndedConnected(self,assgraph,etid[e1],eiid[e2],back) # @UnusedVariable
if linkedby >= minlink and abs(etid[e1]) <= abs(eiid[e2]):
extended.add(etid[e1])
......@@ -2169,7 +2169,7 @@ def path2sam(self,assgraph,path,identifier="contig",minlink=10,nlength=20,back=2
s2 = alledges[path[0]]
sid1=assgraph.getEdgeAttr(*s1)['stemid']
sid2=assgraph.getEdgeAttr(*s2)['stemid']
connected,ml,sl,delta = pairEndedConnected(self,assgraph,sid1,sid2,back)
connected,ml,sl,delta = pairEndedConnected(self,assgraph,sid1,sid2,back) # @UnusedVariable
if s1[1]==s2[0]:
print("Path is circular",file=sys.stderr)
......
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