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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091259‑3091349                                               |                                             gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            >  1:6487453/1‑49 (MQ=255)   ttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                          >  1:4422440/1‑49 (MQ=255)    tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                         <  1:1275377/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       >  1:5776048/1‑49 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       <  1:2007235/49‑1 (MQ=255)      aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       <  1:1836740/49‑1 (MQ=255)        ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    >  1:7017385/1‑50 (MQ=255)         gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                   >  1:3468006/1‑50 (MQ=255)         gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                   <  1:2746507/50‑1 (MQ=255)          tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                  >  1:898063/1‑50 (MQ=255)          tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                  >  1:4584500/1‑50 (MQ=255)           aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 >  1:2168032/1‑50 (MQ=255)           aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 >  1:5980551/1‑50 (MQ=255)            gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 <  1:6356312/49‑1 (MQ=255)            gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 <  1:2956804/49‑1 (MQ=255)            gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                >  1:7105710/1‑50 (MQ=255)             cAACCAGGCCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                >  1:7102081/1‑49 (MQ=38)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:5435808/1‑49 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:4733760/1‑49 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:1321790/1‑49 (MQ=255)              aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:4287152/1‑49 (MQ=255)               aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:81413/1‑48 (MQ=255)               aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                             <  1:6941924/50‑1 (MQ=255)                ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                            >  1:3008974/1‑50 (MQ=255)                 cAGGTCTTTGAGGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           <  1:5184693/50‑1 (MQ=38)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:6282284/1‑50 (MQ=255)                 cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:2431391/1‑50 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:4953533/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:207424/1‑49 (MQ=255)                  aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:1816312/1‑49 (MQ=255)                   ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          >  1:2437228/1‑49 (MQ=255)                   ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                         <  1:4070354/50‑1 (MQ=255)                   ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                         <  1:2690434/50‑1 (MQ=255)                   ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                         <  1:3164826/50‑1 (MQ=255)                     tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                       <  1:3789042/50‑1 (MQ=255)                     tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                       >  1:971019/1‑50 (MQ=255)                      cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGCAAGAGAAgg                      >  1:1637473/1‑50 (MQ=38)                       tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     >  1:5417541/1‑50 (MQ=255)                        ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                    >  1:4480144/1‑50 (MQ=255)                         tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     <  1:385117/48‑1 (MQ=255)                         tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                    <  1:1874548/49‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   <  1:2118743/49‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   <  1:1057596/49‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   <  1:1173391/49‑1 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                  >  1:1573032/1‑50 (MQ=255)                          gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                  >  1:2433299/1‑50 (MQ=255)                           aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                  <  1:4322680/49‑1 (MQ=255)                           aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 >  1:4848721/1‑50 (MQ=255)                            cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                <  1:2868300/50‑1 (MQ=255)                             gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  1:2896170/1‑49 (MQ=255)                             gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  1:608579/1‑49 (MQ=255)                              ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  1:387836/1‑48 (MQ=255)                              ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  1:5911530/1‑48 (MQ=255)                              ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  1:5935199/1‑48 (MQ=255)                               cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             <  1:2410455/50‑1 (MQ=255)                               cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             >  1:950701/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:4723177/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:2395985/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:4952840/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:3370046/1‑50 (MQ=255)                                gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:2002170/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:5682324/1‑50 (MQ=255)                                 aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:3620401/1‑50 (MQ=255)                                  aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc          <  1:5887176/50‑1 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:2349107/1‑48 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:1512149/1‑48 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:232808/1‑48 (MQ=255)                                   aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:4162313/1‑48 (MQ=255)                                    ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:3073412/1‑47 (MQ=255)                                    ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg        <  1:6928518/50‑1 (MQ=255)                                     tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:85971/1‑46 (MQ=255)                                     tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:6332364/1‑46 (MQ=255)                                      aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:4691700/1‑45 (MQ=255)                                      aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa       <  1:1930769/49‑1 (MQ=255)                                      aCCAGCGCCGGC‑AAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt      <  1:2323472/49‑1 (MQ=38)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:6416439/1‑44 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt      <  1:7104922/49‑1 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:6396145/50‑1 (MQ=255)                                       ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:28669/50‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg    <  1:785894/50‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg    <  1:4077602/50‑1 (MQ=255)                                        cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     >  1:680606/1‑49 (MQ=255)                                         aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   <  1:20673/50‑1 (MQ=255)                                         aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   >  1:3510313/1‑50 (MQ=255)                                         aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   >  1:103028/1‑50 (MQ=255)                                          gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:2288984/50‑1 (MQ=255)                                          gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:4969227/50‑1 (MQ=255)                                          gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:4584984/50‑1 (MQ=255)                                          gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:3617956/50‑1 (MQ=255)                                          gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:3484134/50‑1 (MQ=255)                                               |                                             GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091259‑3091349

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