Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,755,770	G→A	D25N (GAC→AAC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,770	0	G	A	100.0%	153.1 / NA	42	D25N (GAC→AAC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/18);  total (24/18)

GGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGG  >  NC_000913/1755625‑1755908                                                                                                                                                  |                                                                                                                                           ggATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCt                                                                                                                                       <  1:942200/151‑1 (MQ=255) ggATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCt                                                                                                                                       >  2:721084/1‑151 (MQ=255)   aTTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTgg                                                                                                                                     <  1:222976/151‑1 (MQ=255)         tttCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGa                                                                                                                                >  2:601929/1‑150 (MQ=255)                  gTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCg                                                                                                                      <  1:360589/151‑1 (MQ=255)                    gcgcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGt                                                                                                                     <  1:878406/150‑1 (MQ=255)                      gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTg                                                                                                                  <  2:863051/151‑1 (MQ=255)                      gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTg                                                                                                                  >  2:558111/1‑151 (MQ=255)                       cccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTg                                                                                                                  >  1:186402/1‑150 (MQ=255)                        ccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGa                                                                                                                 >  2:860040/1‑150 (MQ=255)                        ccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGa                                                                                                                 >  1:577102/1‑150 (MQ=255)                          aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACt                                                                                                              >  1:289880/1‑151 (MQ=255)                             aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTtct                                                                                                           >  1:525178/1‑151 (MQ=255)                                      tctcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTg                                                                                                  >  1:557080/1‑151 (MQ=255)                                          ccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACt                                                                                               >  1:491266/1‑150 (MQ=255)                                           cATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTAt                                                                                             <  2:620766/151‑1 (MQ=255)                                                          aaaGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAg                                                                              <  2:338691/151‑1 (MQ=255)                                                            agactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAgcg                                                                            >  2:593475/1‑151 (MQ=255)                                                             gactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAgcgc                                                                           >  2:146976/1‑151 (MQ=255)                                                               ctaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCAt                                                                         >  2:21548/1‑151 (MQ=255)                                                                 aagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCAtt                                                                        <  2:787228/150‑1 (MQ=255)                                                                 aagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCAt                                                                         <  1:509327/149‑1 (MQ=255)                                                                    actGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGa                                                                    <  1:34173/151‑1 (MQ=255)                                                                            gAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTgcgc                                                            <  2:285567/151‑1 (MQ=255)                                                                             aaaaaGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTgcgc                                                            >  1:516234/1‑150 (MQ=255)                                                                              aaaaGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCa                                                           >  2:835566/1‑150 (MQ=255)                                                                                   aCCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtg                                                      >  2:382154/1‑150 (MQ=255)                                                                                           tGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCaa                                              <  2:498257/150‑1 (MQ=255)                                                                                            gTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCaaa                                             >  2:424842/1‑150 (MQ=255)                                                                                            gTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCaaa                                             >  2:783133/1‑150 (MQ=255)                                                                                              ttGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACt                                          >  2:787750/1‑151 (MQ=255)                                                                                                 cACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGt                                       <  1:360109/151‑1 (MQ=255)                                                                                                   ccATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaa                                     <  2:899171/151‑1 (MQ=255)                                                                                                    cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaa                                    <  1:820708/151‑1 (MQ=255)                                                                                                    cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaa                                     >  1:825446/1‑150 (MQ=255)                                                                                                                      aTCTGAAGAGATGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCtt                  >  1:28695/1‑151 (MQ=255)                                                                                                                            agagATGTTAGATAAAATGGTGAACGCTGGAATAAACTTTATGCGTCTGATCTTCTCTCTGGGTAAATAGAAAAAACGCGGTCGGCGAATGCAGATTCTGCGCATGGTAAGGAAAAAATTTGGAAAAACCGCCGCTATCCTGCTTGATa              <  1:460322/149‑1 (MQ=255)                                                                                                                             gagaTGTTAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATAcc            >  2:915259/1‑150 (MQ=255)                                                                                                                                 tgttaaaaaaaaTGCTGAACGCTGGCATGACCTTTATGCGTCTAAACTTCTCTCATGGTGACTAGGCAGAACCCGGTCAGCGCTTTCAAAATCTGCGCACCGTGATGAGAAAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCaa          <  2:307031/140‑1 (MQ=255)                                                                                                                                    tAGCTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTcc    >  1:591210/1‑151 (MQ=255)                                                                                                                                     agaaaaaaTGCTGAACGCTGGAATGAACTTTATGCGTCTGACTTTCTCTCATGGTGACTAGGAAGAACACGGTCAGCGAATTAAGAATCTGCACAACGTGAGGACAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCg   <  2:820903/147‑1 (MQ=255)                                                                                                                                       cTAAAATGCTGAACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCgg  >  1:33510/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                           GGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGG  >  NC_000913/1755625‑1755908

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