[Bioperl-l] Error parsing genbank file

Koul, Sunita sunita at pathology.wustl.edu
Tue Aug 13 17:12:47 UTC 2013


All,
I am trying to retrieve some gene information from genbank file.
Part of the code is to ignore the pseudo genes


if( $featobj->primary_tag eq 'gene' )
                                             {
                                               print "$geneID\n";


                                                            if(scalar($featobj->get_tag_values('pseudo')) > 0)
                                    {
                                                                           $pseudo{$geneID} = 1;

                                                                           next;
                                                            }
}

This throws me an error

What is the reason for this error



MSG: asking for tag value that does not exist pseudo
STACK: Error::throw
STACK: Bio::Root::Root::throw /usr/local/share/perl5/Bio/Root/Root.pm:472
STACK: Bio::SeqFeature::Generic::get_tag_values /usr/local/share/perl5/Bio/SeqFeature/Generic.pm:522



If I comment out the above code -I get another error for this part of the error with the message that
MSG: asking for tag value that does not exist 'db_xref'


                              foreach my $featobj ($seqobj->all_SeqFeatures() )
                              {
                                             my $geneID;

                                             for my $value ($featobj->get_tag_values('db_xref'))
      {
               if ($value =~ /GeneID:(.*)/)
               {
                              $geneID = $1;
               }
      }

Please help me understand what is wrong in this code


The materials in this email are private and may contain Protected Health Information. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this email in error, please immediately notify the sender via telephone or return email.




More information about the Bioperl-l mailing list