Question: Using efetch(biopython) to get the fasta file of an mRNA, as well as the exon locations
0
gravatar for Jacob
2.7 years ago by
Jacob10
Jacob10 wrote:

I want to use some type of e-utility tool to obtain the fasta mRNA sequence file of a gene. First I'm getting the uid from Gene with the code below

from Bio import Entrez
Entrez.email="email@address.com
handle = Entrez.esearch(db="gene", term="Acan[Gene Name] AND Homo sapiens[Organism]", rettype='fasta')"
record = Entrez.read(handle)
id = record["IdList"]

Then I would like to use this uid to get the fasta sequence for the mRNA gene. I would also like to obtain the bp positions of the exons in this gene and the bp positions of the longest off. This is how I'm trying to use efetch

Entrez.efetch(db="nucleotide", term=id[0]+"[uid]", rettype="fasta")

It is not working and I don't know how to specify to get the exon positions either. I know the exon positions are there because I can see them on the website.

https://www.ncbi.nlm.nih.gov/nuccore/NM_000046.3
entrez biopython efetch python • 1.4k views
ADD COMMENTlink modified 2.7 years ago by Ben50 • written 2.7 years ago by Jacob10
0
gravatar for Ben
2.7 years ago by
Ben50
Ben50 wrote:

Firstly, you should have the GTF/GFF and fasta files of the Reference genome; Then you can extract the cDNA sequence of mRNA with custome script; Last, you should convert the cDNA sequence to mRNA sequence.

ADD COMMENTlink written 2.7 years ago by Ben50

Is there a way to get all the exon regions from that website using a command line tool?

ADD REPLYlink written 2.7 years ago by Jacob10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1561 users visited in the last hour