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%	193.1 / NA	60	A437T (GCT→ACT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/31);  total (29/31)

CGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCG  >  NC_000913/1756872‑1757140                                                                                                                                       |                                                                                                                                       cGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGActct                                                                                                                                    >  1:1109387/1‑139 (MQ=255)          ccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCgg                                                                                                                           <  1:1056449/139‑1 (MQ=255)           caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                          <  2:228359/139‑1 (MQ=255)                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                              <  2:739827/122‑1 (MQ=255)                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                              >  1:739827/1‑122 (MQ=255)                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                             >  2:1240508/1‑139 (MQ=255)                            ccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                                         <  1:1240508/139‑1 (MQ=255)                                  aaaaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGt                                                                                                   >  1:1145843/1‑139 (MQ=255)                                           cTCATCAGTTGGTACTGAGCAAAGGGGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAGAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                          <  2:851714/139‑1 (MQ=255)                                                        aCTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAg                                                                                                 <  1:640477/119‑1 (MQ=255)                                                        aCTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAg                                                                                                 >  2:640477/1‑119 (MQ=255)                                                           gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                                    <  1:348227/129‑1 (MQ=255)                                                           gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                                    >  2:348227/1‑129 (MQ=255)                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                            >  2:261482/1‑104 (MQ=255)                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                            <  1:261482/104‑1 (MQ=255)                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTg                                                                                                          <  2:929037/106‑1 (MQ=255)                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTg                                                                                                          >  1:929037/1‑106 (MQ=255)                                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcac                                                                                                                   <  1:1183150/96‑1 (MQ=255)                                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcac                                                                                                                   >  2:1183150/1‑96 (MQ=255)                                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                                >  2:784751/1‑131 (MQ=255)                                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                                <  1:784751/131‑1 (MQ=255)                                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                                        >  1:37166/1‑139 (MQ=255)                                                                 aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                              >  2:864721/1‑81 (MQ=255)                                                                 aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                              <  1:864721/81‑1 (MQ=255)                                                                 aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                                     >  2:730160/1‑138 (MQ=255)                                                                 aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                                     <  1:730160/138‑1 (MQ=255)                                                                             gCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCATTACTAac                                                        >  1:329389/1‑139 (MQ=255)                                                                              cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                          >  2:146718/1‑104 (MQ=255)                                                                              cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                          <  1:146718/104‑1 (MQ=255)                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                                 >  1:184223/1‑114 (MQ=255)                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                                 <  2:184223/114‑1 (MQ=255)                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                        <  1:1153366/139‑1 (MQ=255)                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                        >  2:145747/1‑139 (MQ=255)                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                        <  1:145747/139‑1 (MQ=255)                                                                                              tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGt                                       >  1:1318359/1‑139 (MQ=255)                                                                                               cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTc                                      <  2:528869/139‑1 (MQ=255)                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                          <  1:156066/44‑1 (MQ=255)                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                          >  2:156066/1‑44 (MQ=255)                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                                      >  1:608213/1‑64 (MQ=255)                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                                      <  2:608213/64‑1 (MQ=255)                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTg                           >  1:649709/1‑139 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTc                                                                                                                         <  1:1024626/43‑1 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTc                                                                                                                         >  2:1024626/1‑43 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                            >  1:1042815/1‑56 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                            <  1:1158193/56‑1 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                            >  2:1158193/1‑56 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                            <  2:1042815/56‑1 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCactac                                                            >  2:614665/1‑104 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCactac                                                            <  1:614665/104‑1 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                           >  2:847671/1‑121 (MQ=255)                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                           <  1:847671/121‑1 (MQ=255)                                                                                                                    ccGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaa                 <  2:649709/139‑1 (MQ=255)                                                                                                                      gTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCAc                                                                               >  2:44436/1‑75 (MQ=255)                                                                                                                      gTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCAc                                                                               <  1:44436/75‑1 (MQ=255)                                                                                                                       tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaatta              >  2:465249/1‑139 (MQ=255)                                                                                                                         tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat            <  1:465249/139‑1 (MQ=255)                                                                                                                           ggTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAttt          <  2:1109387/139‑1 (MQ=255)                                                                                                                                aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGtt                                            <  2:24559/100‑1 (MQ=255)                                                                                                                                aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGtt                                            >  1:24559/1‑100 (MQ=255)                                                                                                                                   aCTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCg  <  2:329389/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       CGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCG  >  NC_000913/1756872‑1757140

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