[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