Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,882,059	A→C	I88S (ATT→AGT)	dnaN ←	DNA polymerase III, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,882,059	0	A	C	100.0%	178.8 / NA	47	I88S (ATT→AGT)	dnaN	DNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/28);  total (19/28)

ACGCTTCATCGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCAATTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTA  >  NC_000913/3881911‑3882202                                                                                                                                                     |                                                                                                                                                aCGCTTCATCGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACtt                                                                                                                                               <  1:337926/151‑1 (MQ=255)      tCATCGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTcc                                                                                                                                            <  1:394371/149‑1 (MQ=255)       cATCGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCcgcg                                                                                                                                         >  1:585863/1‑151 (MQ=255)          cGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGcc                                                                                                                                       <  1:239306/150‑1 (MQ=255)           gTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCtt                                                                                                                                     >  2:157825/1‑151 (MQ=255)              gCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCa                                                                                                                                   <  1:162432/150‑1 (MQ=255)                 tGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGc                                                                                                                                >  2:149693/1‑150 (MQ=255)                  gCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCa                                                                                                                               <  1:644826/150‑1 (MQ=255)                   cGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGc                                                                                                                                >  2:251381/1‑148 (MQ=255)                        gggTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACg                                                                                                                        <  1:579258/151‑1 (MQ=255)                              aTTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGat                                                                                                                  <  2:320757/151‑1 (MQ=255)                               ttCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGat                                                                                                                  >  1:355426/1‑150 (MQ=255)                                tCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGatat                                                                                                                >  1:89412/1‑151 (MQ=255)                                tCGACTGCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGat                                                                                                                  <  1:780258/149‑1 (MQ=255)                                           ctGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTg                                                                                                     >  1:158865/1‑151 (MQ=255)                                            tGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTgc                                                                                                    >  1:366963/1‑151 (MQ=255)                                             gCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTgc                                                                                                    <  2:281267/150‑1 (MQ=255)                                             gCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTg                                                                                                     <  1:623112/149‑1 (MQ=255)                                               cAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTgcgcg                                                                                                 <  1:89705/151‑1 (MQ=255)                                                       gAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACg                                                                                         >  1:690658/1‑151 (MQ=255)                                                         ggTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGgt                                                                                       <  2:70672/151‑1 (MQ=255)                                                              gggAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGtcgtcg                                                                                  <  2:109028/151‑1 (MQ=255)                                                                  aaTCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTcc                                                                              <  1:664397/151‑1 (MQ=255)                                                                             ggCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGtgtg                                                                   >  2:515759/1‑151 (MQ=255)                                                                              gCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTgg                                                                  >  1:420399/1‑151 (MQ=255)                                                                              gCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTgg                                                                  <  1:880534/151‑1 (MQ=255)                                                                                aGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCt                                                                <  2:148076/151‑1 (MQ=255)                                                                                 gggTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTg                                                               >  1:80309/1‑151 (MQ=255)                                                                                  ggaaGAAAGCGAAAAACGGCTACGCCCGGAGCGGACCAGCATCCGTTCACCTTCCAGATGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCCCCTGGCTCGTGTGGCt                                                                >  1:482581/4‑149 (MQ=255)                                                                                  ggTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGa                                                              <  1:620653/151‑1 (MQ=255)                                                                                     aGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAAc                                                            <  1:797325/150‑1 (MQ=255)                                                                                        cAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAgc                                                        <  1:795304/151‑1 (MQ=255)                                                                                        cAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAgc                                                        <  1:79268/151‑1 (MQ=255)                                                                                          gCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCTCGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAgcg                                                       <  1:520031/150‑1 (MQ=255)                                                                                              aaaaCAGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAac                                                   >  1:23274/1‑150 (MQ=255)                                                                                                  cGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTg                                              >  2:531068/1‑151 (MQ=255)                                                                                                        cGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGccacc                                         >  2:737007/1‑150 (MQ=255)                                                                                                           ccGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCAttt                                     <  2:556488/151‑1 (MQ=255)                                                                                                           ccGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCAttt                                     >  2:785824/1‑151 (MQ=255)                                                                                                             ggAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTcc                                   <  1:54961/151‑1 (MQ=255)                                                                                                                gCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCAtc                                <  1:616488/151‑1 (MQ=255)                                                                                                                 cGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCAtct                               <  2:753250/151‑1 (MQ=255)                                                                                                                              gTTCACCTTCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTAcc                  >  2:457248/1‑151 (MQ=255)                                                                                                                                      tCCAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAg            <  2:412928/149‑1 (MQ=255)                                                                                                                                        cAGCTGCACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGAc        >  2:357740/1‑151 (MQ=255)                                                                                                                                              cACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGCCAACGTa  <  1:58253/151‑1 (MQ=255)                                                                                                                                              cACGGCACTTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTa  <  2:210871/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                ACGCTTCATCGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTCGGGAAATCCGCCGCTGGCAGGGTAGACAGCGAAAAACGGCTACGCCCGGAGCGTACCAGCATCCGTTCACCTTCCAGCTGCACGGCAATTTCCGCGCCTTCAGGCAGACCACGGCAGATATCAAAGAATTTGCGCGCCGGAACGGTCGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTA  >  NC_000913/3881911‑3882202

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