[Bioperl-l] about the iupac - It works...

[hanfos]

The result are:

Original sequence: AYTCGUTGR
ACTCGUTGA
ATTCGUTGA
ACTCGUTGG
ATTCGUTGG

I am sorry why you could not do that!

Sincerely yours,

Sang Chul Choi

----- Original Message ----- 
From: "Magic Fang" <fangl at genomics.org.cn>
To: <bioperl-l at bioperl.org>
Sent: Thursday, March 06, 2003 10:10 AM
Subject: [Bioperl-l] about the iupac


> 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:-)
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>