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.C97.6% 141.7 / ‑3.6 42coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (21/20);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GG‑CCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250490‑1250754
                                                                                                                                        |                                                                                                                                  
tCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GG‑cccc                                                                                                                                >  1:370603/1‑139 (MQ=255)
                 aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCTCGGACAGGG‑CCCCCGAAAGCATCCCCTGa                                                                                                                <  2:655326/139‑1 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                         >  2:978937/1‑139 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                         >  1:483833/1‑139 (MQ=255)
                        cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                         >  1:831145/1‑139 (MQ=255)
                                tAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAg                                                                                                 <  2:8433/139‑1 (MQ=255)
                                                   gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  2:623000/107‑1 (MQ=255)
                                                   gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              >  1:623000/1‑107 (MQ=255)
                                                       acCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  1:1028747/103‑1 (MQ=255)
                                                       acCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              >  2:1028747/1‑103 (MQ=255)
                                                       acCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgcc                                                                          >  1:350526/1‑139 (MQ=255)
                                                           tCGCCGCCATGAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAg                                                                                                                          <  2:814185/87‑1 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAg                                                                                                                          >  1:814185/1‑87 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGa                                                                                                                >  1:573156/1‑97 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGa                                                                                                                <  2:573156/97‑1 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                      >  2:851130/1‑139 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GG‑c                                                                                                                                   >  1:743369/1‑77 (MQ=255)
                                                           tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GG‑c                                                                                                                                   <  2:743369/77‑1 (MQ=255)
                                                            cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAc                                                                     >  2:477859/1‑139 (MQ=255)
                                                              ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAc                                                                                                               >  1:254105/1‑95 (MQ=255)
                                                              ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAc                                                                                                               <  2:254105/95‑1 (MQ=255)
                                                              ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGGTGCCGCCa                                                                      >  1:741291/1‑136 (MQ=255)
                                                              ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                      <  2:741291/136‑1 (MQ=255)
                                                                       aCTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAgg                                                          <  2:483833/139‑1 (MQ=255)
                                                                       aCTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAgg                                                          <  2:370603/139‑1 (MQ=255)
                                                                           gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAg                                                                                                                          <  1:339949/71‑1 (MQ=255)
                                                                           gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAg                                                                                                                          >  2:339949/1‑71 (MQ=255)
                                                                              tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              >  2:730054/1‑80 (MQ=255)
                                                                              tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAcc                                                                                                              <  1:730054/80‑1 (MQ=255)
                                                                                gcgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACtttt                                                 >  1:805700/1‑139 (MQ=255)
                                                                                 cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAc                                                                                                               >  2:545118/1‑76 (MQ=255)
                                                                                 cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGAc                                                                                                               <  1:545118/76‑1 (MQ=255)
                                                                                  gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACtttcct                                               >  2:34381/1‑136 (MQ=255)
                                                                                  gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               >  2:884449/1‑139 (MQ=255)
                                                                                   cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                    >  1:836518/1‑133 (MQ=255)
                                                                                   cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                    <  2:836518/133‑1 (MQ=255)
                                                                                   cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCt                                              <  1:1005246/139‑1 (MQ=255)
                                                                                              tttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGc                                   <  1:851130/139‑1 (MQ=255)
                                                                                                tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                 >  2:196722/1‑139 (MQ=255)
                                                                                                   cGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAg                              <  2:805700/139‑1 (MQ=255)
                                                                                                         tCGGTGAGGTGAAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCggg                        <  1:196722/139‑1 (MQ=255)
                                                                                                           ggtgaggtgaAAATTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTg                      <  2:831145/139‑1 (MQ=255)
                                                                                                                     aaaTTTCGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg            <  1:978937/139‑1 (MQ=255)
                                                                                                                           cGCCCGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCt      >  1:749066/1‑139 (MQ=255)
                                                                                                                               cGGACA‑GGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATg  <  2:350526/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
TCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACA‑GG‑CCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250490‑1250754

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