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

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091260‑3091349                                              |                                             tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                          <  1:184341/50‑1 (MQ=255)    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                        >  1:1824811/1‑49 (MQ=255)      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                     <  1:1424583/50‑1 (MQ=255)         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                  >  1:1325081/1‑50 (MQ=255)           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 <  1:2018406/49‑1 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:237154/1‑49 (MQ=255)                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           <  1:1599357/50‑1 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                    <  1:774050/49‑1 (MQ=255)                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  1:649439/50‑1 (MQ=255)                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                <  1:1217773/50‑1 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             <  1:1963294/50‑1 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             <  1:2319092/50‑1 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             <  1:815305/50‑1 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:19816/1‑50 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc          <  1:219444/50‑1 (MQ=255)                                    tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:302878/1‑46 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:923913/1‑44 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:1541311/50‑1 (MQ=255)                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:335703/50‑1 (MQ=255)                                              |                                             TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091260‑3091349

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