I need to retrieve a set of protein and mRNA sequences
1 answer
Using EntrezDirect:
$ esearch -db gene -query '7157' | elink -db gene -target nuccore -name gene_nuccore_refseqrna | efetch -format fasta | grep ">"
>NR_176326.1 Homo sapiens tumor protein p53 (TP53), transcript variant 14, non-coding RNA
>NM_001407264.1 Homo sapiens tumor protein p53 (TP53), transcript variant 10, mRNA
>NM_001407263.1 Homo sapiens tumor protein p53 (TP53), transcript variant 9, mRNA
>NM_001407268.1 Homo sapiens tumor protein p53 (TP53), transcript variant 12, mRNA
>NM_001407269.1 Homo sapiens tumor protein p53 (TP53), transcript variant 12, mRNA
>NM_001407271.1 Homo sapiens tumor protein p53 (TP53), transcript variant 13, mRNA
>NM_001407270.1 Homo sapiens tumor protein p53 (TP53), transcript variant 13, mRNA
For proteins
$ esearch -db gene -query '7157' | elink -db gene -target protein -name gene_protein_refseq | efetch -format fasta | grep ">"
>NP_001394193.1 cellular tumor antigen p53 isoform a [Homo sapiens]
>NP_001394192.1 cellular tumor antigen p53 isoform g [Homo sapiens]
>NP_001394197.1 cellular tumor antigen p53 isoform b [Homo sapiens]
>NP_001394198.1 cellular tumor antigen p53 isoform i [Homo sapiens]
>NP_001394200.1 cellular tumor antigen p53 isoform i [Homo sapiens]
>NP_001394199.1 cellular tumor antigen p53 isoform b [Homo sapiens]
Option 2:
Use datasets: https://www.ncbi.nlm.nih.gov/datasets/gene/
Click on the search by identifiers tab and paste your identifiers in (may need to do in batches).
On results page, select gene, click Downloads --> Download package --> Choose gene/protein/transcript sequence --> Download data package.
Can also be done using command line datasets.
Log in to answer this question.
More posts like this
-
How to do KEGG pathway analysis when I have a gene with multiple entrez IDs?
written by bioinfo •
-
Retrieve subsequences from a list of gene sequences
written by keki •
-
Retrieve nucleotide ID from protein ID
written by kmears •
-
Retrieving gene ID using transcripts ID from Entrez database using CLI or Batch Entrez
written by hans •
-
Extract Consensus Sequence from DIAMOND
written by c.pouchon •
-
DAVID gene enrichment - multiple Entrez IDs per gene symbol, what to do?
written by BioBing
-
Using Batch Entrez to extract .gb files with specific coordinates
written by vykintas.jauniskis •
-
Get fastas for genes by gene symbol
written by acorella •
-
Batch Retrieval of Coding Sequences
written by fire_water
-
Promoter Sequences For Non Model Organism From Old Gene Build
written by Sudeep
What have you tried?
I tried using curl because I saw a similar post here but that didn't go well so now I am trying the answer of GenoMax bellow
curl is a command line utility that can do a range of operations. What exactly did you do, what did you expect and what errors did you run into? Always include these details when you ask for help.