Skip to content
Snippets Groups Projects
Verified Commit 1903be3c authored by Etienne MORICE's avatar Etienne MORICE
Browse files

Patch missing id in multiple alignment

parent 578d2d27
No related branches found
No related tags found
No related merge requests found
...@@ -136,7 +136,7 @@ def join_alignments(aleft, aright, _mscore=_mscore, gap=-8): ...@@ -136,7 +136,7 @@ def join_alignments(aleft, aright, _mscore=_mscore, gap=-8):
aseqs[len(aleft)+jseq] += aright[jseq][nextj] aseqs[len(aleft)+jseq] += aright[jseq][nextj]
i,j = nexti,nextj i,j = nexti,nextj
ret = [ SeqRecord(s) for s in aleft ] + [ SeqRecord(s) for s in aright ] ret = [ SeqRecord(id=s.id, seq="") for s in aleft ] + [ SeqRecord(id=s.id, seq="") for s in aright ]
for rec, seq in zip(ret, aseqs): for rec, seq in zip(ret, aseqs):
rec.seq = Seq(seq[::-1]) rec.seq = Seq(seq[::-1])
return ret return ret
......
...@@ -15,7 +15,7 @@ class MultipleAlignTestCase(unittest.TestCase): ...@@ -15,7 +15,7 @@ class MultipleAlignTestCase(unittest.TestCase):
def test_align_from_tree(self): def test_align_from_tree(self):
aseqs = [ s.seq for s in aseqs = [ s.seq for s in
alignementseq_multiple.multiple_align_from_tree( alignementseq_multiple.multiple_align_from_tree(
(("CHAT", "CAT"), "HER") ((SeqRecord("CHAT"), SeqRecord("CAT")), SeqRecord("HER"))
)] )]
aseqs.sort() aseqs.sort()
...@@ -28,7 +28,7 @@ class MultipleAlignTestCase(unittest.TestCase): ...@@ -28,7 +28,7 @@ class MultipleAlignTestCase(unittest.TestCase):
import tree_generator import tree_generator
import alignementseq import alignementseq
seqs = [ SeqRecord(s) for s in ("CHAT", "CAT", "HER") ] seqs = [ SeqRecord(seq=s) for s in ("CHAT", "CAT", "HER") ]
tree = tree_generator.tree_build(seqs, alignementseq.vec_align) tree = tree_generator.tree_build(seqs, alignementseq.vec_align)
aseqs = alignementseq_multiple.multiple_align_from_linkage_matrix( aseqs = alignementseq_multiple.multiple_align_from_linkage_matrix(
seqs, seqs,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment