[Biopython-dev] [Biopython - Bug #3312] Failing to parse fasta-m10 format generated by lalign36

[redmine at redmine.open-bio.org]

Issue #3312 has been updated by gahoo lee.

File 3Seqs.zip added

Well, I'm not on the FASTA mailing list. In fact I found a small bug in mshowalign2.c which a colon is missing on line 616, just don't know how to join the mailing list.
Here's the FASTA output with 3 sequences alignment, I hope these file would help. The odd lines changed in the output.
----------------------------------------
Bug #3312: Failing to parse fasta-m10 format generated by lalign36
https://redmine.open-bio.org/issues/3312

Author: gahoo lee
Status: New
Priority: Normal
Assignee: Biopython Dev Mailing List
Category: 
Target version: 
URL: 


When I parse an alignment created by lalign which is included in FASTA36, I got errors. We got two sequences in each fasta file now, but if one sequence each, there's no error. Here are the codes and error.

@lalign36 -m 10 at.fasta os.fasta >test.aln@

@from Bio import AlignIO
handle = open('test.aln')
for a in AlignIO.parse(handle, "fasta-m10"):
    assert len(a) == 2, "Should be pairwise!"
    print "Alignment length %i" % a.get_alignment_length()
    for record in a:
        print record.seq, record.name, record.id
@

@Traceback (most recent call last):
  File "R:\Untitled 4.py", line 5, in <module>
    for a in AlignIO.parse(handle, "fasta-m10"):
  File "D:\Program Files\Python\lib\site-packages\Bio\AlignIO\__init__.py", line 371, in parse
    for a in i:
  File "D:\Program Files\Python\lib\site-packages\Bio\AlignIO\FastaIO.py", line 242, in FastaM10Iterator
    yield build_hsp()
  File "D:\Program Files\Python\lib\site-packages\Bio\AlignIO\FastaIO.py", line 106, in build_hsp
    assert query_tags, query_tags
AssertionError: {}@


-- 
You have received this notification because you have either subscribed to it, or are involved in it.
To change your notification preferences, please click here and login: http://redmine.open-bio.org