This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using Web Services To Get Fasta Files For Refseq Nm And Nr Identifiers

I would like to be able to fetch fasta files for given refseq cdna or rna sequences programatically

I can do this for EMBL cdnas using a wget command like

wget  -q -O - "http://www.ebi.ac.uk/ena/data/view/BX648399.1&display=fasta" > BX648399.1.fasta

Does anyone know if there is an equivalent service at the NCBI

thanks

refseq ncbi fasta
• 4,090 views • link
updated by Biostar Biostar • written by Laura Laura

2 answers

Use NCBI Efetch:

$ curl -s "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nucleotide&id=BX648399&rettype=fasta&retmode=txt" | head
>gi|34367561|emb|BX648399.1| Homo sapiens mRNA; cDNA DKFZp686M11224 (from clone DKFZp686M11224)
GGCCACGGCCGACATGTTTTTTTTTTTTTTTTTTCTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT
TTTTTTTTTTTTTTTTTTTTTTTTAACACAATGGCCCTGCCTCCCACCGCTTTATTTCTTTCGGTTTCGG
ATGCAAAACAAAAAATTTTAAAAGAAAATGTGACTTCAAAGGAAAAGAACAAATTTCCAAAGACTTGGGG
GAGTGAAGGCAGAGCCTGGTGCAGATGGACGAGGTCTGCAGACGGAGGGCAGAGGTGGTGGAAGGGGCCA
GGGGCCTGCAGGCCTCCCCCTGGAACTGGGACTGGTCTCGGTCTGCTGACGTCAGGGTCAGCTCCCCCGC
GGAGCTGACTTCAGCAGCCCACAGCTGTGGGGCTTCAGCAGCCACACCAGCCCAGCCCAGCCCAGCTCTC
GATACGTTTGGTCTTTCATGCTGAAAAATAA
• 116 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

fantastic thanks

• 0 views • link
written by Laura Laura

Yes you can use eFetch: http://www.ncbi.nlm.nih.gov/books/NBK25498/#chapter3.ESearch__ESummaryEFetch

• 162 views • link
written by Damian Kao Damian Kao

Log in to answer this question.

More posts like this
  • Synteny Diagram In Biopython
    written by Pappu Pappu
    <p>I am trying to compare the synteny of a homologous gene of a set of embl files. I wrote the following script in Biopython which …
  • Mpileup Command Line
    written by Giffredo • Giffredo
    <p>Hello,</p> <p>I have a problem to use mpileup. I used this command:</p> <pre><code>samtools mpileup -f mi.fas sorted.bam &gt; X.bcf &amp; bcftools view X.bcf &gt; X.vcf …
  • Customize Link To Ncbi Blast Search
    written by Alex Alex
    <p>I want to insert a link to NCBI Blast search into web page with given GI and coordinates. Link like this:</p> <pre><code>&lt;a href="http://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Nucleotides&amp;PROGRAM=blastn&amp;QUERY={{ GI }}&amp;DATABASE=nr&amp;MEGABLAST=on&amp;BLAST_PROGRAMS=megaBlast&amp;LINK_LOC=nuccore&amp;PAGE_TYPE=BlastSearch&amp;QUERY_FROM={{ …
  • Printing Non Coding Sequence Preceeding A Cds Of Embl File
    written by Pappu Pappu
    <p>I have a huge embl file containing many sequences. I want to print non coding sequence before a CDS. For example I want to print …
  • What Better Way To Get Consensus Sequence Of Mapping Against Reference Genome Using Samtools Mpileu…
    written by thiagomafra thiagomafra
    <p>Hi people,</p> <p>I have a set of reads mapped against a reference genome and need take consensus sequence using <a href='http://samtools.sourceforge.net/'>samtools</a> mpileup. I used the …
  • Parsing Non-Coding Region Around Protein Of Interest From Embl File
    written by Pappu Pappu
    <p>I want to parse noncoding DNA sequences around a protein (P08707) from an embl file: <a href="http://www.ebi.ac.uk/ena/data/view/U32222&amp;display=txt&amp;expanded=true">http://www.ebi.ac.uk/ena/data/view/U32222&amp;display=txt&amp;expanded=true</a> </p> <p>I could grep '^ CDS' and then …
  • Nominal Length And Tophat Inner Mate Distance
    written by narges • narges
    <p>Hi, I wanted to know whether I can use the nominal length provided in EBI website as the tophat inner mate distance or not?like the …
  • Explanation Of Low Mapping Quality!
    written by khikho khikho
    <p>This is my pipeline for the human reads to map with human reference: (I use default solid2fastq.pl script which is in <a href='http://bio-bwa.sourceforge.net/'>BWA</a> package)</p> <pre><code>#! …
  • Translating A Mitochondrial Cdna
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>Here is the cDNA for the UCSC knowngene <strong>uc011mfh.1</strong> , a protein said to be on the <strong>chrM</strong>:5,855-7,427. <a href="http://genome.ucsc.edu/cgi-bin/hgTracks?org=Human&amp;db=hg19&amp;position=uc011mfh.1">http://genome.ucsc.edu/cgi-bin/hgTracks?org=Human&amp;db=hg19&amp;position=uc011mfh.1</a></p> <pre><code>ATGATCTGCTGCAGTGCTCTGAGCCCTAGGATTCATCTTTCTTTTCACCG TAGGTGGCCTGACTGGCATTGTATTAGCAAACTCATCACTAGACATCGTA CTACACGACACGTACTACGTTGTAGCTCACTTCCACTATGTCCTATCAAT AGGAGCTGTATTTGCCATCATAGGAGGCTTCATTCACTGATTTCCCCTAT TCTCAGGCTACACCCTAGACCAAACCTACGCCAAAATCCATTTCACTATC …
  • Fasta Format Sequence Parsing With Bioperl
    written by Raghul Raghul
    <p>Hi I am having a fasta format sequence but the header line is not like NCBI description(please see below). </p> <pre><code>&gt;contig07303 length=480 numreads=24 gene=isogroup00008 status=isotig …