Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,434	A→G	D246G (GAC→GGC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,434	0	A	G	100.0%	208.9 / NA	54	D246G (GAC→GGC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/29);  total (25/29)

TCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAG  >  NC_000913/1756289‑1756580                                                                                                                                                  |                                                                                                                                                   tCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTg                                                                                                                                                <  1:999483/150‑1 (MQ=255)    gACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGgtg                                                                                                                                            <  2:333323/151‑1 (MQ=255)     aCGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGgtgt                                                                                                                                           >  2:212463/1‑151 (MQ=255)          aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGaaa                                                                                                                                      <  1:122786/151‑1 (MQ=255)           tCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAAt                                                                                                                                     >  2:568043/1‑151 (MQ=255)              aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATcct                                                                                                                                  >  2:10147/1‑150 (MQ=255)              aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATccc                                                                                                                                  <  1:201710/151‑1 (MQ=255)               aaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCg                                                                                                                                 >  1:488948/1‑151 (MQ=255)                aTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCgg                                                                                                                                >  1:37541/1‑151 (MQ=255)                aTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCgg                                                                                                                                >  1:861656/1‑151 (MQ=255)                  ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGt                                                                                                                               <  2:242378/150‑1 (MQ=255)                  ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTa                                                                                                                              <  2:791019/151‑1 (MQ=255)                  ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACACCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTa                                                                                                                              <  1:564742/151‑1 (MQ=255)                   cGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTag                                                                                                                             <  1:508771/151‑1 (MQ=255)                    gTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGt                                                                                                                               >  2:1034698/1‑148 (MQ=255)                         cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTagaaga                                                                                                                         >  2:735416/1‑149 (MQ=255)                         cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGaagaag                                                                                                                       <  1:506441/151‑1 (MQ=255)                              gAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATc                                                                                                                  <  2:551414/151‑1 (MQ=255)                               aaaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCt                                                                                                                 >  1:567030/1‑151 (MQ=255)                               aaaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCt                                                                                                                 >  1:729878/1‑151 (MQ=255)                                aaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCtt                                                                                                                >  1:411420/1‑151 (MQ=255)                                    gcACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGcc                                                                                                            <  2:513251/151‑1 (MQ=255)                                          cggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCagaag                                                                                                      <  2:644070/151‑1 (MQ=255)                                           ggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCaga                                                                                                        >  1:822988/1‑148 (MQ=255)                                             cgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatg                                                                                                   <  1:526697/151‑1 (MQ=255)                                                 aaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatga                                                                                               <  2:604049/151‑1 (MQ=255)                                                   cATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATc                                                                                             >  2:830856/1‑151 (MQ=255)                                                         catcatCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaa                                                                                       <  2:806305/151‑1 (MQ=255)                                                          atcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt                                                                                      <  2:706217/151‑1 (MQ=255)                                                          atcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGAGGATGATCGAAAAAt                                                                                      >  2:253550/1‑151 (MQ=255)                                                                   aaaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTg                                                                             >  1:48927/1‑151 (MQ=255)                                                                      aTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCAc                                                                          <  1:895916/151‑1 (MQ=255)                                                                      aTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCAc                                                                          <  1:762284/151‑1 (MQ=255)                                                                        cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGt                                                                        >  2:703744/1‑151 (MQ=255)                                                                             aCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGt                                                                   <  1:170865/151‑1 (MQ=255)                                                                               cAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCg                                                                 <  2:397991/151‑1 (MQ=255)                                                                                 ggAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGtt                                                               <  1:984404/151‑1 (MQ=255)                                                                                       ccTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACt                                                         <  2:727767/151‑1 (MQ=255)                                                                                        cTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                        >  2:542389/1‑151 (MQ=255)                                                                                         tCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGc                                                       >  1:169706/1‑151 (MQ=255)                                                                                          caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGc                                                       <  1:628504/150‑1 (MQ=255)                                                                                                      cgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTg                                          <  1:877273/151‑1 (MQ=255)                                                                                                        aaaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTgg                                         >  2:673631/1‑150 (MQ=255)                                                                                                           tCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTc                                     >  2:421115/1‑151 (MQ=255)                                                                                                             cTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTcc                                    >  2:635249/1‑150 (MQ=255)                                                                                                                            ggCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCAc                    <  1:173695/151‑1 (MQ=255)                                                                                                                              catcatGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGc                  >  2:664237/1‑151 (MQ=255)                                                                                                                                 catGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGGCCCAGACGCTGGATTCCATGATCAAAAACCCACGCCCg               <  1:364371/151‑1 (MQ=255)                                                                                                                                  atGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg               >  2:850423/1‑150 (MQ=255)                                                                                                                                    ggTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc             <  2:300064/150‑1 (MQ=255)                                                                                                                                    ggTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc             >  2:13880/1‑150 (MQ=255)                                                                                                                                    ggTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt            >  1:340867/1‑151 (MQ=255)                                                                                                                                     gTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTc           <  1:212427/151‑1 (MQ=255)                                                                                                                                            gTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGa    <  2:920376/151‑1 (MQ=255)                                                                                                                                              ggcggcCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGaa   <  1:146325/150‑1 (MQ=255)                                                                                                                                              ggcggcCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAg  <  2:5566/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAG  >  NC_000913/1756289‑1756580

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