[Bioperl-l] seq to fasta
    stephan rosecker 
    stephan.rosecker at ish.de
       
    Sun May 23 18:17:30 EDT 2004
    
    
  
Hi,
is it possible to get the fasta format from a seqobject-sequence in a 
variable without using filhandles like below ?
use Bio::SeqIO;
   $in  = Bio::SeqIO->new(-file => "inputfilename",
                          -format => 'EMBL');
   $out = Bio::SeqIO->new(-file => ">outputfilename",
                          -format => 'Fasta');
   while ( my $seq = $in->next_seq() ) {$out->write_seq($seq); }
e.g:
I've a seqobject ($seq).
and this is what i want:
my $fastavar = $seq->seq_fastaformat();
print $fastavar;
result:
 >sp|P32234|128U_DROME GTP-binding protein 128UP - Drosophila 
melanogaster (Fruit fly).
MITILEKISAIESEMARTQKNKATSAHLGLLKANVAKLRRELISPKGGGGGTGEAGFEVA
KTGDARVGFVGFPSVGKSTLLSNLAGVYSEVAAYEFTTLTTVPGCIKYKGAKIQLLDLPG
IIEGAKDGKGRGRQVIAVARTCNLIFMVLDCLKPLGHKKLLEHELEGFGIRLNKKPPNIY
YKRKDKGGINLNSMVPQSELDTDLVKTILSEYKIHNADITLRYDATSDDLIDVIEGNRIY
IPCIYLLNKIDQISIEELDVIYKIPHCVPISAHHHWNFDDLLELMWEYLRLQRIYTKPKG
QLPDYNSPVVLHNERTSIEDFCNKLHRSIAKEFKYALVWGSSVKHQPQKVGIEHVLNDED
VVQIVKKV
regards,
stephan
    
    
More information about the Bioperl-l
mailing list