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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351                                               |                                               gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              >  1:5935283/1‑49 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              <  1:4906576/49‑1 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              >  1:5389332/1‑49 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:895351/49‑1 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:4974249/1‑49 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:2823484/1‑49 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:843854/1‑49 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:4495083/1‑49 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:5014512/50‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:208273/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:4164901/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:127521/1‑49 (MQ=255)       cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:5685206/1‑48 (MQ=255)        ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:1797312/1‑50 (MQ=255)        ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:2804857/1‑50 (MQ=255)          tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    <  1:2913728/50‑1 (MQ=255)          tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  1:3891037/1‑50 (MQ=255)           aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   >  1:883991/1‑50 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:2278342/1‑49 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:3584130/1‑50 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:2534225/1‑50 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:1582393/50‑1 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:160159/50‑1 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:5009845/1‑50 (MQ=255)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:2195885/1‑50 (MQ=255)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:1033040/1‑50 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:4943692/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:3616962/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:179529/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:1286450/50‑1 (MQ=255)                   ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:4969747/1‑49 (MQ=255)                     tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:869427/50‑1 (MQ=255)                     tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         >  1:1561544/1‑50 (MQ=255)                     tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:4419460/50‑1 (MQ=255)                      cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:3546739/1‑50 (MQ=255)                       tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:847008/1‑50 (MQ=255)                       tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:4541525/1‑50 (MQ=255)                        ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  1:1318009/1‑50 (MQ=255)                         tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:795340/48‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    >  1:1521594/1‑50 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    >  1:5892695/1‑50 (MQ=255)                           aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:5806563/50‑1 (MQ=255)                           aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:2935259/50‑1 (MQ=255)                            cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  1:5442551/50‑1 (MQ=255)                             gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  1:5348470/50‑1 (MQ=255)                               cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               <  1:3279799/50‑1 (MQ=255)                               cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:593430/1‑50 (MQ=255)                               cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:1632310/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:2163551/1‑50 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:1104733/50‑1 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5365909/1‑49 (MQ=255)                                   aTTACCAGCGCCGGCCAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:4683376/50‑1 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:1816801/50‑1 (MQ=255)                                    ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4184768/1‑47 (MQ=255)                                     tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:2036541/1‑46 (MQ=255)                                     tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3217810/1‑46 (MQ=255)                                      aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:2703061/1‑45 (MQ=255)                                      aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:232283/1‑45 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:3939214/49‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:5653212/50‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:2532816/50‑1 (MQ=255)                                            gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:2353024/1‑50 (MQ=255)                                               |                                               GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351

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: ‑