[Biopython] Problem using efetch

Peter biopython at maubp.freeserve.co.uk
Tue Jul 14 18:08:50 EDT 2009


On Tue, Jul 14, 2009 at 4:39 PM, bar tomas<bartomas at gmail.com> wrote:
> Hi,
>
> I’m using BioPython to access Entrez databases.  I’m following
> the BioPython tutorial. I’ve tried retrieving all record ids from
> pcassay database with esearch and then retrieving the first full
> record on the list with efetch:
>
> handle = Entrez.esearch(db="pcassay", term="ALL[filt]")
>
> print record["IdList"]
>
> # This prints the following list of ids:
>
> # ['1866', '1865', '1864', '1863', '1862', '1861', '1033', '1860', etc.
>
>
> But when I then try to retrieve the first record:
>
> handle2 = Entrez.efetch(db="pcassay", id="1866")
>
> I get the following error :
>
> <html>
> <body>
> <br/><h2>Error occurred: Report 'ASN1' not found in 'pcassay'
> presentation</h2><br/><ul title="some params from request:">
> <li>db=pcassay</li>
> ...
>
> Do you have an idea of what I’m doing wrong?

This isn't anything wrong with Biopython - this is the sort of
slightly cryptic error the NCBI gives when the return type
and/or return mode isn't supported. Apparently the default
(ASN1) isn't supported for this database. The NCBI efetch
documentation is a little vague or simply missing for the
less main-stream databases. You can make some
guesses from playing with the Entrez website, e.g.

>>> print Entrez.efetch(db="pcassay", id="1866", rettype="uilist").read()
<html><head><title>PmFetch response</title></head><body>
<pre>
1866
</pre></body></html>
>>> print Entrez.efetch(db="pcassay", id="1866", rettype="uilist", retmode="text").read()
1866
>>> print Entrez.efetch(db="pcassay", id="1866", rettype="abstract", retmode="text").read()

1: AID: 1866
Name:  Epi-absorbance-based counterscreen assay for selective VIM-2 inhibitors:
biochemical high throughput screening assay to identify inhibitors of TEM-1
serine-beta-lactamase.
Source:  The Scripps Research Institute Molecular Screening Center
Description:   Source (MLPCN Center Name): The Scripps Research Institute
...

You could also try emailing the NCBI for advice.

Peter



More information about the Biopython mailing list