This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Can't get all coding sequences from list of protein IDs [Entrez Direct]

I have a list of protein ids in number format (e.g. 25121878) that I want to retrieve the coding sequences for.

This E-direct command is working for some, but for some it's giving me 'NO RESULT':

efetch -db protein -format fasta_cds_na -id **id**

Am I doing something wrong or is there really no way to get a CDS in these instances?

Reproducible example where I get no result:

efetch -db protein -format fasta_cds_na -id 25121878
entrez entrez-direct
• 986 views • link
updated by GenoMax GenoMax • written by charlieclark1ee • charlieclark1ee

1 answer

Problem with 25121878 (I assume this is a gi ID which is deprecated for end-user use now) is that this is a conceptual translation so there is no nucleotide record associated with it.

COMMENT     PROVISIONAL REFSEQ: This record has not yet been subject to final
            NCBI review. The reference sequence was derived from
            NP_040350:647-922.
            Related sequence: M20562.
            Method: conceptual translation.
• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this