[Bioperl-l] Bio::DB::Query::GenBank retrieves fewer sequences
	than Webbrowser query
    Paulo Almeida 
    paulo.david at netvisao.pt
       
    Tue Mar 23 17:06:19 EST 2004
    
    
  
I think you need to escape the brackets in Organism. Try replacing your 
query line with this:
(-query   =>'Mus\[Organism\] AND exon NOT mRNA NOT cDNA',   
Hope it works,
-Paulo Almeida
Jürgen Schmitz wrote:
The following query  retrieves less than a hundred sequences, much less 
than I get via Entrez:
my $query = Bio::DB::Query::GenBank->new
        (-query   =>'Mus[Organism] AND exon NOT mRNA NOT cDNA',
         -db      => 'nucleotide');
$seqio = $gb->get_Stream_by_query($query);
while( my $seq =  $seqio->next_seq ) {...}
but entering "Mus[Organism] AND exon NOT mRNA NOT cDNA" at
http://www.ncbi.nlm.nih.gov:80/entrez/query.fcgi?CMD=search&DB=nucleotide
retrieves more than 5000 sequences.
What could be the problem ?
There are no error messages printed to STDERR.
I'm using Bioperl 1.4, Perl v5.8.1-RC3 on a MAC (OSX)
    
    
More information about the Bioperl-l
mailing list