Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,248,410	A→G	A235A (GCT→GCC)	dhaM ←	putative dihydroxyacetone‑specific PTS enzymes: HPr, EI components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,248,410	0	A	G	100.0%	194.9 / NA	54	A235A (GCT→GCC)	dhaM	putative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/26);  total (28/26)

CTTCCACGGTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCAT  >  NC_000913/1248277‑1248533                                                                                                                                      |                                                                                                                            cttcCACGGTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGTCTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        >  2:288145/1‑139 (MQ=255) cttcCACGGTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        <  2:1286618/139‑1 (MQ=255) cttcCACGGTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        <  2:138244/139‑1 (MQ=255)      aCGGTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGAc                                                                                                                   >  1:1095937/1‑139 (MQ=255)        ggTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACgg                                                                                                                 <  1:525884/139‑1 (MQ=255)              ggTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACgg                                                                                                                 <  2:1091566/133‑1 (MQ=255)              ggTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACgg                                                                                                                 >  1:1091566/1‑133 (MQ=255)                             tACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        >  1:530800/1‑111 (MQ=255)                             tACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        <  2:530800/111‑1 (MQ=255)                                  tACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGaaa                                                                                                              <  1:526853/116‑1 (MQ=255)                                  tACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGaaa                                                                                                              >  2:526853/1‑116 (MQ=255)                                       aaaaCTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGcc                                                                                  >  2:944697/1‑139 (MQ=255)                                         aaCTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCg                                                                                                           >  1:228874/1‑112 (MQ=255)                                         aaCTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCg                                                                                                           <  2:228874/112‑1 (MQ=255)                                              gTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        >  1:251355/1‑94 (MQ=255)                                              gTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        <  2:251355/94‑1 (MQ=255)                                              gTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGc                                                                           >  1:564400/1‑139 (MQ=255)                                                      taaAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTc                                                                                             <  2:1196382/113‑1 (MQ=255)                                                      taaAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTc                                                                                             >  1:1196382/1‑113 (MQ=255)                                                         aaGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTg                                                                                                    <  2:331935/103‑1 (MQ=255)                                                         aaGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTg                                                                                                    >  1:331935/1‑103 (MQ=255)                                                         aaGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCa                                                                <  1:944697/139‑1 (MQ=255)                                                             cTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATc                                                                                                        >  1:656963/1‑95 (MQ=255)                                                             cTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATc                                                                                                        <  2:656963/95‑1 (MQ=255)                                                             cTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                                            <  1:432218/139‑1 (MQ=255)                                                             cTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                                            >  1:763437/1‑139 (MQ=255)                                                              tttACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGt                                                                                                     <  2:605950/97‑1 (MQ=255)                                                              tttACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGt                                                                                                     >  1:605950/1‑97 (MQ=255)                                                              tttACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCg                                                                            <  2:37152/122‑1 (MQ=255)                                                              tttACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAAt                                                                        <  2:1334013/126‑1 (MQ=255)                                                              tttACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAAt                                                                        >  1:1334013/1‑126 (MQ=255)                                                              tttACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTGCTGGCCCTTTCg                                                                            >  1:37152/1‑122 (MQ=255)                                                                tACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCa                                                                                            <  1:597318/104‑1 (MQ=255)                                                                tACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCa                                                                                            >  2:597318/1‑104 (MQ=255)                                                                tACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                                            <  2:971624/136‑1 (MQ=255)                                                                tACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGt                                                            >  1:971624/1‑136 (MQ=255)                                                                      aaaCAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTc                                                                                             >  2:1341363/1‑97 (MQ=255)                                                                      aaaCAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTc                                                                                             <  1:1341363/97‑1 (MQ=255)                                                                        aCAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGc                                                 <  2:1095937/139‑1 (MQ=255)                                                                          aGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        <  2:1127191/66‑1 (MQ=255)                                                                          aGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTcggccag                                                                                                                        >  1:1127191/1‑66 (MQ=255)                                                                            gCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTgg                                                                                    <  2:215798/100‑1 (MQ=255)                                                                            gCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTgg                                                                                    >  1:215798/1‑100 (MQ=255)                                                                            gCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAAc                                             >  1:414053/1‑139 (MQ=255)                                                                                     gACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTaa                                    >  2:18697/1‑139 (MQ=255)                                                                                          agagTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCtt                                                                                              >  2:870447/1‑76 (MQ=255)                                                                                          agagTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCtt                                                                                              <  1:870447/76‑1 (MQ=255)                                                                                          agagTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGgcg                                                                  >  2:126963/1‑104 (MQ=255)                                                                                          agagTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGgcg                                                                  <  1:126963/104‑1 (MQ=255)                                                                                           gagTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAAt                              <  2:1332090/139‑1 (MQ=255)                                                                                           gagTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAAt                              <  1:18697/139‑1 (MQ=255)                                                                                                    gAGCGTCTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAAt                     <  1:288145/139‑1 (MQ=255)                                                                                                     aGCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATa                    <  2:763437/139‑1 (MQ=255)                                                                                                      gCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACgg                                                                                                                 >  1:1090684/1‑45 (MQ=255)                                                                                                      gCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACgg                                                                                                                 <  2:1090684/45‑1 (MQ=255)                                                                                                      gCGACTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATAc                   >  1:495740/1‑139 (MQ=255)                                                                                                         aCTTCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATActct                >  1:179690/1‑139 (MQ=255)                                                                                                           ttCCTCCGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATActctct              >  1:134162/1‑139 (MQ=255)                                                                                                            tcctccGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGc                                                 >  1:100310/1‑103 (MQ=255)                                                                                                            tcctccGTTTCACCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGc                                                 <  2:100310/103‑1 (MQ=255)                                                                                                                       aCCAAAGTTATCTTCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCAt  >  2:224496/1‑139 (MQ=255)                                                                                                                                   ttCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtag                                       >  1:1039928/1‑90 (MQ=255)                                                                                                                                   ttCGGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTgtag                                       <  2:1039928/90‑1 (MQ=255)                                                                                                                                      |                                                                                                                            CTTCCACGGTCAGGGTTGATTTTGCCTGTACCGTACATAAAACTGGTTGATAATAAAAGGCTTTACCCGAAACAGGCGGAACGGGACGCAGAGTAGGTGGAGCGACTTCCTCCGTTTCACCAAAGTTATCTTCAGCCAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCAT  >  NC_000913/1248277‑1248533

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