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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACG  >  NC_000913/3091259‑3091350                                               |                                              gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                             <  1:102270/49‑1 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                             >  1:2408163/1‑49 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                             >  1:2071646/1‑49 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                          >  1:350289/1‑49 (MQ=255)        ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                     <  1:533279/50‑1 (MQ=255)            gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                 <  1:2362181/50‑1 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                <  1:2378335/50‑1 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                            >  1:535112/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                           <  1:1517369/50‑1 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                           <  1:1936001/50‑1 (MQ=255)                       tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                      >  1:2062122/1‑50 (MQ=255)                        ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                     >  1:1744563/1‑50 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                   >  1:2591160/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg             >  1:1296872/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg             <  1:517583/50‑1 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc            >  1:739396/1‑50 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa          <  1:953300/50‑1 (MQ=255)                                      aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa        <  1:931103/49‑1 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt       <  1:1256382/49‑1 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc      <  1:352425/50‑1 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc      <  1:1402539/50‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg     <  1:1549165/50‑1 (MQ=255)                                         aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga    >  1:604819/1‑50 (MQ=255)                                           cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg  <  1:556942/50‑1 (MQ=255)                                               |                                              GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACG  >  NC_000913/3091259‑3091350

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