[Biopython-dev] [Bug 2531] Nexus and fasta parsers have a problem with identical taxa names

Frank Kauff fkauff at biologie.uni-kl.de
Mon Jun 30 17:10:15 UTC 2008



bugzilla-daemon at portal.open-bio.org wrote:
>
>
> In the Bio.SeqIO code that calls Bio.Nexus, I hadn't realized that Bio.Nexus
> kept the un-edited taxon names around.  It is this list of the non-unique
> original identifiers that Bio.SeqIO was using, which explains why you end up
> with two copies of HI99.Line5.
>
> Sorry Frank - I was pointing fingers when it was my own bug after all!
>
>
> Looking back, the reason I was using the original_taxon_order list was I wanted
> to get the sequences in their original order.  I see now that I can't use the
> elements in this list as keys to the matrix because the matrix keys are the
> modified taxon names.
>
> Is there any way to get the modified taxon names in the original order?  Other
> than looping over original_taxon_order and repeating your naming algorithm?
>   
Actually -this *IS* a bug. All fingers were pointing correctly... 
Original_taxon labels was just kept just for compatibility, and is the 
same as taxlabels. Taxlabels is supposed to have the unique identifiers 
- it just doesn't work correctly with non-unique ids in interleaved data 
sets.
Fix following soon

Frank



More information about the Biopython-dev mailing list