# Silva.bacteria.fasta file

**URL:** <https://forum.mothur.org/t/silva-bacteria-fasta-file/22341>\
**Category:** Commands in mothur\
**Created:** [June 17, 2025, 4:37pm UTC](https://forum.mothur.org/t/silva-bacteria-fasta-file/22341 "2025-06-17T16:37:15Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![James146](https://avatars.discourse-cdn.com/v4/letter/j/6f9a4e/32.png) [@James146](https://forum.mothur.org/u/James146)\
**Post date:** [June 17, 2025, 4:37pm UTC](https://forum.mothur.org/t/silva-bacteria-fasta-file/22341/1 "2025-06-17T16:37:15Z")

</div>

Hi

I’m struggling to obtain the newest version of the SILVA reference file (v 138.2) from the [SILVA reference files](https://mothur.org/wiki/silva_reference_files/#release-1381) link. I can’t download ARB to my PC. Has anyone managed to follow the README file to obtain the silva.bacteria.fasta file required to run the MiSeq SOP?

Is there any reason why each new silva.bacteria.fasta file with the updated SILVA version isn’t readily provided (like v102 is provided at the top of the MiSeq SOP workflow)? Perhaps they are and I just missed it…?

Thanks in advance for any help!

James

---

<div class="post-metadata">

**Author:** ![pschloss](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/pschloss/32/4_2.png) [@pschloss](https://forum.mothur.org/u/pschloss)\
**Post date:** [June 24, 2025, 12:50pm UTC](https://forum.mothur.org/t/silva-bacteria-fasta-file/22341/2 "2025-06-24T12:50:23Z")

</div>

Hi James,

There’s no need for you to use ARB unless you’re doing some type of further customization. The links at [Silva reference files](https://mothur.org/wiki/silva_reference_files/) are the fasta files you need to run the alignment. They are what the code in the README files generate. Does this help?

Pat

---

<div class="post-metadata">

**Author:** ![James146](https://avatars.discourse-cdn.com/v4/letter/j/6f9a4e/32.png) [@James146](https://forum.mothur.org/u/James146)\
**Post date:** [June 24, 2025, 3:18pm UTC](https://forum.mothur.org/t/silva-bacteria-fasta-file/22341/3 "2025-06-24T15:18:59Z")

</div>

Hi Pat,

Thanks for your response.

I’m struggling to get my head around the difference between the **silva.bacteria.fasta** file from “[SILVA-based bacterial reference alignment (v102)](https://mothur.s3.us-east-2.amazonaws.com/wiki/silva.bacteria.zip)” and the [release 138.2 files](https://mothur.s3.us-east-2.amazonaws.com/wiki/silva.seed_v138_2.tgz).

The release 138.2 download gives **silva.seed\_v138\_2.align** and **silva.seed\_v138\_2.tax**.

Is the **silva.bacteria.fasta** used as the example in the SOP a reformatted/trimmed version of a larger **silva.seed\_vxyz.align** file? Therefore, to obtain the V4 region for v138.2, must I preprocess the **silva.seed\_v138\_2.align** file, or can I directly substitute it for the **silva.bacteria.fasta** file in the `pcr.seqs` command?

E.g. should I now run:  
`mothur > pcr.seqs(fasta=silva.seed_v138_2.align, start=11895, end=25318, keepdots=F)`  
compared to my previous work where I ran:  
`mothur > pcr.seqs(fasta=silva.bacteria.fasta, start=11895, end=25318, keepdots=F)`  
(former v138.2, latter v102)  
Should the start/end coordinates be changed for v138.2?

Then, later on we run `mothur > remove.lineage( ... taxon=Chloroplast-Mitochondria-unknown-Archaea-Eukaryota)` to remove the non bacterial stuff from the aligned and trimmed **silva.seed\_v138\_2.align**.

James

EDIT -  
I have now run:  
`pcr.seqs(fasta=sequence.fasta, oligos=pcrTest.oligos)`  
`align.seqs(fasta=sequence.pcr.fasta, reference=silva.seed_v138_2.align)`  
`summary.seqs(fasta=sequence.pcr.align)`  
using an _E. coli_ sequence and primers 515f (GTGCCAGCMGCCGCGGTAA) and 806r (GGACTACHVGGGTWTCTAAT). The coordinates return as 13862 and 23444, as per the SOP when using silva.bacteria.fasta, so I will keep the coordinates as the original 11895 and 25318 for consistency.  
My question still stands as to whether I directly substitute **silva.bacteria.fasta** for **silva.seed\_v138\_2.align** in `pcr.seqs` - I now think yes and that I am overthinking it, but confirmation would be great - I’ll get there in the end!

Thanks,  
James

---

<div class="post-metadata">

**Author:** ![pschloss](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/pschloss/32/4_2.png) [@pschloss](https://forum.mothur.org/u/pschloss)\
**Post date:** [June 25, 2025, 12:46pm UTC](https://forum.mothur.org/t/silva-bacteria-fasta-file/22341/4 "2025-06-25T12:46:44Z")

</div>

Hi James,

Yeah, the silva.\*.align file is a drop in replacement for silva.bacteria.fasta. Initially we separated the three domains into separate files. But now they’re all in one file. I’d suggest not using remove.lineage on the silva files. Instead, run it as shown in the SOP. Also, we find that using the numeric coordinates works better than the primers in pcr.seqs since no primer set is 100% for all of the sequences in the database.

Hope this helps!  
Pat

---

<div class="post-metadata">

**Author:** ![system](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/system/32/2_2.png) [@system](https://forum.mothur.org/u/system)\
**Post date:** [July 5, 2025, 12:47pm UTC](https://forum.mothur.org/t/silva-bacteria-fasta-file/22341/5 "2025-07-05T12:47:00Z")

</div>

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