This is a test version of Biostars. For the public version, visit https://www.biostars.org.
From NCBI Taxonomy to 16s RNA (automated)

I have a list of NCBI Taxonomy IDs and I want to retrieve the corresponding 16S RNA sequence. I can do this for every single element using SILVA https://www.arb-silva.de/. But I need an automated way. As far as I know, SILVA does not have an REST API, but maybe other databases are accessible in a programmatic way?

16s rna taxonomy

Can you post a few examples of the NCBI IDs that you have?

Are you able to use these IDs to locate the 16s RNA sequences of interest on NCBI website?

You could use NCBI eutils to get all rRNA sequences and then filter out 16S sequences based on the descriptions:

esearch -db nuccore -query "txid243164[Organism:noexp] AND biomol_rrna[PROP] AND 16S"|efetch -format fasta

However, this query would fail if 16S is written as 16s in the description.

0 answers

No answers yet.

Log in to answer this question.