Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,152	T→C	L152P (CTG→CCG)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,152	0	T	C	100.0%	222.7 / NA	63	L152P (CTG→CCG)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/32);  total (31/32)

CGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAG  >  NC_000913/1756023‑1756285                                                                                                                                  |                                                                                                                                      cGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGaaaa                                                                                                                              >  2:1312585/1‑139 (MQ=255)   tATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACa                                                                                                                            <  2:1128881/139‑1 (MQ=255)    aTGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACaa                                                                                                                           <  2:564268/139‑1 (MQ=255)                       ccTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGc                                                                                                        <  2:1302590/139‑1 (MQ=255)                                tGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCAtt                                                                                               >  1:934280/1‑139 (MQ=255)                                   cAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGct                                                                                            <  1:1312585/139‑1 (MQ=255)                                     acacCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGct                                                                                            >  2:517779/1‑137 (MQ=255)                                     acacCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGct                                                                                            <  1:517779/137‑1 (MQ=255)                                                gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACa                                                                                                                            <  2:1864175/94‑1 (MQ=255)                                                gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACa                                                                                                                            >  1:1864175/1‑94 (MQ=255)                                                gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGct                                                                                            >  2:269099/1‑126 (MQ=255)                                                gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGct                                                                                            <  1:269099/126‑1 (MQ=255)                                                   gACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTg                                                                            >  1:953863/1‑139 (MQ=255)                                                   gACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTg                                                                            >  1:1729229/1‑139 (MQ=255)                                                          gTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGACTGGCGTTTCCATTGCTCTGCCAGCACTGGACGaaagaca                                                                     >  1:479601/1‑135 (MQ=255)                                                           tCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAg                                                                      >  2:788510/1‑137 (MQ=255)                                                           tCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAg                                                                      <  1:788510/137‑1 (MQ=255)                                                               aTCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGATCAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAAc                                                                <  2:934280/139‑1 (MQ=255)                                                                    tATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAc                                                           >  1:338293/1‑139 (MQ=255)                                                                    tATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAc                                                           <  1:1537081/139‑1 (MQ=255)                                                                    tATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAc                                                           <  1:728812/139‑1 (MQ=255)                                                                     aTGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAcc                                                          <  2:338293/139‑1 (MQ=255)                                                                     aTGGAAGTGACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTCTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCATTGGCTGAAAAAGACAAACAGGAcc                                                          <  1:1376532/139‑1 (MQ=255)                                                                           gTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCt                                                    >  2:913081/1‑139 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGaaaa                                                                                                                              <  2:727572/62‑1 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGaaaa                                                                                                                              >  1:727572/1‑62 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGgtg                                                                                                                      <  2:1646432/70‑1 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGgtg                                                                                                                      >  1:1646432/1‑70 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAcc                                                          <  2:1622748/130‑1 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAcc                                                          >  1:1622748/1‑130 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTg                                                 >  2:1638633/1‑139 (MQ=255)                                                                              aCCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTg                                                 >  2:938762/1‑139 (MQ=255)                                                                                 gCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGtt                                              <  1:573843/139‑1 (MQ=255)                                                                                       gAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAGACAAAGGTGTGAACCTGCCTGGCCTTTCCATTGCTCTGCCAGCACTGGCTAAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAAc                                        <  2:1195218/139‑1 (MQ=255)                                                                                                  aGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGt                                                                                                      <  2:533298/66‑1 (MQ=255)                                                                                                  aGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGt                                                                                                      >  1:533298/1‑66 (MQ=255)                                                                                                  aGTTATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGAc                             <  2:953863/139‑1 (MQ=255)                                                                                                     tATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGaaaa                                                                                                                              >  2:775788/1‑39 (MQ=255)                                                                                                     tATCTGTAAAGTGCTGAACAACGGTGACCCGGGCGaaaa                                                                                                                              <  1:775788/39‑1 (MQ=255)                                                                                                      aTCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACaa                                                                                                                           <  1:11385/41‑1 (MQ=255)                                                                                                      aTCTGTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACaa                                                                                                                           >  2:11385/1‑41 (MQ=255)                                                                                                          gTAAAGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGCTTGCGAACAAGGCGTAGACTTTGTtgc                     >  2:127483/1‑139 (MQ=255)                                                                                                             aaGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTtgctgc                  <  1:1638633/139‑1 (MQ=255)                                                                                                              aGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGt                                                                                                      <  1:82738/54‑1 (MQ=255)                                                                                                              aGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGt                                                                                                      >  2:82738/1‑54 (MQ=255)                                                                                                              aGTGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTgctgct                 >  2:2732/1‑139 (MQ=255)                                                                                                                tGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAcc                                                          >  2:501977/1‑96 (MQ=255)                                                                                                                tGCTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGAcc                                                          <  1:501977/96‑1 (MQ=255)                                                                                                                  cTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCa                                                                                   <  2:231331/69‑1 (MQ=255)                                                                                                                  cTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCa                                                                                   >  1:231331/1‑69 (MQ=255)                                                                                                                  cTGAACAACGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCt             >  2:716334/1‑139 (MQ=255)                                                                                                                     aacaacGGTGACCCGGGCGAAAACAAAGGTGTGaa                                                                                                                  >  2:1798558/1‑35 (MQ=255)                                                                                                                     aacaacGGTGACCCGGGCGAAAACAAAGGTGTGaa                                                                                                                  <  1:1798558/35‑1 (MQ=255)                                                                                                                     aacaacGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCa                                                                                   <  2:74723/66‑1 (MQ=255)                                                                                                                     aacaacGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCa                                                                                   >  1:74723/1‑66 (MQ=255)                                                                                                                     aacaacGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAAc                                        >  1:1824012/1‑109 (MQ=255)                                                                                                                     aacaacGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAAc                                        <  2:1824012/109‑1 (MQ=255)                                                                                                                          cGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCtt                                                   <  2:696302/93‑1 (MQ=255)                                                                                                                          cGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCtt                                                   >  1:696302/1‑93 (MQ=255)                                                                                                                          cGGTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGt     >  1:664603/1‑139 (MQ=255)                                                                                                                            gTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAAc                                        <  1:916449/102‑1 (MQ=255)                                                                                                                            gTGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAAc                                        >  2:916449/1‑102 (MQ=255)                                                                                                                             tGACCCGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAg  <  2:1462220/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      CGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAG  >  NC_000913/1756023‑1756285

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