Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,937,048 T→C 7.8% F133L (TTT→CTT)  yebK → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,937,0480TC 7.8% 146.9 / 3.5 51F133L (TTT→CTT) yebKputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGA  >  NC_000913/1936906‑1937195
                                                                                                                                               |                                                                                                                                                   
cAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCg                                                                                                                                               >  4:225156/1‑149 (MQ=255)
     aaGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTa                                                                                                                                          <  1:4942/149‑1 (MQ=255)
            cAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGCCTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCtt                                                                                                                                   >  1:161976/1‑149 (MQ=255)
                 ttGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAgcc                                                                                                                              <  3:225156/149‑1 (MQ=255)
                      tCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGt                                                                                                                         >  8:156464/1‑149 (MQ=255)
                         tACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGc                                                                                                                      <  4:32692/149‑1 (MQ=255)
                            acaGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGccc                                                                                                                    <  7:224357/149‑1 (MQ=255)
                            acaGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTACGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCa                                                                                                                   >  6:13211/1‑149 (MQ=255)
                              aGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAc                                                                                                                  >  1:200855/1‑148 (MQ=255)
                                 gAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCGTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatg                                                                                                              >  1:89680/1‑149 (MQ=255)
                                   aaaTATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatgc                                                                                                             <  2:81624/148‑1 (MQ=255)
                                    aaTATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatgcga                                                                                                           <  2:200855/149‑1 (MQ=255)
                                    aaTATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatgcga                                                                                                           <  2:209706/149‑1 (MQ=255)
                                     atatTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatgcgat                                                                                                          <  1:135113/149‑1 (MQ=255)
                                        tttGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGaa                                                                                                       <  1:135709/149‑1 (MQ=255)
                                              tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGtt                                                                                                 >  5:156531/1‑149 (MQ=255)
                                                   aaTGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTc                                                                                            >  5:104106/1‑149 (MQ=255)
                                                                   cATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCggt                                                                            >  8:24628/1‑149 (MQ=255)
                                                                    aTGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCggtg                                                                           <  5:155503/149‑1 (MQ=255)
                                                                              ttCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTcc                                                                  >  2:212126/1‑148 (MQ=255)
                                                                                 aCTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgatg                                                              <  4:59338/149‑1 (MQ=255)
                                                                                 aCTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGGTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgat                                                               >  1:233711/1‑148 (MQ=255)
                                                                                     gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGatat                                                          >  1:160897/1‑149 (MQ=255)
                                                                                     gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGatat                                                          >  6:152448/1‑149 (MQ=255)
                                                                                     gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGata                                                           <  6:156531/148‑1 (MQ=255)
                                                                                     gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCCTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGatat                                                          >  8:154370/1‑149 (MQ=255)
                                                                                     gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCCTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGatat                                                          >  8:150739/1‑149 (MQ=255)
                                                                                      aTAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTATTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATc                                                         >  6:195250/1‑149 (MQ=255)
                                                                                               ccATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCTGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGtgctgc                                                  <  2:163638/147‑1 (MQ=255)
                                                                                               ccATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGtgctgc                                                  <  2:160172/147‑1 (MQ=255)
                                                                                               ccATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCaa                                                >  7:14386/1‑149 (MQ=255)
                                                                                                cATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAAc                                               <  8:14386/149‑1 (MQ=255)
                                                                                                 aTCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACg                                              >  1:191165/1‑149 (MQ=255)
                                                                                                  tCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACg                                              >  2:15686/1‑148 (MQ=255)
                                                                                                         cgcCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatgcga                                                                                                           <  1:48261/80‑1 (MQ=255)
                                                                                                         cgcCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCAcgatgcga                                                                                                           >  2:48261/1‑80 (MQ=255)
                                                                                                             gTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTATTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTAt                                  <  5:195250/149‑1 (MQ=255)
                                                                                                             gTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCCTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTAt                                  <  7:154370/149‑1 (MQ=255)
                                                                                                             gTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCCTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTAt                                  <  7:150739/149‑1 (MQ=255)
                                                                                                                gACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCTACGCATGAGTTGTATGaa                               >  3:23522/1‑149 (MQ=255)
                                                                                                                 aCTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAAt                              >  4:139273/1‑149 (MQ=255)
                                                                                                                  cTTGCTCACTCAGGC‑AAAAAAAATCTCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACTCATGAGTTGTATGAAtt                             <  2:161976/149‑1 (MQ=255)
                                                                                                                      ctcactcaGGC‑AAAAAATATCGCCTATTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAATTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTa                          <  2:89680/148‑1 (MQ=255)
                                                                                                                      ctcactcaGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTa                          <  2:160897/148‑1 (MQ=255)
                                                                                                                           tcaGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAg                                        >  1:192766/1‑129 (MQ=255)
                                                                                                                           tcaGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAg                                        <  2:192766/129‑1 (MQ=255)
                                                                                                                                  aaaaaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACg              >  1:26717/1‑149 (MQ=255)
                                                                                                                                    aaaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACgtg            >  5:17380/1‑149 (MQ=255)
                                                                                                                                    aaaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACgtg            <  1:15686/149‑1 (MQ=255)
                                                                                                                                     aaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCga                                                                <  3:183924/96‑1 (MQ=255)
                                                                                                                                              cTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGa  >  5:29446/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAATGAAGATGACAGCGTTGAATCATACACAGGGAAAATATTTGAGTCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGC‑AAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGA  >  NC_000913/1936906‑1937195

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