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

redmine at redmine.open-bio.org redmine at redmine.open-bio.org
Tue Nov 8 10:17:00 UTC 2011

Issue #3312 has been reported by gahoo lee.

Bug #3312: Failing to parse fasta-m10 format generated by lalign36

Author: gahoo lee
Status: New
Priority: Normal
Target version: 

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 this email was added to the New Issue Alert plugin

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

More information about the Biopython-dev mailing list