Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

archaea and bacteria 16S duplicate #35

Open
chloelulu opened this issue Feb 11, 2019 · 1 comment
Open

archaea and bacteria 16S duplicate #35

chloelulu opened this issue Feb 11, 2019 · 1 comment
Assignees

Comments

@chloelulu
Copy link

Hi, developer,
Thanks for creating such efficient software. I have used it to find the 16S rRNA hits in my de-novo assembled genome bins. My purpose is to search for archaea and bacteria, so I run the result separately with -k bac and -k arc.
However, the result is so confusing. For example, one of the bin found two 16S hits of archaea and also two hits of bacteria. The header of the hits are >16S_rRNA::NODE_2_length_100533_cov_5.789665:250-1687(-) and >16S_rRNA::NODE_8_length_10807_cov_5.393508:10362-10807(-) in bacteria output. The header of the hits are >16S_rRNA::NODE_2_length_100533_cov_5.789665:251-1678(-) and >16S_rRNA::NODE_8_length_10807_cov_5.393508:10363-10803(-)
And I blast both fasta hits to RDP classifier, and the archaea hits outputs are 16S_rRNA::NODE_2_length_100533_cov_5.789665:251-1678(-);+;Bacteria;100%;"Bacteroidetes";98%;"Bacteroidia";96%;"Bacteroidales";96%;"Rikenellaceae";38%;Mucinivorans;33% 16S_rRNA::NODE_8_length_10807_cov_5.393508:10363-10803(-);+;Bacteria;99%;Firmicutes;70%;Clostridia;61%;Clostridiales;61%;Ruminococcaceae;43%;Hydrogenoanaerobacterium;14%
Also bacteria hits outputs are 16S_rRNA::NODE_2_length_100533_cov_5.789665:250-1687(-);+;Bacteria;100%;"Bacteroidetes";98%;"Bacteroidia";94%;"Bacteroidales";94%;"Rikenellaceae";34%;Mucinivorans;24% 16S_rRNA::NODE_8_length_10807_cov_5.393508:10362-10807(-);+;Bacteria;99%;Firmicutes;78%;Clostridia;53%;Clostridiales;53%;Ruminococcaceae;40%;Hydrogenoanaerobacterium;14%
So my question are -
(1) The result of bacteria and archaea are the same, both are bacteria. Why they are classified into two parts, bacteria and archaea?
(2) The two hits came from one genome bin, why they can be predicted and have two 16S with different taxonomy classification?

Thanks so much for your patience!
Best.

@tseemann tseemann self-assigned this Oct 3, 2019
@tseemann
Copy link
Owner

tseemann commented Oct 3, 2019

Yes, both bacteria and archaea share a 16S model from RFAM.

NAME  16S_rRNA
ACC   RF00177

Barrnap is designed for bacterial isolates. It was not designed to predict kingdom of MAGs.

I do not know why And I blast both fasta hits to RDP classifier gives different answers for the same identical (?) sequences.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants