Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,091,305	A→G	S58G (AGC→GGC)	rsmE →	16S rRNA m(3)U1498 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,091,305	0	A	G	100.0%	151.0 / NA	38	S58G (AGC→GGC)	rsmE	16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/18);  total (20/18)

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCG  >  NC_000913/3091259‑3091347                                               |                                           gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                          <  1:703768/49‑1 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                          >  1:226164/1‑49 (MQ=255)  tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                        <  1:3222010/50‑1 (MQ=255)  tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                         <  1:1023511/49‑1 (MQ=255)  tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                         >  1:1923441/1‑49 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                       <  1:2722156/49‑1 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                      <  1:1102247/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                     <  1:35318/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                     >  1:4203101/1‑49 (MQ=255)        ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                  >  1:840261/1‑50 (MQ=255)         gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                 <  1:1207372/50‑1 (MQ=255)         gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                 >  1:3586211/1‑50 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                             <  1:2334266/50‑1 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                             <  1:2663101/50‑1 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                             >  1:822243/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                         >  1:752176/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                        <  1:2720887/50‑1 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                        <  1:3760115/50‑1 (MQ=255)                    gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                        >  1:1191827/1‑48 (MQ=255)                    gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                      <  1:736900/50‑1 (MQ=255)                         tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                   <  1:2939067/48‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                 <  1:2115193/49‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                 <  1:504545/49‑1 (MQ=255)                            cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg              <  1:1092498/50‑1 (MQ=255)                             gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg              >  1:3341423/1‑49 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg          >  1:2880528/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:2683805/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:2268414/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:1882975/1‑50 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:548350/1‑49 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:623782/1‑49 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:3006745/1‑48 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:1037210/1‑48 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa       <  1:2659826/50‑1 (MQ=255)                                     tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:2766138/1‑46 (MQ=255)                                     tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc         >  1:811371/1‑46 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg  <  1:291195/50‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc   >  1:2168806/1‑49 (MQ=255)                                               |                                           GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCG  >  NC_000913/3091259‑3091347

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