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

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091260‑3091349                                              |                                             tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            >  1:712960/1‑48 (MQ=255) tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            >  1:2275348/1‑48 (MQ=255) tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            <  1:893862/48‑1 (MQ=255) tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            >  1:361366/1‑48 (MQ=255)  ttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                           >  1:1301074/1‑48 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       <  1:2391847/49‑1 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                        >  1:529788/1‑48 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                        <  1:1683410/48‑1 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                        >  1:415116/1‑48 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                        >  1:1314100/1‑48 (MQ=255)      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                      <  1:2622651/49‑1 (MQ=255)       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       >  1:1349990/1‑47 (MQ=255)        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    <  1:795706/49‑1 (MQ=255)        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    >  1:932193/1‑49 (MQ=255)         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                   >  1:830887/1‑49 (MQ=255)         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                   <  1:621549/49‑1 (MQ=255)          aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                  >  1:897177/1‑49 (MQ=255)           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                  <  1:114759/48‑1 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               <  1:1731530/49‑1 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                >  1:1124823/1‑48 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                >  1:1926489/1‑48 (MQ=255)              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:474235/1‑48 (MQ=255)               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                             <  1:2568653/49‑1 (MQ=255)                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                            >  1:1069141/1‑49 (MQ=255)                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                            >  1:1563886/1‑49 (MQ=255)                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:734847/1‑49 (MQ=255)                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          <  1:588886/49‑1 (MQ=255)                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                       >  1:1855969/1‑49 (MQ=255)                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                       >  1:2006393/1‑49 (MQ=255)                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     >  1:1705104/1‑49 (MQ=255)                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     >  1:310286/1‑49 (MQ=255)                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     >  1:1337940/1‑49 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                      <  1:1922670/47‑1 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     <  1:52213/48‑1 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                    <  1:1461800/48‑1 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   >  1:1017780/1‑49 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   >  1:777728/1‑49 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   >  1:796313/1‑49 (MQ=255)                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                  <  1:2322179/49‑1 (MQ=255)                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                  <  1:2413402/49‑1 (MQ=255)                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  1:1221816/49‑1 (MQ=255)                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  1:1916740/49‑1 (MQ=255)                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 >  1:680305/1‑49 (MQ=255)                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 >  1:154889/1‑48 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  1:1169423/1‑47 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa              <  1:1781690/49‑1 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:484855/1‑49 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            <  1:1716388/49‑1 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:33037/1‑49 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           <  1:1486616/49‑1 (MQ=255)                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:2068605/1‑48 (MQ=255)                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:799454/1‑48 (MQ=255)                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:2520083/1‑48 (MQ=255)                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:2530351/1‑47 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa       <  1:180319/48‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa       <  1:1431290/48‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt      <  1:754683/49‑1 (MQ=255)                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:667648/1‑43 (MQ=255)                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:584696/49‑1 (MQ=255)                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:445366/49‑1 (MQ=255)                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg    <  1:1175475/49‑1 (MQ=255)                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     >  1:1520496/1‑48 (MQ=255)                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   <  1:1480979/49‑1 (MQ=255)                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   >  1:990499/1‑49 (MQ=255)                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:2225088/49‑1 (MQ=255)                                              |                                             TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091260‑3091349

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