Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,495	(A)5→4	coding (798/1413 nt)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,491	0	A	.	100.0%	147.9 / NA	35	coding (794/1413 nt)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/16);  total (19/16)

AAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGA  >  NC_000913/1756356‑1756617                                                                                                                                        |                                                                                                                               aaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaa                                                                                                                              >  1:51072/1‑138 (MQ=255) aaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaa                                                                                                                              <  2:51072/138‑1 (MQ=255)     cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaatgta                                                                                                                         >  1:1093231/1‑135 (MQ=255)     cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaatgta                                                                                                                         >  1:889574/1‑135 (MQ=255)                   gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTaaa                                                                                                          <  2:923347/139‑1 (MQ=255)                     cTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGt                                                                                                        <  2:323211/139‑1 (MQ=255)                             cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCa                                                                                                >  1:368037/1‑139 (MQ=255)                             cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCa                                                                                                >  2:524239/1‑139 (MQ=255)                              ttCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCAc                                                                                               >  1:909116/1‑139 (MQ=255)                                cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                             <  2:368037/139‑1 (MQ=255)                                cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                             <  2:117217/139‑1 (MQ=255)                                            cGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGc                                                                                 <  1:761066/139‑1 (MQ=255)                                                           catcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCa                                                                  >  2:805680/1‑139 (MQ=255)                                                             tcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaaa                                                                >  2:714582/1‑139 (MQ=255)                                                                 ggTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                                            >  1:168202/1‑139 (MQ=255)                                                                 ggTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                                            <  2:193658/139‑1 (MQ=255)                                                                      cgcgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                       >  1:79894/1‑139 (MQ=255)                                                                                     gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTg                                                                                                                  <  1:537609/65‑1 (MQ=255)                                                                                     gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTg                                                                                                                  >  2:537609/1‑65 (MQ=255)                                                                                       gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc                                      <  2:1093231/139‑1 (MQ=255)                                                                                          gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                   <  2:339577/139‑1 (MQ=255)                                                                                              tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                               >  1:364160/1‑139 (MQ=255)                                                                                              tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                               >  2:648790/1‑139 (MQ=255)                                                                                                ccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCGCTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                             >  1:636795/1‑139 (MQ=255)                                                                                                  ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                             >  1:350490/1‑73 (MQ=255)                                                                                                  ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                             <  2:350490/73‑1 (MQ=255)                                                                                                  ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                           >  2:1020096/1‑139 (MQ=255)                                                                                                   gTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                          >  2:329803/1‑139 (MQ=255)                                                                                                         gaagTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTAt                                                                                                  >  1:565465/1‑61 (MQ=255)                                                                                                         gaagTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTAt                                                                                                  <  2:565465/61‑1 (MQ=255)                                                                                                          aagTTATCTTCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                   <  2:528491/139‑1 (MQ=255)                                                                                                                   tCGCCCAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc          <  2:990997/139‑1 (MQ=255)                                                                                                                      cccAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt       >  2:534451/1‑139 (MQ=255)                                                                                                                      cccAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt       >  1:580962/1‑139 (MQ=255)                                                                                                                        cAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTAc     <  2:909116/139‑1 (MQ=255)                                                                                                                        cAGAAGATGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTAc     <  1:302804/139‑1 (MQ=255)                                                                                                                          gaagaTGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg   >  1:665183/1‑139 (MQ=255)                                                                                                                           aagaTGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGa  <  2:636795/139‑1 (MQ=255)                                                                                                                           aagaTGATGATCG‑AAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGAGACGTTGCAAACGCCATCCTCGACGGTACTGa  >  2:828580/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                               AAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGA  >  NC_000913/1756356‑1756617

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