How To Identify The 16S Rrna Gene Sequence In A Genome
3
0
Entering edit mode
11.4 years ago
Ontheway ▴ 10

Hi, i have a question about the identification of 16S rRNA gene. One of the common methods is to resequence the 16S rRNA gene and the other is to search the database, such as greengenes. However, if i only have the genome sequence of a bacterial strain and i don't want to resequence the gene, how could i identify the 16S rRNA gene sequence of this strain quickly?

• 14k views
ADD COMMENT
0
Entering edit mode

"Only" have the genome sequence?!

ADD REPLY
6
Entering edit mode
11.4 years ago
cts ★ 1.7k

Use RNAmmer either through their web server or by downloading the code and running it locally

ADD COMMENT
3
Entering edit mode
5.4 years ago

For people finding this now, there is also Barrnap . Not sure how it compares to RNAmmer but it's super fast and easy to use from a command line.

ADD COMMENT
0
Entering edit mode

@gabriel.jabud Just wondering, for Barrnap and RNAmmer, is the required input an aligned genome? The barrnap website says "fasta-formatted file with the genome sequence" and RNAmmer website says "full genome sequences". I just wanted to verify, if I had raw read files, this would not be an appropriate input to these programs, right? I would have to align these raw read files first? Thank you.

ADD REPLY
0
Entering edit mode

It's also worth mentioning Barrnap is in the bacterial annotation system Prokka.

Prokka can be installed on Debian systems with sudo apt install prokka to make your life nice and easy.

ADD REPLY
1
Entering edit mode
5.4 years ago
Mensur Dlakic ★ 28k

There is a long history of identifying "common" non-coding RNAs (tRNA and rRNA) and it can be done very reliably from sequence. In addition to the two excellent suggestions above (RNAmmer and Barrnap), you may want to try SSU-align:

http://eddylab.org/software/ssu-align/

This program will not only score small-subunit rRNAs from all three kingdoms, but it will also align them.

ADD COMMENT

Login before adding your answer.

Traffic: 2435 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6