[Bioperl-l] Bio::DB::GenBank problem?
James Wasmuth
james.wasmuth at ed.ac.uk
Tue Sep 9 05:04:05 EDT 2003
Hi Wayne,
I identified this a few months back, and from what I can remember Heikki
fixed it and is now in the CVS.
james
wxu wrote:
>Hello,
>I wrote a simple bioperl program to retrieve genBank sequences by genBank
>accession numbers using $gb->get_Stream_by_id function.
>But it seems to me it can only work on GI number, not work on accession
>numbers. I use $gb->get_Stream_by_acc, the result is the same. In my
>following code, it only give me 2981014 sequences. I search NCBI web
>manually, all "J00522","AF303112","2981014","BG065024","BG065143" work.
>And I tried $gb->get_Seq_by_acc("J00522"); It complains "No such acc error".
>
>What is the problem?
>Your help would be greatly appreciated.
>
>Wayne
>--
>
>----------------------------------------------------------
>#!/usr/local/bin/perl
>use lib "/usr/local/bioperl/bioperl-1.2";
>use lib qw(/home/wxu/bioperl);
>
>use Bio::Perl;
>use Bio::DB::GenBank;
>use Bio::SeqIO;
>my $gb = new Bio::DB::GenBank();
>
>my $seqio =
>$gb->get_Stream_by_id(["J00522","AF303112","2981014","BG065024","BG065143"])
>;
>
>while( my $seq = $seqio->next_seq ) {
> #print "seq length is ", $seq->length,"\n";
> write_sequence('>>outputfile', 'fasta', $seq)
> }
>----------------------------------------------------------------------------
>--
>Wayne Xu
>Computational Genomics Specialist www.msi.umn.edu/user_support/compgen
>
>Supercomputing Institute
>550 Walter Library
>117 Pleasant Street SE
>University of Minnesota
>Minneapolis, Minnesota 55455
>email: wxu at msi.umn.edu help email: help at msi.umn.edu
>phone: 612-624-1447 help phone: 612-626-0802
>fax: 612-624-8861
>
>
>_______________________________________________
>Bioperl-l mailing list
>Bioperl-l at portal.open-bio.org
>http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
>
--
Nematode Bioinformatics
Blaxter Nematode Genomics Group
Institute of Cell, Animal and Population Biology
Ashworth Labs
University of Edinburgh
King's Buildings
Edinburgh
EH9 3JT
UK
(+44)(0)131 650 7403
More information about the Bioperl-l
mailing list