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

TCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGA  >  NC_000913/1250486‑1250752
                                                                                                                                           |                                                                                                                                 
tCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAgg                                                                                                                                    <  1:246725/139‑1 (MQ=255)
 cGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGc                                                                                                                                  >  2:407242/1‑139 (MQ=255)
   gTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGAGAGGTGAAAATTTCGCCCGGACAGGccc                                                                                                                                >  1:729196/1‑139 (MQ=255)
           ttCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAg                                                                                                                        <  2:251218/139‑1 (MQ=255)
                     aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                              <  2:444290/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            <  2:26808/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            <  2:680199/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            <  1:728784/139‑1 (MQ=255)
                       aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTAGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            <  2:611623/139‑1 (MQ=255)
                              ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            >  1:667732/1‑132 (MQ=255)
                              ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            <  2:667732/132‑1 (MQ=255)
                                   gtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACcaca                                                                                                <  2:260449/139‑1 (MQ=255)
                                                      aGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCg                                                                             <  2:24862/139‑1 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAc                                                                                                             <  2:488614/100‑1 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAc                                                                                                             >  1:488614/1‑100 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                      >  2:192937/1‑139 (MQ=255)
                                                             cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                      <  1:799981/139‑1 (MQ=255)
                                                               tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCTAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                    >  2:567229/1‑139 (MQ=255)
                                                                                                    tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                               >  1:291161/1‑139 (MQ=255)
                                                                                                          gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            <  2:155602/56‑1 (MQ=255)
                                                                                                          gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                            >  1:155602/1‑56 (MQ=255)
                                                                                                              cGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGt                     <  1:376871/139‑1 (MQ=255)
                                                                                                                   aggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGACTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACAt                <  1:459194/139‑1 (MQ=255)
                                                                                                                        aaaaTTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  2:78391/139‑1 (MQ=255)
                                                                                                                         tgttttCGACCGGGCAGGCCCCCCGAAAGCATCTCCTGACCGATATAACCACAGTGCATCGGGTCGTGTCCGCGGCCGCCTCACGACAGCAGGGCGACTTTTCCTGCAACATGGACAGCAGCGCGGGTGACATATACCg          <  1:716784/136‑1 (MQ=255)
                                                                                                                         aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  2:454933/139‑1 (MQ=255)
                                                                                                                         aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  1:192937/139‑1 (MQ=255)
                                                                                                                         aaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAAAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg          <  1:784279/139‑1 (MQ=255)
                                                                                                                             ttCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATc      <  2:395894/139‑1 (MQ=255)
                                                                                                                             ttCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATc      <  2:700435/139‑1 (MQ=255)
                                                                                                                                 cctggACAGGCCCGCCGAAAGCAACCCCTGAAAGATACAACCACAGGGCATCGGCTCGTCTCCGCTGCCGCCACCCGACAGAAGAGCGACTTTTCCTGCAACAGGGGCAACAGCTCGGGTGACAGAAAACGGATCCGGa  <  1:567229/136‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                 
TCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGA  >  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: