Question: Extracting an information about conservation of 16S rRNA genes of bacteria from multiple sequence alignment
gravatar for Peter Leontev
5.5 years ago by
Peter Leontev90 wrote:


I am trying to get any useful information from 16S rRNA bacterial genes multiple sequence alignment (MSA), i.e. for sake of simplicity, average length of conservative and non-conservative regions. I think it may be absolutely incorrect to build MSA using, for example, Clustal, MAFFT or T-Coffee to further look at length of those regions and make some assumptions. Could you please tell me about some approach one should try?

To be more clear, I am trying to extract some info from these MSAs to use it during assembly of 16S rRNA from the NGS reads (potentially, even from metagenomic reads).

I am also looking for some conservation estimates of 16S rRNA genes and their dependence on evolutionary distance in relation to different species. For example, if we consider 16S rRNA genes MSA  from only one genus G then that estimate will have different value than if we deal with the 16S rRNA genes from G and some other genus. What potential ways exist to explore this topic?

Thank you. 



ADD COMMENTlink written 5.5 years ago by Peter Leontev90

I kinda feel the same way -- i think that knowing the (estimated) length of the conserved regions vs. the hypervariable regions should result, for example, in a more straightforward and faster way to look at alignments between reads and the gene sequences.

ADD REPLYlink written 5.2 years ago by Lynxoid220
Please log in to add an answer.


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