Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,757,108	T→C	*471Q (TAA→CAA)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,108	0	T	C	100.0%	204.9 / NA	53	*471Q (TAA→CAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/27);  total (26/27)

CTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGA  >  NC_000913/1756974‑1757235                                                                                                                                       |                                                                                                                                cTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAata                                                                                                                             <  1:1497173/139‑1 (MQ=255) cTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAata                                                                                                                             >  1:235835/1‑139 (MQ=255)    atgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAata                                                                                                                             <  1:1114895/136‑1 (MQ=255)    atgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAata                                                                                                                             >  2:1114895/1‑136 (MQ=255)      gatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCt                                                                                                                        >  2:846683/1‑139 (MQ=255)         ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCtttt                                                                                                                     <  2:804349/139‑1 (MQ=255)                 tCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaatta                                                                                                             >  1:1256015/1‑139 (MQ=255)                   tGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           <  2:948127/139‑1 (MQ=255)                   tGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           <  1:157574/139‑1 (MQ=255)                               tGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGaa                                                                                               >  1:761519/1‑139 (MQ=255)                                            cGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATggg                                                                                  >  2:764895/1‑139 (MQ=255)                                             ggTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGgc                                                                                 <  1:1118913/139‑1 (MQ=255)                                                    cacacaAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttt                                                                          <  1:764895/139‑1 (MQ=255)                                                         aaaGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           >  1:535944/1‑101 (MQ=255)                                                         aaaGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           <  2:535944/101‑1 (MQ=255)                                                             gTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                                                                 >  1:1442378/1‑139 (MQ=255)                                                                aCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaa                                                                                                                <  2:931968/89‑1 (MQ=255)                                                                aCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaa                                                                                                                >  1:931968/1‑89 (MQ=255)                                                                aCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              <  2:235835/139‑1 (MQ=255)                                                                           aTGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGc                                                   >  2:1445694/1‑139 (MQ=255)                                                                            tGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattta                                                                   <  2:502498/122‑1 (MQ=255)                                                                            tGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattta                                                                   >  1:502498/1‑122 (MQ=255)                                                                                             gTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           >  1:1410806/1‑65 (MQ=255)                                                                                             gTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           <  2:1410806/65‑1 (MQ=255)                                                                                              tACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                <  2:1442378/139‑1 (MQ=255)                                                                                              tACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                <  2:1256015/139‑1 (MQ=255)                                                                                                 cGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGAttt                             <  2:761519/139‑1 (MQ=255)                                                                                                   aGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAaaatcaaat                                    <  1:36990/130‑1 (MQ=255)                                                                                                   aGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAaaatcaaat                                    >  2:36990/1‑130 (MQ=255)                                                                                                   aGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCa                           <  2:644948/139‑1 (MQ=255)                                                                                                   aGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCa                           >  2:955024/1‑139 (MQ=255)                                                                                                   aGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCa                           >  1:478079/1‑139 (MQ=255)                                                                                                   aGCGGCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCa                           <  1:846683/139‑1 (MQ=255)                                                                                                       gCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           >  2:1158899/1‑55 (MQ=255)                                                                                                       gCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGaattaat                                                                                                           <  1:1158899/55‑1 (MQ=255)                                                                                                       gCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              >  1:804358/1‑100 (MQ=255)                                                                                                       gCACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              <  2:804358/100‑1 (MQ=255)                                                                                                        cACTACTAACACCGCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACtata                      >  2:344669/1‑139 (MQ=255)                                                                                                                     gCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              >  1:1278067/1‑86 (MQ=255)                                                                                                                     gCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              <  2:1278067/86‑1 (MQ=255)                                                                                                                     gCATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTa         >  2:664744/1‑139 (MQ=255)                                                                                                                      cATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              <  2:765867/85‑1 (MQ=255)                                                                                                                      cATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              >  1:765867/1‑85 (MQ=255)                                                                                                                      cATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                <  2:634693/115‑1 (MQ=255)                                                                                                                      cATCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGa                                >  1:634693/1‑115 (MQ=255)                                                                                                                       aTCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattt                                                                    >  1:217056/1‑78 (MQ=255)                                                                                                                       aTCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattt                                                                    <  2:217056/78‑1 (MQ=255)                                                                                                                       aTCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              >  2:748326/1‑84 (MQ=255)                                                                                                                       aTCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                                                              <  1:748326/84‑1 (MQ=255)                                                                                                                        tCTGTTCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCt      <  1:1445694/139‑1 (MQ=255)                                                                                                                            ttCACGTCCTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGa  >  1:184454/1‑139 (MQ=255)                                                                                                                                   ccTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATa                                               >  2:570888/1‑87 (MQ=255)                                                                                                                                   ccTGCAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATa                                               <  1:570888/87‑1 (MQ=255)                                                                                                                                       |                                                                                                                                CTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGA  >  NC_000913/1756974‑1757235

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