# silva database question

**URL:** <https://forum.qiime2.org/t/silva-database-question/18904>\
**Category:** General Discussion\
**Created:** [March 18, 2021, 1:35am UTC](https://forum.qiime2.org/t/silva-database-question/18904 "2021-03-18T01:35:46Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![kimshinseung](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/k/e274bd/32.png) [@kimshinseung](https://forum.qiime2.org/u/kimshinseung)\
**Post date:** [March 18, 2021, 1:35am UTC](https://forum.qiime2.org/t/silva-database-question/18904/1 "2021-03-18T01:35:47Z")

</div>

Hi

I am wondering why the silva database is so long. There are over 2000 bacteria, so why?

And I am wondering if there are no start and stop codons in 16srRNA.

Thank you.

---

<div class="post-metadata">

**Author:** ![colinbrislawn](https://forum.qiime2.org/user_avatar/forum.qiime2.org/colinbrislawn/32/6221_2.png) [@colinbrislawn](https://forum.qiime2.org/u/colinbrislawn)\
**Post date:** [March 18, 2021, 3:44pm UTC](https://forum.qiime2.org/t/silva-database-question/18904/2 "2021-03-18T15:44:16Z")

</div>

Hello @kimshinseung,

> [@kimshinseung](#):
>
> I am wondering why the silva database is so long. There are over 2000 bacteria, so why?

The newest SILVA release when I write this is [v138.1](https://www.arb-silva.de/documentation/release-1381/). The SSU Ref contains 2 million reads. After clustering this at 99% identity to make SSU Ref NR 99, there are still **510,508 reads**.

That's just a lot of bacteria! That's probably why the database is so big... There are even more reads in the LSU Ref.

> [@kimshinseung](#):
>
> And I am wondering if there are no start and stop codons in 16srRNA.

Start and stop codons show where protein-coding genes will be translated into proteins in messenger-RNA. Ribosomal-RNA is part of the ribosome and is not protein-coding, so rRNA does not have any codons at all!

---

<div class="post-metadata">

**Author:** ![kimshinseung](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/k/e274bd/32.png) [@kimshinseung](https://forum.qiime2.org/u/kimshinseung)\
**Post date:** [March 19, 2021, 1:52am UTC](https://forum.qiime2.org/t/silva-database-question/18904/3 "2021-03-19T01:52:32Z")

</div>

Thank you for answer.

The first question was about the length, not the number of bacteria in the silva database.

The 16s rrna is about 1500 bp, but most of the silva is over 2000 bp.

Why is it so long?

Thank you.

---

<div class="post-metadata">

**Author:** ![colinbrislawn](https://forum.qiime2.org/user_avatar/forum.qiime2.org/colinbrislawn/32/6221_2.png) [@colinbrislawn](https://forum.qiime2.org/u/colinbrislawn)\
**Post date:** [March 19, 2021, 3:53am UTC](https://forum.qiime2.org/t/silva-database-question/18904/4 "2021-03-19T03:53:37Z")

</div>

I think I'm a little confused...

Figure 1 from the [2007 Silva paper](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2175337/) show that most sequences in SSU Ref are _under_ 1500 in length:  
[https://oup.silverchair-cdn.com/oup/backfile/Content\_public/Journal/nar/35/21/10.1093\_nar\_gkm864/2/gkm864f1.jpeg?Expires=1619149446&Signature=TL6U8rLB8uxNCWrkV6hKTpNtmLwgsjwc10ClJYKjlh0f04hy-yZhvuYry8gkPklYjORKQEuViBme3USXZT8gFL~jn9gqhfE7Q4l3XO2aDXdgeS71rMQr9ogFSm6AKIS085x~jWOQ~oJ9nmU7zes2oAKtaO4cu6YtrtDWn6ecOPf7dpvx5ZFmQe7bNRPeFo~cqFsQPcWqtrfqjGY-xqjobLmfhmrcDCKAW3emytuAzbEtN1q0dOWDwkiAPAPALe4Ueug2r45nYLo57aWPVkcz84LonGr04n1KyS4mBl3E-Ei2GGLPfglcsYZPnrP0AQaSdyN13MVPZH6F4IdY9Z3jcQ\_\_&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA](https://oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/nar/35/21/10.1093_nar_gkm864/2/gkm864f1.jpeg?Expires=1619149446&Signature=TL6U8rLB8uxNCWrkV6hKTpNtmLwgsjwc10ClJYKjlh0f04hy-yZhvuYry8gkPklYjORKQEuViBme3USXZT8gFL~jn9gqhfE7Q4l3XO2aDXdgeS71rMQr9ogFSm6AKIS085x~jWOQ~oJ9nmU7zes2oAKtaO4cu6YtrtDWn6ecOPf7dpvx5ZFmQe7bNRPeFo~cqFsQPcWqtrfqjGY-xqjobLmfhmrcDCKAW3emytuAzbEtN1q0dOWDwkiAPAPALe4Ueug2r45nYLo57aWPVkcz84LonGr04n1KyS4mBl3E-Ei2GGLPfglcsYZPnrP0AQaSdyN13MVPZH6F4IdY9Z3jcQ__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)

Sequences in the LSU Ref database are longer, but I'm not sure what database you are using...

What is the median length you are seeing in Silva? I'm not an expert on the Silva database, so I could be missing something!

Let me know what you see.

---

<div class="post-metadata">

**Author:** ![kimshinseung](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/k/e274bd/32.png) [@kimshinseung](https://forum.qiime2.org/u/kimshinseung)\
**Post date:** [March 19, 2021, 8:18am UTC](https://forum.qiime2.org/t/silva-database-question/18904/5 "2021-03-19T08:18:22Z")

</div>

Thank you for answer.

I have extracted the 2 qza files downloaded here.

qiime rescript get-silva-data   
--p-version '138'   
--p-target 'SSURef\_NR99'   
--p-include-species-labels   
--o-silva-sequences silva-138-ssu-nr99-seqs.qza   
--o-silva-taxonomy silva-138-ssu-nr99-tax.qza

Would there be a problem if I write this data as it is?

Thank you.

---

<div class="post-metadata">

**Author:** ![SoilRotifer](https://forum.qiime2.org/user_avatar/forum.qiime2.org/soilrotifer/32/21071_2.png) [@SoilRotifer](https://forum.qiime2.org/u/SoilRotifer)\
**Post date:** [March 19, 2021, 4:16pm UTC](https://forum.qiime2.org/t/silva-database-question/18904/6 "2021-03-19T16:16:21Z")

</div>

> [@kimshinseung](#):
>
> The 16s rrna is about 1500 bp, but most of the silva is over 2000 bp.

Remember the SILVA SSU contains not only Bacterial, Acrachaeal, and Eukaryal 16S rRNA gene sequences, but also contains the Eukaryal 18S rRNA gene sequences too (which are longer). That is, the 18S rRNA gene is a cystolic homologue of the 16S rRNA gene.

-Mike

---

<div class="post-metadata">

**Author:** ![system](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/2/1/21af5fe23cb6f4579467c66a9ed94e55274ca7bd.svg) [@system](https://forum.qiime2.org/u/system)\
**Post date:** [April 19, 2021, 10:16pm UTC](https://forum.qiime2.org/t/silva-database-question/18904/7 "2021-04-19T22:16:40Z")

</div>

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.
