[Bioperl-l] forcing alphabet in Bio::AlignIO

Jason Stajich jason at bioperl.org
Fri Jan 22 19:18:30 UTC 2010


Done, as of r16739. Look forward to the refactor work too.

-jason
On Jan 22, 2010, at 5:34 AM, Chris Fields wrote:

> Sounds good to me.  The warnings are a bit too tight on this module  
> anyway.
>
> I still think we have plans towards refactoring some of this, not  
> sure how far along they are:
>
> http://www.bioperl.org/wiki/Align_Refactor
>
> chris
>
> On Jan 22, 2010, at 12:17 AM, Jason Stajich wrote:
>
>> I'm considering putting in allowable initialization parameter (and  
>> get/set) for Bio::AlignIO that would allow setting of the  
>> alphabet.  This is then passed to Bio::LocatableSeq creation so  
>> that _guess_alphabet isn't called. This will allow removal of  
>> warnings about empty sequences because _guess_alphabet won't be  
>> called on a sequence if we have explictly set the alphabet.
>>
>> This worked great on my local install and tests pass.  Any  
>> objections or concerns?
>>
>> basically it means when you make an AlignIO you can specify the  
>> alphabet i.e.
>>
>> my $in = Bio::AlignIO->new(-format => 'fasta', -alphabet => 'dna', - 
>> file => 'genome.fasaln');
>>
>> I have some alignments with empty sequences and I think turning off  
>> the warnings is appropriate where I force the alphabet choice. It  
>> should also have a very modest speedup benefit too.
>>
>> -jason
>> --
>> Jason Stajich
>> jason.stajich at gmail.com
>> jason at bioperl.org
>> http://fungalgenomes.org/
>> http://twitter.com/hyphaltip
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>

--
Jason Stajich
jason.stajich at gmail.com
jason at bioperl.org
http://fungalgenomes.org/
http://twitter.com/hyphaltip




More information about the Bioperl-l mailing list