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

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091260‑3091351                                              |                                               tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                            <  1:8091745/50‑1 (MQ=255) tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                            <  1:7573058/50‑1 (MQ=255)     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:4992925/49‑1 (MQ=255)       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:4467290/1‑50 (MQ=255)       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:5567613/50‑1 (MQ=255)       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:5697219/1‑50 (MQ=255)         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  1:4527219/1‑50 (MQ=255)           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   <  1:7728968/49‑1 (MQ=255)           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1744348/1‑50 (MQ=255)            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  1:952985/50‑1 (MQ=255)            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:563063/1‑49 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:1633088/1‑49 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:5283144/50‑1 (MQ=255)             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:4568560/50‑1 (MQ=255)              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:1102095/50‑1 (MQ=255)              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:3848248/50‑1 (MQ=255)               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:2800038/1‑50 (MQ=255)               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:4002735/50‑1 (MQ=255)               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:3711592/1‑50 (MQ=255)                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:6972041/1‑50 (MQ=255)                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:7517596/1‑50 (MQ=255)                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:105797/50‑1 (MQ=255)                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:5255196/50‑1 (MQ=255)                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:3869993/1‑48 (MQ=255)                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:1328967/1‑48 (MQ=255)                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                          <  1:2327790/50‑1 (MQ=255)                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                          <  1:6742550/50‑1 (MQ=255)                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:5337983/50‑1 (MQ=255)                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:5526984/50‑1 (MQ=255)                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:1077053/1‑50 (MQ=255)                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:5789326/1‑50 (MQ=255)                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:4985736/1‑50 (MQ=255)                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:5912676/1‑50 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:808890/48‑1 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:6673896/48‑1 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:5461561/48‑1 (MQ=255)                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:4364183/49‑1 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     <  1:8194141/49‑1 (MQ=255)                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     <  1:3468072/49‑1 (MQ=255)                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:1501735/50‑1 (MQ=255)                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  1:4783317/50‑1 (MQ=255)                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa                <  1:6637452/50‑1 (MQ=255)                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa                <  1:3056259/50‑1 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:6273403/1‑50 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:6772835/1‑50 (MQ=255)                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:1700840/1‑50 (MQ=255)                              cGAAAGTAACAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               <  1:2810432/50‑1 (MQ=38)                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  1:2557546/50‑1 (MQ=255)                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  1:5409127/50‑1 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:7842570/1‑50 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:658664/1‑50 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             <  1:3414927/50‑1 (MQ=255)                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3780384/1‑50 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:7279370/50‑1 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:3775067/50‑1 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:6258792/1‑49 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:6184444/1‑49 (MQ=255)                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:206738/1‑49 (MQ=255)                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:5646145/50‑1 (MQ=255)                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:2099327/50‑1 (MQ=255)                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:4406890/50‑1 (MQ=255)                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:1861517/1‑45 (MQ=255)                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:8164116/1‑45 (MQ=255)                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:391444/1‑45 (MQ=255)                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:6184978/49‑1 (MQ=255)                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:7649935/49‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:8212564/1‑44 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:676245/49‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:4236612/49‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:6816013/49‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:7020634/50‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:855220/50‑1 (MQ=255)                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:1464791/50‑1 (MQ=255)                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:3606291/50‑1 (MQ=255)                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  1:3064027/1‑50 (MQ=255)                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     <  1:7229812/50‑1 (MQ=255)                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:5121744/50‑1 (MQ=255)                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:5934817/50‑1 (MQ=255)                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg   <  1:7984217/50‑1 (MQ=255)                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    >  1:3306080/1‑49 (MQ=255)                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    >  1:1335393/1‑49 (MQ=255)                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    >  1:3882223/1‑49 (MQ=255)                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:5220694/50‑1 (MQ=255)                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:6763838/1‑50 (MQ=255)                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:6308888/50‑1 (MQ=255)                                              |                                               TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091260‑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: ‑