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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351                                               |                                               gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              <  1:5195435/49‑1 (MQ=255)  tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                            <  1:4028420/50‑1 (MQ=255)  tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                            <  1:36827/50‑1 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:1407277/49‑1 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:5256492/1‑49 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:1945329/49‑1 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          <  1:5006962/49‑1 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          <  1:1606333/49‑1 (MQ=255)     gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          <  1:3114476/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:2749126/1‑49 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:1091040/1‑49 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:2870165/1‑49 (MQ=255)       cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:4127192/50‑1 (MQ=255)       cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:2031389/50‑1 (MQ=255)        ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:1628845/1‑50 (MQ=255)          tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  1:642585/1‑50 (MQ=255)          tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  1:2551926/1‑50 (MQ=255)           aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   >  1:6077361/1‑50 (MQ=255)            gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   <  1:3875591/49‑1 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:2939309/1‑49 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:5111324/1‑49 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1169834/1‑49 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1384089/1‑49 (MQ=255)             cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1602946/1‑49 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:4620353/50‑1 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:2118892/50‑1 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:4303623/50‑1 (MQ=255)               aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:1200142/50‑1 (MQ=255)               aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:2874083/50‑1 (MQ=255)               aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:3714978/50‑1 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:1964609/1‑50 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:5200451/50‑1 (MQ=255)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:271190/1‑50 (MQ=255)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             <  1:1975671/50‑1 (MQ=255)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:2374179/1‑50 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:1475206/50‑1 (MQ=255)                   ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:2437357/50‑1 (MQ=255)                    gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:3019927/1‑48 (MQ=255)                      cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:3146935/1‑50 (MQ=255)                       tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:1776863/1‑50 (MQ=255)                       tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:2719256/1‑50 (MQ=255)                         tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:1166859/48‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    >  1:4841165/1‑50 (MQ=255)                           aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:1216618/49‑1 (MQ=255)                           aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:4480759/49‑1 (MQ=255)                               cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:1524386/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:719380/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4127633/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3137461/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:6125897/1‑50 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:322523/1‑49 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5372504/1‑49 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3107202/1‑48 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:4372599/50‑1 (MQ=255)                                    ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5424676/1‑47 (MQ=255)                                    ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:4845462/50‑1 (MQ=255)                                       ccAGCGGCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:3638233/49‑1 (MQ=38)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3329549/1‑44 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       >  1:6639499/1‑49 (MQ=255)                                         aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     <  1:4776374/50‑1 (MQ=255)                                           cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg   <  1:6574629/50‑1 (MQ=255)                                           cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    >  1:2236954/1‑49 (MQ=255)                                            gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:5366387/50‑1 (MQ=255)                                               |                                               GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351

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