[Bioperl-l] Beginners HOWTO query a range of lengths 0:3000[SLEN]
    Alexander Ptok 
    alexanderptok at web.de
       
    Wed Jan  9 15:34:56 UTC 2008
    
    
  
Hi,
I am a beginner to BioPerl and working through the Beginners HOWTO
Version of BioPerl is 1.4-1 running on Debian etch
In the Howto everything worked fine until the section
Retrieving multiple sequences from a database
from where i copied the following script:
use Bio::DB::GenBank;
use Bio::DB::Query::GenBank;
 
$query = "Arabidopsis[ORGN] AND topoisomerase[TITL] and 0:3000[SLEN]";
$query_obj = Bio::DB::Query::GenBank->new(-db    => 'nucleotide',  -query => $query );
 
$gb_obj = Bio::DB::GenBank->new;
 
$stream_obj = $gb_obj->get_Stream_by_query($query_obj);
 
while ($seq_obj = $stream_obj->next_seq) {    
    # do something with the sequence object    
    print $seq_obj->display_id, "\t", $seq_obj->length, "\n";
}
If i cut the 0:3000[SLEN] query it works and returns a lot of sequences, when i alter the query to e.g. 1830[SLEN] it
finds the one sequence that has the length 1830, but i was not able to query a range of lengths.
Please, does anyone know what i am doing wrong.
Greetings
A. Ptok
_________________________________________________________________________
In 5 Schritten zur eigenen Homepage. Jetzt Domain sichern und gestalten! 
Nur 3,99 EUR/Monat! http://www.maildomain.web.de/?mc=021114
    
    
More information about the Bioperl-l
mailing list