[Bioperl-l] $hit->accession

suhoiy suhoiy@21cn.com
Fri, 1 Nov 2002 11:17:06 GMT


Hello all,

In Bio::SearchIO::blast module, the accession number is obtained by:
        my @pieces = split(/\|/,$id);
        my $acc = pop @pieces;

but some hits are different, for example:
        gi|3024862|sp|P76372|WZZB_ECOLI Chain length determinant protein...   551   e-156
        gi|18266411|gb|AAL67565.1|AF461121_16 (AF461121) O-antigen chain...   538   e-152
        gi|584950|sp|P37792|WZZB_SHIFL Chain length determinant protein ...   535   e-151

so the returned $hit->accession are "WZZB_ECOLI" and so on, while the true 
accessions are "P76372" and so on.

Is it a bug? maybe we can get the accessions by /gi\|\d+\|\w+\|([\w\.]*)/ ?

thanks.

suhoiy