no rRNA sequences were detected in an almost complete genome
Closed this issue · 2 comments
Hi,
I am with great doubt with the accuracy of "RiboDetector" as a new user since I made a test for my metagenomic-assembly -genomes (MAGs). Several MAG from my previous binning and dRep work with the completeness of >95 and contaminant of <5 were selected for the rRNA detection test in "RiboDetector", the scripts are as follow (taking "MAG2809_RMAG1998.fa" as an example):
ribodetector_cpu -t 20 -l 150 -i /data2/yjk/Bison_meta/rRNA/bin_for_test/MAG2809_RMAG1998.fa -r /data2/yjk/Bison_meta/rRNA/bin_for_test/rRNA/MAG2809_RMAG1998_rRNA.fa -e both -o /data2/yjk/Bison_meta/rRNA/bin_for_test/nonrRNA/MAG2809_RMAG1998_nonrRNA.fa 1>/data2/yjk/Bison_meta/rRNA/bin_for_test/MAG2809_RMAG1998.o 2>/data2/yjk/Bison_meta/rRNA/bin_for_test/MAG2809_RMAG1998.e
the error file indicated no rRNA sequence were detected
and the output file is empty too.
Could any give me an interpretation? Many thanks in advance!
Best regards,
Jiangkun Yu
Hi Jiangkun Yu,
Thank you for your interest in RiboDetector. Our method is not designed for rRNA gene annotation for genomes. It is usually for identifying rRNA reads from sequencing data. To annotate rRNA gene in genome, you can try barrnap https://github.com/tseemann/barrnap. I hope this helps.
Best wishes,
Zhi-Luo