[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