[Bioperl-l] Parsing kegg files
    stephen.marshall at novartis.com 
    stephen.marshall at novartis.com
       
    Wed Feb 21 17:01:00 UTC 2007
    
    
  
Hello
I"m trying to parse a Kegg file and I can't seem to get at the pathway 
information... Here's a snippet of my code. I only see dblink and 
description as annotation
use Bio::SeqIO;
my $stream = Bio::SeqIO->new(-file => $filename, -format => 'KEGG');
while ( my $seq = $stream->next_seq() ) {
        # do something with $seq
        my $id = $seq->display_id();
        print "$id:";
        my $ann = $seq->annotation();
        foreach my $key ( $ann->get_all_annotation_keys() ) {
                my @values = $ann->get_Annotations($key);
                foreach my $value ( @values ) {
                        print "Annotation: ",$key," value: 
",$value->as_text,"\n";
                }
        }
}
_________________________
CONFIDENTIALITY NOTICE
The information contained in this e-mail message is intended only for the 
exclusive use of the individual or entity named above and may contain 
information that is privileged, confidential or exempt from disclosure 
under applicable law. If the reader of this message is not the intended 
recipient, or the employee or agent responsible for delivery of the 
message to the intended recipient, you are hereby notified that any 
dissemination, distribution or copying of this communication is strictly 
prohibited. If you have received this communication in error, please 
notify the sender immediately by e-mail and delete the material from any 
computer.  Thank you.
    
    
More information about the Bioperl-l
mailing list