Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,250,629 (C)5→6 coding (210/1071 nt) dhaK ← dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,250,6241.C95.8% 158.6 / ‑0.3 48coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (19/27);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250488‑1250754
                                                                                                                                          |                                                                                                                                   
ggTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGcc                                                                                                                                  >  1:790302/1‑139 (MQ=255)
  tCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGcccc                                                                                                                                >  1:708101/1‑139 (MQ=255)
                   cgAATATCGCCGGTGTAATTTTTGATAATGAACAGTACACCTTCGCCGCCATCAGCGTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                <  1:55292/138‑1 (MQ=255)
                   aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                <  2:137026/139‑1 (MQ=255)
                   aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCTCCCCGAAAGCATCCCCTGa                                                                                                                <  1:792594/139‑1 (MQ=255)
                     aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  2:708101/139‑1 (MQ=255)
                     aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  2:550406/139‑1 (MQ=255)
                         tCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGata                                                                                                          <  2:361605/139‑1 (MQ=255)
                                                      tACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCt                                                                             >  2:825582/1‑139 (MQ=255)
                                                        cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                         <  2:936538/93‑1 (MQ=255)
                                                        cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                         >  1:936538/1‑93 (MQ=255)
                                                        cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                                           >  2:1373144/1‑139 (MQ=255)
                                                        cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                                           >  1:1291608/1‑139 (MQ=255)
                                                        cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                    >  2:1173277/1‑82 (MQ=255)
                                                        cacCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                    <  1:1173277/82‑1 (MQ=255)
                                                         acCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgcc                                                                          >  2:1093357/1‑139 (MQ=255)
                                                          ccTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGACCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccg                                                                         >  1:263408/1‑139 (MQ=255)
                                                           cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              >  2:1114866/1‑101 (MQ=255)
                                                           cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  1:1114866/101‑1 (MQ=255)
                                                             tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACca                                                                                                    <  1:319521/109‑1 (MQ=255)
                                                             tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACca                                                                                                    >  2:319521/1‑109 (MQ=255)
                                                              cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                         >  1:19358/1‑87 (MQ=255)
                                                              cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                         <  2:19358/87‑1 (MQ=255)
                                                                ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc                                                                   >  2:1024121/1‑139 (MQ=255)
                                                                              cATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGAc                                                     >  2:1123212/1‑139 (MQ=255)
                                                                                   cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTc                                                >  1:743453/1‑139 (MQ=255)
                                                                                    gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               <  1:1024121/139‑1 (MQ=255)
                                                                                    gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               >  2:32166/1‑139 (MQ=255)
                                                                                    gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               <  2:261929/139‑1 (MQ=255)
                                                                                     cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  2:565949/75‑1 (MQ=255)
                                                                                     cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              >  1:565949/1‑75 (MQ=255)
                                                                                        tCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                             >  2:1414988/1‑54 (MQ=255)
                                                                                        tCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                             <  1:1414988/56‑3 (MQ=255)
                                                                                             gATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAgg                                      <  1:825582/139‑1 (MQ=255)
                                                                                                  tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                 >  1:382396/1‑139 (MQ=255)
                                                                                                  tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                 >  1:354576/1‑139 (MQ=255)
                                                                                                  tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                    <  1:521525/40‑1 (MQ=255)
                                                                                                  tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                    >  2:1366568/1‑40 (MQ=255)
                                                                                                  tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                    >  2:521525/1‑40 (MQ=255)
                                                                                                  tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                    <  1:1366568/40‑1 (MQ=255)
                                                                                                      gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                             <  1:1093357/139‑1 (MQ=255)
                                                                                                      gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                             <  2:354576/139‑1 (MQ=255)
                                                                                                       ggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCt                            <  2:263408/139‑1 (MQ=255)
                                                                                                        gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGCCCGatat                                                                                                         >  2:1087749/1‑61 (MQ=255)
                                                                                                        gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                         <  1:1087749/61‑1 (MQ=255)
                                                                                                        gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           >  2:600466/1‑139 (MQ=255)
                                                                                                        gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                           <  2:382396/139‑1 (MQ=255)
                                                                                                          gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                         <  1:1373144/139‑1 (MQ=255)
                                                                                                               tgaggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGAc                    <  1:600466/139‑1 (MQ=255)
                                                                                                                  ggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAAc                                                                                                      <  1:117957/54‑1 (MQ=255)
                                                                                                                  ggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAAc                                                                                                      >  2:117957/1‑54 (MQ=255)
                                                                                                                  ggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATa                 <  1:32166/139‑1 (MQ=255)
                                                                                                                    tgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              >  2:1016059/1‑44 (MQ=255)
                                                                                                                    tgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  1:1016059/44‑1 (MQ=255)
                                                                                                                    tgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATaca               <  2:743453/139‑1 (MQ=255)
                                                                                                                       aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg            <  1:59663/139‑1 (MQ=255)
                                                                                                                       aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg            <  2:1137983/139‑1 (MQ=255)
                                                                                                                                 cGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATg  <  2:1240715/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                   
GGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250488‑1250754

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