Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,757,006	G→A	A437T (GCT→ACT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,006	0	G	A	100.0%	126.3 / NA	40	A437T (GCT→ACT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/20);  total (20/20)

CCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756870‑1757130                                                                                                                                         |                                                                                                                             cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGAct                                                                                                                            >  2:264419/1‑139 (MQ=255)                ccTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  2:461251/1‑139 (MQ=255)                          accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  2:233981/1‑139 (MQ=255)                                        aaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGtt                                                                                     <  1:1000303/139‑1 (MQ=255)                                             cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                >  2:799053/1‑139 (MQ=255)                                                     ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTg                                                                        <  1:233981/139‑1 (MQ=255)                                                            tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTGGTAGTTATGGTTTCTGGTGCACTGGt                                                                 >  2:283627/1‑139 (MQ=255)                                                              aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAc                                                               >  1:412549/1‑139 (MQ=255)                                                               gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                              >  2:748352/1‑139 (MQ=255)                                                                    ggCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCg                                                         <  1:283627/139‑1 (MQ=255)                                                                      cGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       <  2:1573089/139‑1 (MQ=255)                                                                      cGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       <  2:412549/139‑1 (MQ=255)                                                                        ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTa                                                                                    <  2:901907/108‑1 (MQ=255)                                                                        ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTa                                                                                    >  1:901907/1‑108 (MQ=255)                                                                                   cTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                  <  1:414120/83‑1 (MQ=255)                                                                                   cTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                  >  2:414120/1‑83 (MQ=255)                                                                                       ttAAAGAGATCACTGCTACTTATGATTTCTACCGTCTGGGTAAAGAACTGAGTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                      <  1:264419/139‑1 (MQ=255)                                                                                             agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaa                                                                                                     >  1:155790/1‑70 (MQ=255)                                                                                             agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaa                                                                                                     <  2:155790/70‑1 (MQ=255)                                                                                              gaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCAc                               <  1:412916/139‑1 (MQ=255)                                                                                               aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                >  2:881475/1‑57 (MQ=255)                                                                                               aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                <  1:881475/57‑1 (MQ=255)                                                                                               aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  2:1591132/1‑110 (MQ=255)                                                                                               aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  1:1591132/110‑1 (MQ=255)                                                                                               aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  2:483283/1‑139 (MQ=255)                                                                                               aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              <  1:748352/139‑1 (MQ=255)                                                                                                 cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  1:103110/1‑58 (MQ=255)                                                                                                 cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             <  2:103110/58‑1 (MQ=255)                                                                                                           gatgatTTCTACCGTCTGGGTAAAGAACTGACTCTg                                                                                                                         >  2:1044411/1‑36 (MQ=255)                                                                                                           gatgatTTCTACCGTCTGGGTAAAGAACTGACTCTg                                                                                                                         <  1:1044411/36‑1 (MQ=255)                                                                                                              gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  2:1465272/80‑1 (MQ=255)                                                                                                              gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  1:1465272/1‑80 (MQ=255)                                                                                                              gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGa                                                            <  1:809113/94‑1 (MQ=255)                                                                                                              gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGa                                                            >  2:809113/1‑94 (MQ=255)                                                                                                                 ttCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             <  2:357451/42‑1 (MQ=255)                                                                                                                 ttCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  1:357451/1‑42 (MQ=255)                                                                                                                      ccGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            >  2:1381429/1‑54 (MQ=255)                                                                                                                      ccGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            <  1:1381429/54‑1 (MQ=255)                                                                                                                           tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  >  1:983499/1‑139 (MQ=255)                                                                                                                           tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  <  2:556486/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                             CCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756870‑1757130

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