Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,250,964 G→A 12.2% intergenic (‑126/‑102) dhaK ← / → dhaR dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit/dhaKLM operon transcription activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,250,9640GA12.2% 136.1 / 5.7 48intergenic (‑126/‑102)dhaK/dhaRdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit/dhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base A (3/3);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGC  >  NC_000913/1250819‑1251098
                                                                                                                                                 |                                                                                                                                      
tcatTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGa                                                                                                                                     >  1:220115/1‑149 (MQ=255)
       tCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACatat                                                                                                                              >  1:46097/1‑149 (MQ=255)
         aaTTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATg                                                                                                                             <  3:351225/148‑1 (MQ=255)
              tttCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTaa                                                                                                                       <  2:220115/149‑1 (MQ=255)
                tCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACataat                                                                                                                             <  1:130848/141‑3 (MQ=255)
                tCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACataat                                                                                                                             >  2:130848/1‑139 (MQ=255)
                       tttGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTc                                                                                                              <  2:46097/149‑1 (MQ=255)
                         tGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCa                                                                                                            <  1:224448/149‑1 (MQ=255)
                            tCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACataatttgg                                                                                                                         <  1:108203/133‑7 (MQ=255)
                            tCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACataatttgg                                                                                                                         >  2:108203/1‑127 (MQ=255)
                              cAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCata                                                                                                          <  2:27523/146‑1 (MQ=255)
                              cAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCata                                                                                                          >  1:27523/1‑146 (MQ=255)
                                gCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGa                                                                                                      >  4:59554/1‑148 (MQ=255)
                                  aaTTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAAc                                                                                                   >  2:201086/1‑149 (MQ=255)
                                         gTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACg                                                                                            >  2:216551/1‑149 (MQ=255)
                                             ggCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTa                                                                                        >  1:14364/1‑149 (MQ=255)
                                                     tgatgTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACt                                                                                >  1:360159/1‑149 (MQ=255)
                                                             aaCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTg                                                                        >  2:240803/1‑149 (MQ=255)
                                                              aCGACACGGCGACCACTGACTGCCGATGAATCCATTTTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGc                                                                       <  1:23145/149‑1 (MQ=255)
                                                               cGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGc                                                                       >  3:34000/1‑148 (MQ=255)
                                                                gACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCg                                                                      <  3:59554/148‑1 (MQ=255)
                                                                     ggCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCg                                                                 <  3:257417/148‑1 (MQ=255)
                                                                      gCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                               >  1:235611/1‑149 (MQ=255)
                                                                         aCCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGTAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAac                                                             >  4:6123/1‑148 (MQ=255)
                                                                         aCCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAac                                                             >  4:326516/1‑148 (MQ=255)
                                                                           cACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAacac                                                           >  4:173405/1‑148 (MQ=255)
                                                                               gactgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCtttt                                                      <  1:201086/149‑1 (MQ=255)
                                                                                 ctgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTg                                                     <  3:326516/148‑1 (MQ=255)
                                                                                 ctgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGt                                                    <  2:104402/149‑1 (MQ=255)
                                                                                  tgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTc                                                   >  1:167708/1‑149 (MQ=255)
                                                                                  tgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTc                                                   >  1:388963/1‑149 (MQ=255)
                                                                                      gATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATa                                                >  4:241192/1‑148 (MQ=255)
                                                                                        tGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATa                                                <  3:262591/146‑1 (MQ=255)
                                                                                        tGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAg                                              <  4:152267/148‑1 (MQ=255)
                                                                                        tGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTCAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATa                                                >  4:262591/1‑146 (MQ=255)
                                                                                           aTCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATAc                                                                                 <  1:208636/110‑1 (MQ=255)
                                                                                           aTCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATAc                                                                                 >  2:208636/1‑110 (MQ=255)
                                                                                                ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                     >  1:81504/1‑149 (MQ=255)
                                                                                                 tgtgCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                     >  3:132403/1‑148 (MQ=255)
                                                                                                 tgtgCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                     <  4:140102/148‑1 (MQ=255)
                                                                                                      aTCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGa                               >  1:216480/1‑149 (MQ=255)
                                                                                                           gCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGgc                          >  1:54280/1‑149 (MQ=255)
                                                                                                                ggAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGt                             <  3:69684/141‑1 (MQ=255)
                                                                                                                ggAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGt                             >  4:69684/1‑141 (MQ=255)
                                                                                                                                ccAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCg      <  4:132403/148‑1 (MQ=255)
                                                                                                                                   aaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATTGTCGGGGc  >  2:334444/1‑149 (MQ=255)
                                                                                                                                     atATGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTa                                                                                              >  1:176981/2‑55 (MQ=37)
                                                                                                                                     atATGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTa                                                                                              <  2:176981/54‑1 (MQ=37)
                                                                                                                                                 |                                                                                                                                      
TCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGC  >  NC_000913/1250819‑1251098

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