Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,529	A→T	T278S (ACC→TCC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,529	0	A	T	98.6%	242.5 / NA	70	T278S (ACC→TCC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (38/31);  minor base G (1/0);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGA  >  NC_000913/1756431‑1756625                                                                                                   |                                                                                                 gCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGag                                                                                                 >  1:1049129/1‑99 (MQ=255)    aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCa                                                                                              <  1:763107/100‑1 (MQ=255)    aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCa                                                                                              <  1:1786554/100‑1 (MQ=255)      cTGGGTGTAGAAATCCCGGTAGAAGAATTTATCTTTGCCCAGAAGTTGATGATCGAATACTGTATCGGTGCACGTTATGTTTTTATCTCTTCGTCCCaa                                                                                             <  2:933531/99‑2 (MQ=255)       tGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGAt                                                                                           >  2:280501/1‑100 (MQ=255)          gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCt                                                                                        >  1:1269863/1‑100 (MQ=255)           tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                       <  2:1662152/100‑1 (MQ=255)            gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTgg                                                                                      <  1:1349095/100‑1 (MQ=255)            gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTgg                                                                                      <  1:1657975/100‑1 (MQ=255)             tAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGa                                                                                     >  2:1086027/1‑100 (MQ=255)              aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                    >  2:100880/1‑100 (MQ=255)               gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTGGTTATCACTGCGTCGCAGATGCTGGAtt                                                                                   >  2:1104636/1‑100 (MQ=255)               gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAtt                                                                                   <  1:1720507/100‑1 (MQ=255)                    cccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATg                                                                              >  2:1225562/1‑100 (MQ=255)                      cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGAt                                                                            <  2:1269860/100‑1 (MQ=255)                      cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGAt                                                                            <  1:1086029/100‑1 (MQ=255)                       ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATc                                                                           <  2:1710546/100‑1 (MQ=255)                        gTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTGAAATCGTTAACACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                          >  2:1617980/1‑100 (MQ=255)                            aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCaaaaa                                                                      <  1:100879/100‑1 (MQ=255)                            aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGGCCCAGATGCTGGATTCCATGATCaaaaa                                                                      >  2:1376878/1‑100 (MQ=255)                              gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAcc                                                                    <  1:238823/100‑1 (MQ=255)                              gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAcc                                                                    >  1:342240/1‑100 (MQ=255)                               nagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAccc                                                                   >  2:1083247/2‑100 (MQ=255)                               aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAccc                                                                   >  1:911934/1‑100 (MQ=255)                                agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCa                                                                  <  2:857040/100‑1 (MQ=255)                                agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCa                                                                  >  2:1225500/1‑100 (MQ=255)                                   tATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                               <  1:1332775/100‑1 (MQ=255)                                     tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGCTTCCATGATCAAAAACCCACGcac                                                             >  2:180395/1‑98 (MQ=255)                                     tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGccc                                                             >  2:1773797/1‑100 (MQ=255)                                      cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                            >  2:2269554/1‑100 (MQ=255)                                      cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                            >  2:1327080/1‑100 (MQ=255)                                      cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                            >  2:1166939/1‑100 (MQ=255)                                        tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc                                                          <  1:1773802/100‑1 (MQ=255)                                         cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                         <  2:2177980/100‑1 (MQ=255)                                          gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTc                                                        <  1:1617985/100‑1 (MQ=255)                                            ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgc                                                      >  1:513004/1‑100 (MQ=255)                                              agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCgcgc                                                    <  2:1491963/100‑1 (MQ=255)                                              agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCgcgc                                                    >  2:610558/1‑100 (MQ=255)                                                aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAg                                                  <  1:1144444/100‑1 (MQ=255)                                                aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAg                                                  <  2:891288/100‑1 (MQ=255)                                                aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAg                                                  <  1:1376881/100‑1 (MQ=255)                                                 agaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGa                                                 <  1:1225503/100‑1 (MQ=255)                                                    tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc                                              >  2:641572/1‑100 (MQ=255)                                                     gatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                             <  1:641573/100‑1 (MQ=255)                                                       cgatCTAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                           <  1:461095/99‑1 (MQ=255)                                                         atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                          >  2:2178824/1‑99 (MQ=255)                                                         atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                         >  2:852764/1‑100 (MQ=255)                                                            gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                      >  1:1523878/1‑100 (MQ=255)                                                             aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGgagatcg                                     >  2:461097/1‑97 (MQ=255)                                                             aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                     <  1:852766/100‑1 (MQ=255)                                                                aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                  >  2:1814451/1‑100 (MQ=255)                                                                    tATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                              >  1:809735/1‑100 (MQ=255)                                                                     ccgatctcACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                             <  1:614739/93‑1 (MQ=255)                                                                        cGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                          >  2:2048502/1‑100 (MQ=255)                                                                        cGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                          <  2:809733/100‑1 (MQ=255)                                                                            cACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACgagatcgg                      >  2:614739/1‑93 (MQ=255)                                                                            cACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                      <  2:1688478/100‑1 (MQ=255)                                                                              cGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                    >  1:635658/1‑100 (MQ=255)                                                                               gTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAAGTTGCAAACGCAATCCTCGa                   >  2:2009205/1‑100 (MQ=255)                                                                                tAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                  >  1:731592/1‑100 (MQ=255)                                                                                tAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGACCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                  <  2:911932/100‑1 (MQ=255)                                                                                   aGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt               >  1:1184602/1‑100 (MQ=255)                                                                                    gTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTa              >  1:2224926/1‑100 (MQ=255)                                                                                      cGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt            >  2:691555/1‑100 (MQ=255)                                                                                          aTCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg        <  2:2224916/100‑1 (MQ=255)                                                                                            cACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa      >  1:1638536/1‑100 (MQ=255)                                                                                            cACTGCATCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa      >  2:1101541/1‑100 (MQ=255)                                                                                             aCTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg     >  1:920254/1‑100 (MQ=255)                                                                                             aCTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg     <  1:1594896/100‑1 (MQ=255)                                                                                              cTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt    >  2:1373604/1‑100 (MQ=255)                                                                                                gCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGa  >  1:2045563/1‑100 (MQ=255)                                                                                                   |                                                                                                 GCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGA  >  NC_000913/1756431‑1756625

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