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, 857🔥, 0💬
Popular Posts:
Can I load JSME JavaScript code in HTML "body" element instead of "head"? Yes, you can load JSME Jav...
Molecule Summary: ID: FYI-1005620 Names: InChIKey: NIKHGUQULKYIGE-XSBFNCLBS A-NSMILES: C=C3C[C@@]24C...
Molecule Summary: ID: FYI-1002248 Names: InChIKey: VUBWTGOADLGEEC-UHFFFAOYS A-NSMILES: Nc2ccc(c1[nH]...
Molecule Summary: ID: FYI-1000196 SMILES: O=C(N(C(C(OC(=O)CCCCC(=O )OC)([H])[H])([H])[H])[H ])C2=C([N]...
Molecule Summary: ID: FYI-1001810 SMILES: [2H]C1=C(C(=C(C(=C1[2H]) [2H])C(=C([2H])[2H])[2H] )[2H])[2H]...