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.C96.2% 83.7 / ‑3.2 26coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (6/19);  total (6/20)
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.88e-01

TCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGA  >  NC_000913/1250486‑1250752
                                                                                                                                           |                                                                                                                                 
tCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAgg                                                                                                                                    <  1:246725/139‑1 (MQ=255)
 cGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑c                                                                                                                                  >  2:407242/1‑139 (MQ=255)
   gTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGAGAGGTGAAAATTTCGCCCGGACAGG‑ccc                                                                                                                                >  1:729196/1‑139 (MQ=255)
           ttCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAg                                                                                                                        <  2:251218/139‑1 (MQ=255)
                     aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGa                                                                                                              <  2:444290/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            <  2:26808/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            <  2:680199/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            <  1:728784/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTAGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            <  2:611623/139‑1 (MQ=255)
                              ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            >  1:667732/1‑132 (MQ=255)
                              ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            <  2:667732/132‑1 (MQ=255)
                                   gtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACcaca                                                                                                <  2:260449/139‑1 (MQ=255)
                                                      aGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCg                                                                             <  2:24862/139‑1 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAc                                                                                                             <  2:488614/100‑1 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAc                                                                                                             >  1:488614/1‑100 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                      >  2:192937/1‑139 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                      <  1:799981/139‑1 (MQ=255)
                                                               tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCTAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                    >  2:567229/1‑139 (MQ=255)
                                                                                                    tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                               >  1:291161/1‑139 (MQ=255)
                                                                                                          gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            <  2:155602/56‑1 (MQ=255)
                                                                                                          gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGAcc                                                                                                            >  1:155602/1‑56 (MQ=255)
                                                                                                              cGGTGAGGTGAAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGt                     <  1:376871/139‑1 (MQ=255)
                                                                                                                   aggtgaAAATTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGACTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACAt                <  1:459194/139‑1 (MQ=255)
                                                                                                                        aaaaTTTCGCCCGGACAGG‑CCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  2:78391/139‑1 (MQ=255)
                                                                                                                         tgttttCGACCGGGCAGGCCCC‑CCGAAAGCATCTCCTGACCGATATAACCACAGTGCATCGGGTCGTGTCCGCGGCCGCCTCACGACAGCAGGGCGACTTTTCCTGCAACATGGACAGCAGCGCGGGTGACATATACCg          <  1:716784/136‑1 (MQ=255)
                                                                                                                         aaaTTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  2:454933/139‑1 (MQ=255)
                                                                                                                         aaaTTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  1:192937/139‑1 (MQ=255)
                                                                                                                         aaaTTTCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAAAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  1:784279/139‑1 (MQ=255)
                                                                                                                             ttCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATc      <  2:395894/139‑1 (MQ=255)
                                                                                                                             ttCGCCCGGACAGGCCCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATc      <  2:700435/139‑1 (MQ=255)
                                                                                                                                 cctggACAGG‑CCCGCCGAAAGCAACCCCTGAAAGATACAACCACAGGGCATCGGCTCGTCTCCGCTGCCGCCACCCGACAGAAGAGCGACTTTTCCTGCAACAGGGGCAACAGCTCGGGTGACAGAAAACGGATCCGGa  <  1:567229/136‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                 
TCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CCC‑CCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGA  >  NC_000913/1250486‑1250752

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