[Bioperl-l] getting gi number

[Jason Stajich]

#!/usr/bin/perl -w
use strict;
use Bio::DB::GenBank;
use Bio::DB::GenPept;
# if sequence is protein accession
# use Bio::DB::GenBank if it is nt accession
my $db = Bio::DB::GenPept->new();

my $seq = $db->get_Seq_by_acc('AAA83776');

print $seq->display_id, " ", $seq->primary_id, "\n";

-jason

On Wed, 9 Jun 2004, Sucheta Tripathy wrote:

> Hello all,
>
> I have the accession number in this format gb|AAA83776.1| and would like to
> get the gi number for that. Is there a method in bioperl that takes care of
> that?
>
> Thanks
>
> Sucheta
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>

--
Jason Stajich
Duke University
jason at cgt.mc.duke.edu