Collections:
Fetch Sequences from NCBI with Bio.Entrez.efetch()
How to Fetch Sequences from NCBI with the Bio.Entrez.efetch() function?
✍: FYIcenter.com
Bio.Entrez.efetch() function allows you to fetch DNA or protein sequences from NCBI databases: PubMed, GenBank, GEO, and many others. It uses the Entrez Web services provided by www.ncbi.nlm.nih.gov.
Here is an example on how to Fetch Sequences in FASTA format from NCBI with the Bio.Entrez.efetch() function.
fyicenter$ python >>> from Bio import Entrez >>> from Bio import SeqIO >>> Entrez.email = "A.N.Other@example.com" >>> >>> with Entrez.efetch( ... db="nucleotide", rettype="fasta", retmode="text", id="6273291" ... ) as handle: ... seq_record = SeqIO.read(handle, "fasta") ... >>> print(seq_record) ID: AF191665.1 Name: AF191665.1 Description: AF191665.1 Opuntia marenae rpl16 gene; chloroplast gene for chloroplast product, partial intron sequence Number of features: 0 Seq('TATACATTAAAGGAGGGGGATGCGGATAAATGGAAAGGCGAAAGAAAGAAAAAA...AGA')
You can also fetch the same sequence in GenBank format, which contains more properties.
>>> Entrez.email = "A.N.Other@example.com" >>> with Entrez.efetch( ... db="nucleotide", rettype="gb", retmode="text", id="6273291" ... ) as handle: ... seq_record = SeqIO.read(handle, "gb") ... >>> print(seq_record) ID: AF191665.1 Name: AF191665 Description: Opuntia marenae rpl16 gene; chloroplast gene for chloroplast product, partial intron sequence Number of features: 3 /molecule_type=DNA /topology=linear /data_file_division=PLN /date=07-NOV-1999 /accessions=['AF191665'] /sequence_version=1 /keywords=[''] /source=chloroplast Grusonia marenae /organism=Grusonia marenae /taxonomy=['Eukaryota', 'Viridiplantae', 'Streptophyta', 'Embryophyta', ...] /references=[Reference(title='Phylogeny of the subfamily Opuntioideae (Cactaceae)', ...), ...] Seq('TATACATTAAAGGAGGGGGATGCGGATAAATGGAAAGGCGAAAGAAAGAAAAAA...AGA') >>> print(seq_record.features[0]) type: source location: [0:902](+) qualifiers: Key: db_xref, Value: ['taxon:106980'] Key: mol_type, Value: ['genomic DNA'] Key: note, Value: ['subfamily Opuntioideae; synonym: Marenopuntia marenae, Grusonia marenae'] Key: organelle, Value: ['plastid:chloroplast'] Key: organism, Value: ['Grusonia marenae']
⇒ Search for Related Items in NCBI with Bio.Entrez.elink()
⇐ Retrieve Record Summary with Bio.Entrez.esummary()
2023-08-25, 287🔥, 0💬
Popular Posts:
Why am I getting no results with "babel -i fs ... -s query_file" command? It should perform a substr...
Molecule Summary: ID: FYI-1000994 SMILES: Received at FYIcenter.com on: 2021-09-13
Molecule Summary: ID: FYI-1002013 Names: InChIKey: RZTAMFZIAATZDJ-UHFFFAOYS A-NSMILES: CCOC(=O)/C1=C...
Molecule Summary: ID: FYI-1000306 SMILES: COCCNC(=O)C(=O)N/N=C/c1c ccc(c1OC)OCReceived at FYIcenter....
Molecule Summary: ID: FYI-1001986 SMILES: NC(C1=C[N+]([C@@H]2O[C@@ H]([C@H]([C@H]2O)O)COP([ O-])(OP(O)...