@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.
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?
3 answers
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.
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.
Log in to answer this question.
"Only" have the genome sequence?!