Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,709	G→A	V338M (GTG→ATG)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,709	0	G	A	100.0%	193.6 / NA	49	V338M (GTG→ATG)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/31);  total (18/31)

GCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGC  >  NC_000913/1756563‑1756852                                                                                                                                                   |                                                                                                                                                gCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCAtgat                                                                                                                                             <  1:462931/151‑1 (MQ=255)    cGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCAtgatga                                                                                                                                           >  1:260502/1‑150 (MQ=255)     gACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAAc                                                                                                                                         <  1:340797/151‑1 (MQ=255)         cgcgCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAg                                                                                                                                       <  1:257127/149‑1 (MQ=255)            gcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGtc                                                                                                                                  <  1:202229/151‑1 (MQ=255)                         aCGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTcaacaa                                                                                                                     <  2:665685/151‑1 (MQ=255)                            ttGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATg                                                                                                                   >  1:710268/1‑150 (MQ=255)                                     ccATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTaa                                                                                                         <  2:102214/151‑1 (MQ=255)                                      cATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGGATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTaa                                                                                                         >  1:696824/1‑150 (MQ=255)                                         ccTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACt                                                                                                      <  1:272525/150‑1 (MQ=255)                                         ccTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTg                                                                                                     <  1:361070/151‑1 (MQ=255)                                         ccTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTg                                                                                                     >  1:676887/1‑151 (MQ=255)                                              aCGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAt                                                                                                <  2:400152/151‑1 (MQ=255)                                                ggTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTa                                                                                              <  2:335105/151‑1 (MQ=255)                                                 gTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTAc                                                                                             <  1:729525/151‑1 (MQ=255)                                                       aCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATGCCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGaa                                                                                         <  1:604437/149‑1 (MQ=255)                                                       aCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCTAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGc                                                                                       <  2:185475/151‑1 (MQ=255)                                                       aCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGc                                                                                       <  1:712671/151‑1 (MQ=255)                                                           aGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTa                                                                                   <  2:708954/151‑1 (MQ=255)                                                             tGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATg                                                                                 >  1:745537/1‑151 (MQ=255)                                                                 gCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgt                                                                             >  1:150734/1‑151 (MQ=255)                                                                 gCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCg                                                                              <  2:607700/150‑1 (MQ=255)                                                                  cTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgtg                                                                            >  1:3026/1‑151 (MQ=255)                                                                    gTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgtggt                                                                          <  1:748876/151‑1 (MQ=255)                                                                       tGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGTCGTGGTGcc                                                                       >  1:704844/1‑151 (MQ=255)                                                                        ggTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgtggtg                                                                         <  2:3026/148‑1 (MQ=255)                                                                        ggTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGGGCCg                                                                      >  2:539884/1‑151 (MQ=255)                                                                               ccGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAAc                                                               <  1:144655/151‑1 (MQ=255)                                                                                cGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAAct                                                              >  1:19544/1‑151 (MQ=255)                                                                                   aaaaGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActgct                                                           <  2:493503/151‑1 (MQ=255)                                                                                     aaGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGa                                                         <  2:557990/151‑1 (MQ=255)                                                                                           aaTACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACt                                                   >  2:483930/1‑151 (MQ=255)                                                                                              aCCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTgg                                                 >  2:323417/1‑150 (MQ=255)                                                                                                ccGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATg                                              >  2:682406/1‑151 (MQ=255)                                                                                                    tGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTcc                                          >  1:424762/1‑151 (MQ=255)                                                                                                     ggAAGCGGTTTCTATCATGGCGACCAGCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGAGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCg                                         <  2:38231/151‑1 (MQ=255)                                                                                                        aGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTg                                      <  1:233606/151‑1 (MQ=255)                                                                                                        aGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTg                                      <  1:376049/151‑1 (MQ=255)                                                                                                          cGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGAt                                    <  2:375790/151‑1 (MQ=255)                                                                                                                    cATGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCt                          >  1:658454/1‑151 (MQ=255)                                                                                                                     aTGGCGACCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTa                         <  2:713751/151‑1 (MQ=255)                                                                                                                           aCCATCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTGCCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAgg                   <  2:136320/151‑1 (MQ=255)                                                                                                                               tCTGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCgg               >  2:227777/1‑151 (MQ=255)                                                                                                                                 tGCGAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTCCCGAAGCGGTATGCCGTGGTGCCGCTGAACCTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCGACTCAGGGc                 <  2:282986/147‑1 (MQ=255)                                                                                                                                    gAACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAAt          >  1:411678/1‑151 (MQ=255)                                                                                                                                      aCGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTCCTCAGGGCGGTAAATCt        <  2:390424/151‑1 (MQ=255)                                                                                                                                       cGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTg       <  2:424741/151‑1 (MQ=255)                                                                                                                                        gTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGc      >  2:253747/1‑151 (MQ=255)                                                                                                                                            cgaccgCATGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTcgc  <  1:97293/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGC  >  NC_000913/1756563‑1756852

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