This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ensembl REST to download transcript fastas

Sorry in advance if this is a newbie question, but I always get a little tripped up with the Ensembl Rest API.

I want to download fastqs given a stable id for a gene transcript, like so:

https://uswest.ensembl.org/Mus_musculus/Export/Output/Transcript?db=core;flank3_display=0;flank5_display=0;g=ENSMUSG00000029648;output=fasta;r=5:147498414-147662821;strand=feature;t=ENSMUST00000031652;param=cdna;param=coding;param=peptide;param=utr5;param=utr3;param=exon;param=intron;genomic=unmasked;_format=HTML.

Where it returns separate fasta entries for exons and introns and cdna. Is this possible with a ensembl query?

ensembl
• 1,143 views • link
updated by Emily Emily • written by H Miller • H Miller

1 answer

Not ideally. You can use the sequence/id endpoint with the masking option to get the cDNA with UTRs in lower case, eg: https://rest.ensembl.org/sequence/id/ENSMUST00000031652?content-type=text/x-fasta;multiple_sequences=1;type=cdna;mask_feature=1

You can also get the genomic sequence of the transcript with the introns in lower case: https://rest.ensembl.org/sequence/id/ENSMUST00000031652?content-type=text/x-fasta;multiple_sequences=1;type=genomic;mask_feature=1

• 0 views • link
written by Emily Emily

Log in to answer this question.

More posts like this
  • Ensembl IDs
    written by jensen.416 • jensen.416
    I'll start with my end goal: I want to add RNAseq expression data to a whole exome sequencing dataset (VCF file). Both of which are …
  • Searching a DNA sequence for eukaryotic genes
    written by mohammad.shoumik.saad.aarman • mohammad.shoumik.saad.aarman
    How to find eukaryotic genes from [this][1] FASTA format of a sequence in human chromosome 13? [1]: https://asia.ensembl.org/Homo_sapiens/Export/Output/Location?db=core;flank3_display=0;flank5_display=0;output=fasta;r=13:136027-146027;strand=1;intron=yes;genomic=unmasked;cdna=yes;peptide=yes;exon=yes;utr3=yes;utr5=yes;coding=yes;_format=HTML
  • Question about getSequence function of biomaRt in R
    written by Ruby • Ruby
    Hi everyone, I am working on a project about copy number variant calling algorithms. One question I am asking is that if GC content correlated …
  • No Prominent Change in Structure of Mutated Proteins
    written by ammarsabir15 ammarsabir15
    I am working on Parkinson disease causing protein **PINK1**, I am focusing on **26 disease causing variants** for parkinson. These variants are reported [here][1]. There …
  • 200Bp Long Query Returning 201Bp
    written by secretjess secretjess
    <p>In Ensembl when I use this <a href="http://www.ensembl.org/Homo_sapiens/Export/Output/Location?db=core;flank3_display=0;flank5_display=0;output=fasta;r=X:200000-200200;strand=1;coding=yes;cdna=yes;utr3=yes;peptide=yes;intron=yes;exon=yes;genomic=unmasked;utr5=yes;_format=Text">link</a> to return the sequence on the sense strand of X chromosome between 200000 and 200200 I get …
  • Difference Between Biomart Query And The Ensembl Database
    written by Nicolas Rosewick Nicolas Rosewick
    <p>Hi,</p> <p>I'm using the biomaRt R packages (bioconductor) to retrieve the 3'UTR sequences of a list of genes. I've the entrezgene Id for each of …