This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Linking NCBI gene ID's to nucleotide sequence

I've been trying to get the nucleotide sequence for a list of genes on the command line, for example for the gene LOC100136426 I would like to get the following:

https://www.ncbi.nlm.nih.gov/nuccore/NC_027322.1?report=fasta&from=3409007&to=3419329&strand=true

I have been trying to use eutils to do this with the following command:

$ esearch -db nucleotide -query "LOC100136426" | efetch -format fasta

My issue is that this gives three results, where as I only want the nucleotide sequence for the gene so that I can perform a blast.

I'm not quite sure how to implement a filter in the pipe to only give the gene sequence. I'd greatly appreciate any help or even suggestions to alternate methods to find this information.

rna-seq ncbi eutils
• 2,723 views • link
updated by Jean-Karim Heriche Jean-Karim Heriche • written by jared.j.tromp • jared.j.tromp

1 answer

You need to query the gene database not the nucleotide one or find out what is the corresponding identifier in the nucleotide database.

• 0 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

I couldn't get the pipe to return the nucleotide sequence when using the gene database, although this would correctly identify the gene. I've corrected the command using the GI number as you suggested and this returns the correct result, thanks for the help!

Is there a way to use the gene database to return the fasta sequence? I have tried:

$ esearch -db gene query "LOC100136426" | efetch -format fasta

This didn't return the fasta sequence

• 0 views • link
written by jared.j.tromp • jared.j.tromp

Log in to answer this question.

More posts like this