Commit 0b02d753 by Eric Coissac

patch a bug in solexapairend

parent 8177aaef
......@@ -153,7 +153,6 @@ def buildJoinedSequence(ali,options):
d = ali[0].getRoot()
r = ali[1].getRoot()
junction = options.junction
r=r.complement()
......@@ -168,9 +167,7 @@ def buildJoinedSequence(ali,options):
quality = d.quality
else:
quality = [10**-4] * len(d)
quality.extend([nqual] * junction)
if hasattr(r, 'quality'):
quality.extend(r.quality)
else:
......
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