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% 178.8 / ‑0.1 52coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (22/28);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

CTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250493‑1250754
                                                                                                                                    |                                                                                                                                    
cTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGccccccg                                                                                                                              <  1:370846/139‑3 (MQ=255)
     tCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                          <  1:448192/139‑1 (MQ=255)
                     cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGatat                                                                                                          >  2:843912/1‑139 (MQ=255)
                       ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATaa                                                                                                       >  1:897010/1‑139 (MQ=255)
                       ccGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATaa                                                                                                       >  1:473313/1‑139 (MQ=255)
                        cGGTTTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAAc                                                                                                      <  1:969456/139‑1 (MQ=255)
                        cGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTg                                                                                                                  <  1:8556/128‑1 (MQ=255)
                        cGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTg                                                                                                                  >  2:8556/1‑128 (MQ=255)
                          gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAc                                                                                                                >  1:32925/1‑128 (MQ=255)
                          gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAc                                                                                                                <  2:32925/128‑1 (MQ=255)
                             tAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                             <  1:1290516/128‑1 (MQ=255)
                             tAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                             >  2:1290516/1‑128 (MQ=255)
                                     gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                 >  1:599467/1‑116 (MQ=255)
                                     gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                                 <  2:599467/116‑1 (MQ=255)
                                     gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                               <  1:585288/118‑1 (MQ=255)
                                     gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                               >  2:585288/1‑118 (MQ=255)
                                     gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAc                                                                                                                <  1:719022/117‑1 (MQ=255)
                                     gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAc                                                                                                                >  2:719022/1‑117 (MQ=255)
                                                gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                          >  2:782771/1‑96 (MQ=255)
                                                gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                                          <  1:782771/96‑1 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        <  1:644498/139‑1 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        <  1:647621/139‑1 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        >  2:178776/1‑139 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                                        >  2:584273/1‑139 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTATCCGCTgccgc                                                                        >  1:851264/1‑139 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCG‑‑ATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                                      >  2:400709/1‑139 (MQ=255)
                                                      cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCAAAGGATTCCCTTGCCGATTTAAACCACATTCCTCGGGCTGTTTCCGCTgccgcc                                                                       >  1:73429/1‑139 (MQ=255)
                                                         cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCACCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAc                                                                     >  2:939764/1‑139 (MQ=255)
                                                           ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc                                                                   >  2:988859/1‑139 (MQ=255)
                                                                aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcag                                                              <  1:648700/139‑1 (MQ=255)
                                                                      ttGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGc                                                        <  2:108209/139‑1 (MQ=255)
                                                                        gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCG‑‑ATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                                    <  1:400709/139‑1 (MQ=255)
                                                                          aTGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc                                                                                           >  1:176095/1‑100 (MQ=255)
                                                                          aTGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc                                                                                           <  2:176095/100‑1 (MQ=255)
                                                                           tGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACtt                                                   <  1:594368/139‑1 (MQ=255)
                                                                               gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               <  2:1257404/139‑1 (MQ=255)
                                                                               gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                               <  1:843912/139‑1 (MQ=255)
                                                                                cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCt                                                                                       >  1:615913/1‑98 (MQ=255)
                                                                                cATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCt                                                                                       <  2:615913/98‑1 (MQ=255)
                                                                                       aGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAg                                       <  2:509725/139‑1 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCg                                                       <  2:149202/117‑1 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCg                                                       >  1:149202/1‑117 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                 <  2:994559/139‑1 (MQ=255)
                                                                                             tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                                 >  2:1225695/1‑139 (MQ=255)
                                                                                                     gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                         <  2:622478/139‑1 (MQ=255)
                                                                                                      tCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCggg                        <  1:584273/139‑1 (MQ=255)
                                                                                                           gaggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACACGACAGCAGGGCGATTTTTCTTGCAACAGGGGCATCAGCTCGGGTGACa                   <  2:851264/139‑1 (MQ=255)
                                                                                                            aggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACAt                  <  2:92682/139‑1 (MQ=255)
                                                                                                             ggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                         >  2:221536/1‑131 (MQ=255)
                                                                                                             ggtgaAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg                         <  1:221536/131‑1 (MQ=255)
                                                                                                                 aaaaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACAcc             >  2:57158/1‑139 (MQ=255)
                                                                                                                           ccGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGCGCATCGGCTCATGTCCGCTGCCGCCACCCGACAGCAGCGCGACTTTTCCTGCAACAGCGGCATCAGCTCGGGTGACATAAACCGGATCCTGa  >  2:306649/1‑138 (MQ=255)
                                                                                                                            cGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATg  <  2:897010/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATG  >  NC_000913/1250493‑1250754

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