Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,718,546	T→C	T353A (ACT→GCT)	anmK ←	anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,718,546	0	T	C	100.0%	144.5 / NA	48	T353A (ACT→GCT)	anmK	anhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/27);  total (21/27)

GAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCA  >  NC_000913/1718401‑1718692                                                                                                                                                  |                                                                                                                                                   gAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACgg                                                                                                                                               >  1:702292/1‑151 (MQ=255)  aaGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACgg                                                                                                                                               >  2:147141/1‑150 (MQ=255)  aaGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGa                                                                                                                                              <  2:141109/151‑1 (MQ=255)         tCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAg                                                                                                                                        <  1:673534/150‑1 (MQ=255)          ccTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGc                                                                                                                                          >  2:242236/1‑147 (MQ=255)                  ggAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCt                                                                                                                                 <  2:4568/148‑1 (MQ=255)                  ggAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGt                                                                                                                              >  1:72541/1‑151 (MQ=255)                   gAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGt                                                                                                                              >  2:32016/1‑150 (MQ=255)                   gAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGt                                                                                                                              >  2:701573/1‑150 (MQ=255)                          aTACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGc                                                                                                                      >  2:376760/1‑151 (MQ=255)                                       aaCGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGGAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGcc                                                                                                           >  2:338902/1‑149 (MQ=255)                                          gTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGc                                                                                                       >  1:196554/1‑150 (MQ=255)                                             ccAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCCCGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCCGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGc                                                                                                   <  2:603823/151‑1 (MQ=255)                                              cAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAACCCCGCCAGCGTCCGCCAGGCAAGcc                                                                                                  <  1:342785/151‑1 (MQ=255)                                               aaGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCa                                                                                                 <  2:190517/151‑1 (MQ=255)                                                   gtAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCg                                                                                             >  1:651973/1‑151 (MQ=255)                                                       cGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGaaa                                                                                          <  1:398851/150‑1 (MQ=255)                                                                   tAAGAAAAATTCAGTTAACTCTGATATCACGTGCTAGCGGTGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGACGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTcc                                                                              >  1:53076/1‑150 (MQ=255)                                                                    aaGAAAAATTCAGTTAACTCTGATATCCCGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCa                                                                             <  1:591496/150‑1 (MQ=255)                                                                        aaaaTTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCa                                                                        <  2:77610/151‑1 (MQ=255)                                                                         aaaTTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAgg                                                                                               >  2:228892/1‑127 (MQ=255)                                                                         aaaTTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAgg                                                                                               <  1:228859/127‑1 (MQ=255)                                                                           aTTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCg                                                                     <  2:513119/151‑1 (MQ=255)                                                                            ttCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGc                                                                    >  2:306342/1‑151 (MQ=255)                                                                                        tGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCa                                                         >  1:641635/1‑150 (MQ=255)                                                                                         gATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCa                                                         <  2:92647/149‑1 (MQ=255)                                                                                         gATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCCCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGcc                                                          <  2:223915/148‑1 (MQ=255)                                                                                          atatCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCGGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCa                                                         <  1:428421/148‑1 (MQ=255)                                                                                           tatCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACgg                                                     >  1:217234/1‑151 (MQ=255)                                                                                            accgcGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCTTGGTTTGCGCCAGCGACGAAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGc                                                    <  2:454235/147‑1 (MQ=255)                                                                                                    ttAGAGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCGGCGACGGAAGGCAGTTTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtg                                             <  2:447998/150‑1 (MQ=255)                                                                                                      aGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggt                                           <  1:267714/150‑1 (MQ=255)                                                                                                        cGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggtgg                                         >  1:626063/1‑150 (MQ=255)                                                                                                          gggAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGa                                      <  2:532157/151‑1 (MQ=255)                                                                                                           ggAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCACTGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGCGGTGAc                                     >  1:636044/1‑151 (MQ=255)                                                                                                               aaTAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCt                                 <  2:269730/151‑1 (MQ=255)                                                                                                                  aGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCtgt                               <  1:137720/150‑1 (MQ=255)                                                                                                                   gCCCCCAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGc                             >  2:238471/1‑151 (MQ=255)                                                                                                                    cccccATTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCtgtg                              <  2:638578/149‑1 (MQ=255)                                                                                                                       ccAGTACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTgg                         <  2:604432/151‑1 (MQ=255)                                                                                                                           tACCGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCa                       <  1:147525/149‑1 (MQ=255)                                                                                                                              cGTCTCCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAAt                  <  2:87841/151‑1 (MQ=255)                                                                                                                                tctcCTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAAt                  <  1:316182/149‑1 (MQ=255)                                                                                                                                   ccTGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgccg              <  2:373347/150‑1 (MQ=255)                                                                                                                                     tGGCTTGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgccgcc            >  2:577987/1‑150 (MQ=255)                                                                                                                                          tGCGCCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGAcgc      >  1:248919/1‑151 (MQ=255)                                                                                                                                           gcgcCAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGAcgcg     >  2:501869/1‑151 (MQ=255)                                                                                                                                              ccAGCGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCa  <  1:748973/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCA  >  NC_000913/1718401‑1718692

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