Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,445	G→T	E250* (GAA→TAA)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,445	0	G	T	100.0%	148.3 / NA	40	E250* (GAA→TAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/23);  total (17/23)

TATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTG  >  NC_000913/1756297‑1756592                                                                                                                                                     |                                                                                                                                                    tATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATaa                                                                                                                                                   <  2:684590/151‑1 (MQ=255)    cGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATc                                                                                                                                                <  2:766291/151‑1 (MQ=255)      aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATcc                                                                                                                                               <  1:449274/150‑1 (MQ=255)         tCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGt                                                                                                                                           >  1:101530/1‑151 (MQ=255)                        aaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCtt                                                                                                                            <  1:444779/151‑1 (MQ=255)                          gcgcACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCg                                                                                                                          >  1:448039/1‑151 (MQ=255)                             cACGGCGGCGAAAACATCCCCATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTCCGccc                                                                                                                       <  1:554898/151‑1 (MQ=255)                             cACCGCGCCGAACACATCCACTTCATCTCCAAATTCGAAAACCAGGAAGGCCTCACTGGCTTCGACGAAAACCTCGAAGCCTCTGACGGTATCATGGTTGCGCGTGGCGACCTGGGTCTATAAATCCCGGTAGAAGAAGATATCTTCGccc                                                                                                                       <  1:555120/151‑1 (MQ=255)                                  cggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCagaag                                                                                                                  <  2:521924/151‑1 (MQ=255)                                       aaaaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgat                                                                                                             <  2:644346/151‑1 (MQ=255)                                           cATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatgat                                                                                                          >  2:375498/1‑150 (MQ=255)                                                 catcatCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaa                                                                                                   >  1:570322/1‑151 (MQ=255)                                                       ctcCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATc                                                                                             <  2:783585/151‑1 (MQ=255)                                                                          gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCg                                                                             <  2:573469/148‑1 (MQ=255)                                                                          gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCg                                                                             <  2:140140/148‑1 (MQ=255)                                                                          gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCg                                                                             <  2:161549/148‑1 (MQ=255)                                                                              gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATc                                                                        <  1:554593/149‑1 (MQ=255)                                                                                             acgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGc                                                        <  2:740093/150‑1 (MQ=255)                                                                                                      tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCAt                                              <  1:499568/151‑1 (MQ=255)                                                                                                      tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCAt                                              <  2:395112/151‑1 (MQ=255)                                                                                                                  aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                   >  2:562277/1‑150 (MQ=255)                                                                                                                  aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                   <  2:281816/150‑1 (MQ=255)                                                                                                                  aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                   <  2:235014/150‑1 (MQ=255)                                                                                                                  aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                   <  1:463318/150‑1 (MQ=255)                                                                                                                  aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATATATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAccc                                  <  2:283504/151‑1 (MQ=255)                                                                                                                   cGGCATCATGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCa                                 >  2:828760/1‑151 (MQ=255)                                                                                                                        tcatGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGcc                             >  2:307209/1‑150 (MQ=255)                                                                                                                           tGGTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGa                          >  2:194137/1‑150 (MQ=255)                                                                                                                             gTTGCGCGTGGCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                        >  2:397774/1‑150 (MQ=255)                                                                                                                                       gCGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAg              >  1:617582/1‑150 (MQ=255)                                                                                                                                        cGACCTGGGTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc             <  2:78680/150‑1 (MQ=255)                                                                                                                                              gggTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgcg                    >  2:90759/1‑137 (MQ=255)                                                                                                                                              gggTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgcg                    <  1:90759/137‑1 (MQ=255)                                                                                                                                               ggTGTATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg     >  2:647721/1‑151 (MQ=255)                                                                                                                                                gtgtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg     <  1:559987/150‑1 (MQ=255)                                                                                                                                                gtgtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt    >  2:820472/1‑151 (MQ=255)                                                                                                                                                  gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  1:565404/1‑151 (MQ=255)                                                                                                                                                  gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  1:697148/1‑151 (MQ=255)                                                                                                                                                  gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  2:60834/1‑151 (MQ=255)                                                                                                                                                  gtATAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg  >  2:163975/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTG  >  NC_000913/1756297‑1756592

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