Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,910	(A)6→5	coding (1213/1413 nt)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,905	0	A	.	100.0%	254.4 / NA	60	coding (1208/1413 nt)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (30/30);  total (30/30)

GGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGT  >  NC_000913/1756778‑1757035                                                                                                                                |                                                                                                                                   ggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCtca                                                                                                                        >  2:630154/1‑139 (MQ=255)       gTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGtt                                                                                                                  >  1:604613/1‑139 (MQ=255)          gAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAGACGGCTCATCAGTTGGt                                                                                                               <  1:56408/139‑1 (MQ=255)           aaaCTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTa                                                                                                              <  2:903684/139‑1 (MQ=255)                  tgAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGaaaaacgg                                                                                                                             <  1:1080007/118‑4 (MQ=255)                  tgAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGaaaaacgg                                                                                                                             >  2:1080007/1‑115 (MQ=255)                           ggATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTtg                                                                                              >  1:507628/1‑139 (MQ=255)                                  ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  2:604613/139‑1 (MQ=255)                                  ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  2:218273/139‑1 (MQ=255)                                  ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  1:630154/139‑1 (MQ=255)                                          cGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑ATAAACGGCTCATCAGt                                                                                                                   <  2:739241/103‑1 (MQ=255)                                          cGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCAGCAGt                                                                                                                   >  1:739241/1‑103 (MQ=255)                                               ttGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                          <  1:664018/139‑1 (MQ=255)                                               ttGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                          <  1:996649/139‑1 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGaaaaacgg                                                                                                                             >  1:964801/1‑83 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGaaaaacgg                                                                                                                             <  2:964801/86‑4 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       >  1:1005191/1‑123 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  1:124992/123‑1 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  2:1005191/123‑1 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       >  2:124992/1‑123 (MQ=255)                                                  cTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATc                                                                       >  1:508058/1‑139 (MQ=255)                                                    aCTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg                                                                                  >  1:304996/1‑126 (MQ=255)                                                    aCTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg                                                                                  <  2:304996/126‑1 (MQ=255)                                                        aGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTa                                                                                                              >  2:27940/1‑94 (MQ=255)                                                        aGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTa                                                                                                              <  1:27940/94‑1 (MQ=255)                                                          ggCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTAc                                                               >  1:825559/1‑139 (MQ=255)                                                            cGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtg                                                             <  2:761755/139‑1 (MQ=255)                                                             ggTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTa                                                                <  2:483688/135‑1 (MQ=255)                                                             ggTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTa                                                                >  1:483688/1‑135 (MQ=255)                                                               tAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCg                                                                                         <  1:219571/108‑1 (MQ=255)                                                               tAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCg                                                                                         >  2:219571/1‑108 (MQ=255)                                                                                 aCGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtg                                                             >  2:403417/1‑118 (MQ=255)                                                                                 aCGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtg                                                             <  1:403417/118‑1 (MQ=255)                                                                                        tACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  1:793617/85‑1 (MQ=255)                                                                                        tACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       >  2:793617/1‑85 (MQ=255)                                                                                          cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTgg                               >  1:229978/1‑139 (MQ=255)                                                                                          cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCCACTGATGATTTCTACCGTCTGGGTAAAGAACTgg                               >  1:463271/1‑139 (MQ=255)                                                                                            tCCCGGATGCCACCATCCTGTCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGct                             >  2:650087/1‑139 (MQ=255)                                                                                            tCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAAc                                  >  2:967996/1‑134 (MQ=255)                                                                                            tCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAAc                                  <  1:967996/134‑1 (MQ=255)                                                                                            tCCCGGATGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGct                             >  1:1104702/1‑139 (MQ=255)                                                                                                   tGCCACCATCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCaga                      <  2:833233/139‑1 (MQ=255)                                                                                                        ccaTCCTGGCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt                 <  1:182231/139‑1 (MQ=255)                                                                                                            ccTGGCACTGACCACCAACG‑ACAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTgg             <  1:650087/139‑1 (MQ=255)                                                                                                                gCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt                 >  2:319129/1‑131 (MQ=255)                                                                                                                gCACTGACCACCAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt                 <  1:319129/131‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       >  1:750757/1‑55 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                       <  2:750757/55‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGc                                                                                     >  1:983468/1‑57 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGc                                                                                     <  2:983468/57‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                          >  2:623930/1‑68 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                          <  1:623930/68‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAt                                                                        >  2:563707/1‑70 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAt                                                                        <  1:563707/70‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                                     >  2:677793/1‑73 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                                     <  1:677793/73‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtg                                                             >  1:384068/1‑81 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtg                                                             <  2:384068/81‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGt                                               <  1:559380/95‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGt                                               >  2:559380/1‑95 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGt                                               >  1:358759/1‑95 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGt                                               <  2:358759/95‑1 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg   >  2:146128/1‑139 (MQ=255)                                                                                                                      accaccAACG‑AAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCT‑GGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGt  >  1:902646/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                   GGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGT  >  NC_000913/1756778‑1757035

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