Question: which better for the prediction of 16s rRNA
0
gravatar for zhangdengwei
5 weeks ago by
zhangdengwei60
zhangdengwei60 wrote:

Hi all,

Anyone can recommend a better software for the prediction of 16S rRNA? I used barrnap to do this, and the length of 16S rRNA varied a lot, from ~500 bp to ~1400 bp. Of course, the longer, the better. So is there any alternatives to 16S prediction. Any suggestions would be greatly appreciated !!!

ADD COMMENTlink modified 5 weeks ago by Mensur Dlakic5.8k • written 5 weeks ago by zhangdengwei60
0
gravatar for Mensur Dlakic
5 weeks ago by
Mensur Dlakic5.8k
USA
Mensur Dlakic5.8k wrote:

There is nothing wrong with having different 16S rRNA lengths, assuming you are predicting from a metagenomic sample. That simply means that some of your contigs contain incomplete rRNA sequences - no program can predict a complete rRNA if it isn't there. As long as you have the latest version, barrnap is a safe choice.

If you want to convince yourself that barrnap predictions are solid, another program to try is RNAmmer. It has similar speed and sensitivity to barrnap, though the latter is more likely to be up to date. Lastly, if you want to go with heavy artillery, your choice would be SSU-ALIGN. It will have the highest sensitivity because it uses covariance models for scoring, but that's also why it will be slowest by couple orders of magnitude. If you have say 100+ Mb of sequence data, it will take many hours (even a whole day) to complete the prediction, and that's assuming the job is spread over many CPUs. You can give it a try if your dataset is relatively small or if you have time to spare, but keep in mind that it won't give you longer sequences if they are not there. The most you could expect is that it may detect some borderline matches, but usually the result will end up being the same as with barrnap.

ADD COMMENTlink written 5 weeks ago by Mensur Dlakic5.8k

Many thanks. Your comment is useful and helpful to me. Maybe I can have RNAmmer a try, to find whether a difference exists between barrnap and RNAmmer.

ADD REPLYlink written 5 weeks ago by zhangdengwei60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1191 users visited in the last hour