[Bioperl-l] SeqIO fails on masked sequences

Marc Logghe Marc.Logghe at devgen.com
Wed Dec 15 05:27:20 EST 2004

> my $seq_in  = Bio::SeqIO->new(-format=>$format, -file => $infile);
> my $seq_out = Bio::SeqIO->new(-format=>$format, -file => ">$outfile");
> while (my $seq = $seq_in->next_seq()) {

Guess you can do it by setting the alphabet explicitely:
$seq_in->alphabet('dna'); # or 'rna' or 'protein'

You probably also can pass it directly with the constructor, but did not find that immediately.
Indirectly, you can do it also by setting the alphabet for the factory object and passing the factory object with the Bio::SeqIO constructor.


More information about the Bioperl-l mailing list