Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,607,016	A→G	S229P (TCT→CCT)	rsmC ←	16S rRNA m(2)G1207 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,607,016	0	A	G	100.0%	184.0 / NA	53	S229P (TCT→CCT)	rsmC	16S rRNA m(2)G1207 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/21);  total (32/21)

AGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAG  >  NC_000913/4606938‑4607089                                                                               |                                                                          aGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAggga                                                                         >  2:956810/1‑81 (MQ=255)  gACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAggga                                                                         >  1:956822/1‑80 (MQ=255)  gACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAggga                                                                         >  1:956819/1‑80 (MQ=255)  gACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGAc                                                                        <  2:956750/81‑1 (MQ=255)   aCGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACa                                                                       >  1:956804/1‑81 (MQ=255)    cGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACAt                                                                      >  2:956807/1‑81 (MQ=255)        gCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGc                                                                   <  2:956752/80‑1 (MQ=255)        gCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGATGGATTCTACCGCCGGGGCAGGGACATCGc                                                                   <  2:956751/80‑1 (MQ=255)         cTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGcac                                                                 <  1:956753/81‑1 (MQ=255)         cTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCCGGGCAGGGACATCGca                                                                  >  2:956809/1‑80 (MQ=255)            gCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACaa                                                               <  2:956754/80‑1 (MQ=255)             cAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAgg                                                             >  1:956813/1‑81 (MQ=255)              aaaGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGt                                                            >  2:956838/1‑81 (MQ=255)                aGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGt                                                            >  1:956829/1‑79 (MQ=255)                aGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGt                                                            >  1:956811/1‑79 (MQ=255)                aGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTga                                                          <  2:956755/81‑1 (MQ=255)                  aCTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTgaga                                                        <  2:956756/81‑1 (MQ=255)                  aCTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTgaga                                                        >  2:956808/1‑81 (MQ=255)                    ttCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGAc                                                       >  2:956816/1‑80 (MQ=255)                     tCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACg                                                      >  1:956797/1‑80 (MQ=255)                      cACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACg                                                      <  1:956757/79‑1 (MQ=255)                       aCCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAAt                                                   <  2:956758/81‑1 (MQ=255)                        ccTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAAtt                                                  >  1:956846/1‑81 (MQ=255)                         cTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAAttt                                                 <  1:956759/81‑1 (MQ=255)                           tCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTc                                               >  1:956840/1‑81 (MQ=255)                            cAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCg                                              >  2:956826/1‑81 (MQ=255)                             aaCACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCgg                                             >  1:956833/1‑81 (MQ=255)                              acacCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGc                                            <  1:956760/81‑1 (MQ=255)                               cacCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGc                                            >  2:956835/1‑80 (MQ=255)                                   gTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATg                                       >  1:956828/1‑81 (MQ=255)                                    ttGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATgg                                      <  2:956761/81‑1 (MQ=255)                                     tGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGgc                                     >  1:956815/1‑81 (MQ=255)                                      gggcGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGgc                                     <  1:956762/77‑1 (MQ=255)                                       gCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGgcg                                    <  1:956763/80‑1 (MQ=255)                                           cAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCa                                >  2:956841/1‑80 (MQ=255)                                               tgtTGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAgg                            >  1:956836/1‑80 (MQ=255)                                                 ttGCGCGGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCaa                         >  2:956848/1‑81 (MQ=255)                                                     gcgGCTGGCTTGTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGGAAAGGCAACTg                      >  2:956771/1‑80 (MQ=255)                                                     gcgGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTg                      >  1:956821/1‑80 (MQ=255)                                                     gcgGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTg                      >  2:956818/1‑80 (MQ=255)                                                     gcgGCTGGCTTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGa                     <  1:956764/81‑1 (MQ=255)                                                        gctggctTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAg                  >  1:956857/1‑81 (MQ=255)                                                        gctggctTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAg                  >  2:956845/1‑81 (MQ=255)                                                           ggctTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCAc               >  2:956849/1‑81 (MQ=255)                                                           ggctTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCAc               >  1:956832/1‑81 (MQ=255)                                                           ggctTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCAc               >  1:956839/1‑81 (MQ=255)                                                            gctTCTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCAc               <  1:956765/80‑1 (MQ=255)                                                                cTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGccc           <  2:956766/80‑1 (MQ=255)                                                                cTACCGCCGGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGcc            <  1:956767/79‑1 (MQ=255)                                                                   ccgccgGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCgcgc       <  1:956768/81‑1 (MQ=255)                                                                   ccgccgGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCgcgc       >  1:956831/1‑81 (MQ=255)                                                                     gccgGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCca     >  2:956869/1‑81 (MQ=255)                                                                      ccgGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCcac    <  2:956770/81‑1 (MQ=255)                                                                      ccgGGGCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCcac    <  1:956769/81‑1 (MQ=255)                                                                        ggggCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAg  >  2:956847/1‑81 (MQ=255)                                                                        ggggCAGGGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAg  >  1:956843/1‑81 (MQ=255)                                                                               |                                                                          AGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAG  >  NC_000913/4606938‑4607089

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑