I know similar questions have been posted before..
How To Retrieve All Sequences, From Ncbi, That Belong To A Specific Txid And Its Sub Txids?
C: Refseq Proteins For A Given Taxid
But, I am having trouble retrieving the children sequences of a given taxID.
For instance,
from Bio import Entrez
record = Entrez.read(Entrez.esearch(db='protein', term="txid1392[Organism]"))
record['IdList']
Returns just one list of protein UIDs for the Bacillus anthracis species at 1392, not the list for each organism that is below this taxon. Thus, Entrez.efetch only returns one set of protein sequences.
Dropping the [Organism] doesn't change this behavior. Am I missing something?
•
written
by
yarmda •
1 answer
Your code is fetching only the UIDs from the first page. You need to provide retmax= parameter to fetch all records. See below corrected code,
from Bio import Entrez
record = Entrez.read(Entrez.esearch(db='protein', retmax=770094, term="txid1392[Organism]"))
record['IdList']
I have put retmax=770094 as this taxon has 770094 protein records.
Alternative,
You can get the list of sequences for txid1392[Organism] from web NCBI also.
- Go to NCBI Entrez search and type
txid1392[Organism]and chooseproteindatabase from dropdown list (it will fetch all protein sequences fortxid1392[Organism] - Go to
send tobutton and send to file (You can choose FASTA format for downloading the sequences)
Log in to answer this question.
-
How to query NCBI to extract Virus fasta files using BioPython?
written by Bioinformatics_begginner •
-
Retrieving taxonomy from entrez search in biopython
written by pramirez •
-
Collecting CDS from transcripts using Entrez
written by Taylor •
-
Jupyter python cell will not finish to completion how to fix hash(id(seq) issue
written by bcheda •
-
download sequences form ncbi using python
written by HZZ0036 •
-
Using Biopython to search through a list of authors on pubmed
written by flisstee •
-
Downloading Gene Annotation protein.faa files from Genome Accession numbers from NCBI via Entrez
written by ijc2 •
-
Retrieve Genbanks From Taxid
written by JJK •
-
Biopython Entrez.Efetch Problem
written by a11msp
-
Refseq Proteins For A Given Taxid
written by Chris •
Unfortunately, I think you might have to list all the child taxon identifiers explicitly - but try exploring the web interface for building an advanced query first in case that shows a better solution.
Actually, I think it may have been an issue with a default retmax of 20.
Oh good. I should have tried the example myself really to confirm my hunch. Thanks!
hy ,i am currenly doing biopython yet in industry does it have influence and impotance in this era?
Hi, this comment is not appropriate to this (very old) thread.
If you wish to ask a question, please create your own thread. If you do, I strongly encourage you to search the forum first (since questions like this are asked often - and are of dubious usefulness). If you cannot find something that satisfies you, ask a question, but please add much more information and detail and make the question as specific as possible.