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%	126.2 / NA	33	S58G (AGC→GGC)	rsmE	16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/16);  total (17/16)

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091260‑3091351                                              |                                               tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             >  1:2395689/1‑49 (MQ=255)    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:2325498/1‑49 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:1577135/50‑1 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:1765388/50‑1 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:1109942/49‑1 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:1703167/1‑49 (MQ=255)      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:2397092/50‑1 (MQ=255)       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:158622/50‑1 (MQ=255)         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    <  1:2482064/50‑1 (MQ=255)            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  1:655022/50‑1 (MQ=255)            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:854360/1‑49 (MQ=255)            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:2738478/1‑49 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:2562417/1‑49 (MQ=255)              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:1790478/1‑48 (MQ=255)              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:1776817/1‑48 (MQ=255)                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             <  1:2209250/50‑1 (MQ=255)                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:1846915/1‑49 (MQ=255)                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:1031106/1‑49 (MQ=255)                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:1492804/1‑49 (MQ=255)                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:1818115/50‑1 (MQ=255)                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:102650/1‑50 (MQ=255)                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:638630/1‑50 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    >  1:1728054/1‑50 (MQ=255)                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  1:18765/50‑1 (MQ=255)                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa                <  1:2014914/50‑1 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:2516630/50‑1 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:213967/50‑1 (MQ=255)                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:1304935/1‑47 (MQ=255)                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:2799297/50‑1 (MQ=255)                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:879825/50‑1 (MQ=255)                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:1737837/50‑1 (MQ=255)                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:2287654/1‑50 (MQ=255)                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:862883/1‑50 (MQ=255)                                              |                                               TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091260‑3091351

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