[Bioperl-l] about the iupac
Magic Fang
fangl at genomics.org.cn
Thu Mar 6 09:10:09 EST 2003
why my iupac module can not translate seq with ambigious base to the unambigious? the codes are:
#!/usr/bin/perl
use Bio::Seq;
use Bio::Tools::IUPAC;
my $ambiseq = new Bio::Seq (-seq => 'AYTCGUTGR', -type => 'Dna');
my $stream = new Bio::Tools::IUPAC(-seq => $ambiseq);
print $ambiseq->seq(), "\n";
while ($uniqueseq = $stream->next_seq()) {
# process the unique Seq object.
print $uniqueseq->seq, "\n";
}
the out put are the same seq.
does this module will translate the seq to unambigious seq? doubt:-)
More information about the Bioperl-l
mailing list