[Bioperl-l] SimpleAlign - get_seq_by_id
rvos at interchange.ubc.ca
Fri Oct 24 06:40:25 UTC 2008
I would be very hesitant to introduce array contexts - it just yields
subtle bugs and maintenance issues. Perhaps a separate
get_all_seqs_by_id would be better?
On Fri, Oct 24, 2008 at 2:32 PM, Heikki Lehvaslaiho <heikki at sanbi.ac.za> wrote:
> The main reason it has not been Bio::SeqAlign is that sequence ID not
> necessarily a unique identifier in a MSA. Multiple regions of the sequence
> defined by one ID can be in one.
> The current code returns only the more or less randomly selected first
> Bio::LocatebleSeqI object with that ID. Should we make it context sensitive
> and return an array of sequences in array context?
> That brings up an other question: After the change, the get_seq_by_id() will
> behave differently from all other instances of that method, so should it be
> renamed to reflect that?
> On Thursday 23 October 2008 21:29:20 Jason Stajich wrote:
>> I added get_seq_by_id to Bio::SimpleAlign to allow retrieval of a
>> particular sequence from the alignment by ID. Not sure why this didn't
>> exist before.
>> Jason Stajich
>> jason at bioperl.org
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
> ______ _/ _/_____________________________________________________
> _/ _/
> _/ _/ _/ Heikki Lehvaslaiho heikki at_sanbi _ac _za
> _/_/_/_/_/ Senior Scientist skype: heikki_lehvaslaiho
> _/ _/ _/ SANBI, South African National Bioinformatics Institute
> _/ _/ _/ University of Western Cape, South Africa
> _/ Phone: +27 21 959 2096 FAX: +27 21 959 2512
> ___ _/_/_/_/_/________________________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
Dr. Rutger A. Vos
Department of zoology
University of British Columbia
More information about the Bioperl-l