[Bioperl-l] Get nucleotide sequence when expecting protein from genpept

[Qunfeng]

This particular protein record (acc#T16005) was imported from PIR. In other 
words, this is not an original GenBank protein record. When GenBank imports 
protein records from other DB, it keeps their original acc#.

However, gi# should be unique.

Q

At 04:23 PM 7/11/2006, Frederick Partridge wrote:


>I am trying to retrieve various protein sequences from genpept using
>get_Seq_by_acc. All of them work ok, except one T16005:
>
>
>If I try and retrieve it with a reduced program:
>
>
>#!usr/bin/perl -w
>
>use strict;
>
>use Bio::Perl;
>use Bio::SeqIO;
>
>my $genpept = new Bio::DB::GenPept;
>
>my $seq = $genpept->get_Seq_by_acc('T16005');
>
>print ($seq->seq(),'\n');
>
>
>
>I get back a nucleotide sequence, which is another sequence at NCBI with
>the same accession number. (I thought these were meant to be unique? but
>evidently not.)
>
>
>I am using bioperl 1.5.1, perl 5.8.1, Mac OS 10.3
>
>
>Could anyone help me to get this protein sequence with my program?
>
>
>Many thanks,
>
>
>
>Freddie Partridge
>
>University of Oxford
>
>
>_______________________________________________
>Bioperl-l mailing list
>Bioperl-l at lists.open-bio.org
>http://lists.open-bio.org/mailman/listinfo/bioperl-l