How to using GenBank assembly accession to parse out the FASTA file
I have 100 randomly selected genomes from a dataset. Then, I want to use the GenBank assembly accession GCA_001600695.1 to parse out the fasta file for his gene. Below is my code, but did not work. Anyone have a suggestion about what command in biopython I can use?
from Bio import Entrez
handle = Entrez.efetch(db="nucleotide", id="GCA_000021505.1", rettype="fasta", retmode="text")
print(handle.read())
• 1,897 views
•
link
updated
by
Ram
•
written
by
ricordo.yan •
•
written
by
ricordo.yan •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Getting cds files from ncbi using accession numbers
written by lily •
-
Entrez.efetch with retmax=1000000 not retrieving more than 10000 records (python)
written by Aisha •
-
Seqrecord object to only a single line sequence (without ID)
written by voltosj •
-
Having issues with having extra information on my Fasta file.
written by trejomarco6 •
-
How do I make my code go through all through my list in my file.
written by trejomarco6 •
-
parse genome "start point to end point" to find what genes inside? blast in python
written by seok1213neo •
-
How to get all transcript and CDS from gene_id in biopython?
written by ishenweiyan •
-
how to get protein sequence biopython
written by Xylanaser
-
download sequences form ncbi using python
written by HZZ0036 •
-
How To Retrieve Multiple Sequence In One Python Script
written by ahmedakhokhar
Do you see an error message? What does it say? What tells you that the code does not work? Please see: Ten Simple Rules for Getting Help from Online Scientific Communities
That was probably a leftover back-tick from formatting here on the website. Can you try again please?
GCA_001600695.1 is actually assembly accession for Ascoidea asiatica. Actual genome sequence is contained in this link. You may want to try some of the ideas in this thread: Cannot get efetch to download genome - what is wrong?