Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,217	G→T	D174Y (GAC→TAC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,217	0	G	T	100.0%	215.1 / NA	56	D174Y (GAC→TAC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (29/27);  total (29/27)

GGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGA  >  NC_000913/1756069‑1756362                                                                                                                                                     |                                                                                                                                                  ggTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATAc                                                                                                                                                 <  2:697825/151‑1 (MQ=255)  gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACa                                                                                                                                                <  2:321843/151‑1 (MQ=255)  gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACa                                                                                                                                                <  1:412259/151‑1 (MQ=255)   ttGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACa                                                                                                                                                <  2:612274/150‑1 (MQ=255)    tGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACaa                                                                                                                                               >  2:318601/1‑150 (MQ=255)        gATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAgg                                                                                                                                          <  2:678167/151‑1 (MQ=255)        gATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGCCAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACAGGCTGAAAAATACAAACAg                                                                                                                                           >  2:1085415/1‑150 (MQ=255)          tGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGa                                                                                                                                         <  1:438965/150‑1 (MQ=255)            gTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCt                                                                                                                                      <  1:953463/151‑1 (MQ=255)                gATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTg                                                                                                                                     <  2:736983/148‑1 (MQ=255)                gATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGAt                                                                                                                                   <  1:182495/150‑1 (MQ=255)                   cGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCttt                                                                                                                               <  2:525956/151‑1 (MQ=255)                   cGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCttt                                                                                                                               >  1:429865/1‑151 (MQ=255)                     gTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTgg                                                                                                                             >  1:535925/1‑151 (MQ=255)                      tATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGt                                                                                                                            >  1:963559/1‑151 (MQ=255)                       aTGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGtt                                                                                                                           >  2:404781/1‑151 (MQ=255)                       aTGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGtt                                                                                                                           <  1:991916/151‑1 (MQ=255)                         ggAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTTATCTGTGGTTg                                                                                                                          <  1:443995/150‑1 (MQ=255)                          gAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCg                                                                                                                        <  2:117346/151‑1 (MQ=255)                          gAAGTTACCACCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCg                                                                                                                        >  2:631326/1‑151 (MQ=255)                                   gCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCg                                                                                                               >  2:778222/1‑151 (MQ=255)                                      aTTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGt                                                                                                              >  1:579229/1‑149 (MQ=255)                                      aTTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAg                                                                                                            <  1:963656/151‑1 (MQ=255)                                        tGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGAc                                                                                                          >  2:592804/1‑151 (MQ=255)                                         gAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACt                                                                                                         >  1:135131/1‑151 (MQ=255)                                                      ttATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTc                                                                                            >  2:155632/1‑151 (MQ=255)                                                              aaaGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTc                                                                                    >  2:719261/1‑151 (MQ=255)                                                               aaGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCg                                                                                   <  1:648291/151‑1 (MQ=255)                                                                aGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGt                                                                                  <  1:854293/151‑1 (MQ=255)                                                                   gCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTaa                                                                                <  1:208161/150‑1 (MQ=255)                                                                        acaacGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTc                                                                          <  2:584528/151‑1 (MQ=255)                                                                          aacGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTg                                                                        <  1:380542/151‑1 (MQ=255)                                                                                gACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGtt                                                                   <  1:542137/150‑1 (MQ=255)                                                                                 aCCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTAt                                                                 <  1:922536/151‑1 (MQ=255)                                                                                    tGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCg                                                               >  2:221915/1‑150 (MQ=255)                                                                                     gggCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCg                                                               <  1:340427/149‑1 (MQ=255)                                                                                     gggCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGaa                                                             <  2:591825/151‑1 (MQ=255)                                                                                     gggCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGaa                                                             <  1:747862/151‑1 (MQ=255)                                                                                     gggCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGa                                                              >  1:371608/1‑150 (MQ=255)                                                                                         gAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATcc                                                         >  2:251439/1‑151 (MQ=255)                                                                                               aaaGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAg                                                    <  1:318554/150‑1 (MQ=255)                                                                                                 aGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCAc                                                 >  2:710008/1‑151 (MQ=255)                                                                                                  ggTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGCTCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCAcc                                                >  1:1007273/1‑151 (MQ=255)                                                                                                     gtgAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                                              <  2:706859/150‑1 (MQ=255)                                                                                                        aaCCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGaa                                            >  1:255180/1‑149 (MQ=255)                                                                                                        aaCCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGaa                                            >  1:136187/1‑149 (MQ=255)                                                                                                          cctgcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAgcg                                        >  2:568372/1‑151 (MQ=255)                                                                                                            tgcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAgcgc                                       >  1:48205/1‑150 (MQ=255)                                                                                                                tgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCAcggc                                  >  1:696087/1‑151 (MQ=255)                                                                                                                   cGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACggcggc                               <  2:168194/151‑1 (MQ=255)                                                                                                                    gTTTCCATTGCTCTGCCAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGgcggcg                              >  2:455997/1‑151 (MQ=255)                                                                                                                                    cAGCACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCAcatca               >  2:884040/1‑150 (MQ=255)                                                                                                                                       cACTGGCTGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCAtctc           >  1:1036666/1‑151 (MQ=255)                                                                                                                                              tGAAAAATACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATc    >  1:1070073/1‑151 (MQ=255)                                                                                                                                                aaaaaTACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCg   >  1:1105788/1‑150 (MQ=255)                                                                                                                                                 aaaaTACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGa  >  2:249169/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGA  >  NC_000913/1756069‑1756362

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