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.7% 95.5 / ‑3.0 30coding (215/1071 nt)dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (14/15);  total (14/16)
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.91e-01

GCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACC  >  NC_000913/1250492‑1250743
                                                                                                                                     |                                                                                                                        
gCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCTCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGcccccc                                                                                                                   >  1:968501/1‑138 (MQ=255)
gCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGcccccc                                                                                                                   <  2:290443/139‑2 (MQ=255)
           gTTAAGAATATCGCCGGTGGAATTTTAGATAATAAACCGTAAACCTTGGCCGCATTCAACTTGCAAGGCGCATTCAAAGCTTTTATCGGGCGTCGGTAAGGTGACAATTTCGCCCGGACAGGCGCCTCCAAAGCATccc                                                                                                         <  2:968501/139‑1 (MQ=255)
                 aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                   <  1:1138124/139‑1 (MQ=255)
                    atCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGat                                                                                                <  1:560359/139‑1 (MQ=255)
                     tCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGata                                                                                               <  2:442250/139‑1 (MQ=255)
                      cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                              >  1:484696/1‑123 (MQ=255)
                      cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGc                                                                                                              <  2:484696/123‑1 (MQ=255)
                                      gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGat                                                                                                >  1:50525/1‑121 (MQ=255)
                                      gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGat                                                                                                <  2:50525/121‑1 (MQ=255)
                                                 gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGc                                                                   >  2:1273786/1‑139 (MQ=255)
                                                         tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa                                                           >  2:760696/1‑139 (MQ=255)
                                                          cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                             <  1:198644/136‑1 (MQ=255)
                                                          cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                             >  2:198644/1‑136 (MQ=255)
                                                          cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAc                                                          >  2:947628/1‑139 (MQ=255)
                                                          cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAc                                                          >  2:295846/1‑139 (MQ=255)
                                                                 aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcag                                                   >  2:958635/1‑139 (MQ=255)
                                                                         gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                   >  1:1122767/1‑83 (MQ=255)
                                                                         gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGAcc                                                                                                   <  2:1122767/83‑1 (MQ=255)
                                                                               cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                     <  2:1260880/75‑1 (MQ=255)
                                                                               cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGa                                                                                                     >  1:1260880/1‑75 (MQ=255)
                                                                                              tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAg                                                                                                               <  2:41059/50‑1 (MQ=255)
                                                                                              tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAg                                                                                                               >  1:41059/1‑50 (MQ=255)
                                                                                              tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt                      <  1:947628/139‑1 (MQ=255)
                                                                                               aTCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATc                     <  1:958635/139‑1 (MQ=255)
                                                                                                  gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc                  <  1:124020/139‑1 (MQ=255)
                                                                                                    gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc                <  1:760696/139‑1 (MQ=255)
                                                                                                                  aaaaTTTCGGCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAAc                                                                                           >  2:154617/1‑50 (MQ=255)
                                                                                                                  aaaaTTTCGGCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGATATAAc                                                                                           <  1:154617/50‑1 (MQ=255)
                                                                                                                  aaaaTTTCCCCCGGCCAGGCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACACGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACAcc  >  2:305996/1‑139 (MQ=255)
                                                                                                                    aaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                 <  2:464756/42‑1 (MQ=255)
                                                                                                                    aaTTTCGCCCGGACAGGCCCCCCGAAAGCATCCCCTGACCGa                                                                                                 >  1:464756/1‑42 (MQ=255)
                                                                                                                                     |                                                                                                                        
GCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACC  >  NC_000913/1250492‑1250743

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