Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 89,165 A→G 14.1% G148G (GGA→GGG ftsA → ATP‑binding cell division protein involved in recruitment of FtsK to Z ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE89,1650AG14.1% 344.8 / 33.4 196G148G (GGA→GGGftsAATP‑binding cell division protein involved in recruitment of FtsK to Z ring
Reads supporting (aligned to +/- strand):  ref base A (88/80);  new base G (9/19);  total (97/99)
Fisher's exact test for biased strand distribution p-value = 6.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 197 total aligned reads displayed.

GTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAA  >  minE/89098‑89230
                                                                   |                                                                 
gtgtGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAc                                                                  >  1:960498/1‑69 (MQ=255)
gtgtGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAc                                                                  >  1:2900119/1‑69 (MQ=255)
            tgATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTg                                                       >  1:3323842/1‑68 (MQ=255)
            tgATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTg                                                       >  1:1570142/1‑68 (MQ=255)
            tgATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTg                                                       >  1:626653/1‑68 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:2224210/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:251818/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:2696441/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:288312/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:2926583/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:3118558/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:336309/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:835837/68‑1 (MQ=255)
              aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCg                                                     <  1:1195745/68‑1 (MQ=255)
                         gTATGCGATTGACTATCAGGACGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaa                                         <  1:1454895/69‑1 (MQ=255)
                         gTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaa                                         <  1:2677865/69‑1 (MQ=255)
                           aTGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       <  1:474339/69‑1 (MQ=255)
                           aTGCGATTGACTATCAGGAAGAGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       <  1:675530/69‑1 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:513220/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       <  1:2626108/68‑1 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:1089525/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:105926/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:1129169/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:849361/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:1933552/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:1226781/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:339657/1‑68 (MQ=255)
                            tGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       >  1:1563131/1‑68 (MQ=255)
                              cGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACt                                                                 <  1:1212406/40‑1 (MQ=255)
                              cGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACt                                                                 >  1:453152/1‑40 (MQ=255)
                              cGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACt                                                                 >  1:463868/1‑40 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:515636/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:2471252/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:1025521/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:2712725/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:2815705/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:3049187/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTg                               <  1:1426037/69‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaa                                          <  1:869926/58‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaa                                          <  1:23281/58‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaa                                          >  1:1880673/1‑58 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaa                                          <  1:2194620/58‑1 (MQ=255)
                                   gACTATCAGGAAGGGATCAAGAATACGGTAGGACTTTCGGGCGTGCGGATGCAGGCaa                                          >  1:403841/1‑58 (MQ=255)
                                     cTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                              <  1:1296594/68‑1 (MQ=255)
                                     cTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                              <  1:1311497/68‑1 (MQ=255)
                                     cTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                              <  1:3097147/68‑1 (MQ=255)
                                     cTATCAGGAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                              <  1:3059716/68‑1 (MQ=255)
                                      tATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                             <  1:1019079/68‑1 (MQ=255)
                                      tATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                             <  1:2578048/68‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:2606689/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:516494/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:1529327/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:24413/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:3334031/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:3247167/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:3213617/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:3176435/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:3133259/69‑1 (MQ=255)
                                       aTCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           <  1:2150248/69‑1 (MQ=255)
                                          aGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcaca                         <  1:2017883/68‑1 (MQ=255)
                                          aGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcaca                         <  1:3166212/68‑1 (MQ=255)
                                          aGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcaca                         <  1:2208539/68‑1 (MQ=255)
                                           ggAAGGGATCAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                              <  1:2609321/62‑1 (MQ=255)
                                              aGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                     >  1:225888/1‑68 (MQ=255)
                                              aGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                     >  1:1595463/1‑68 (MQ=255)
                                              aGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                     >  1:2426228/1‑68 (MQ=255)
                                              aGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                     >  1:561455/1‑68 (MQ=255)
                                              aGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                     >  1:786209/1‑68 (MQ=255)
                                              aGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTc                     >  1:1176839/1‑68 (MQ=255)
                                              aGGGATAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTcac                   >  1:2286131/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcac                          >  1:593164/1‑57 (MQ=255)
                                                    cAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcac                          >  1:897279/1‑57 (MQ=255)
                                                    cAAGAATCCGGTAGGGCTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATca                           >  1:3039387/1‑56 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGTGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:1573112/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTGACAACGa              >  1:325076/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:743216/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:135849/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:2347354/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2522314/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2855622/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:419127/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:3428710/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2025251/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2914149/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:1912466/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:3016772/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2773166/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2747936/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              >  1:1789367/1‑69 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:1737063/69‑1 (MQ=255)
                                                    cAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa              <  1:2637483/69‑1 (MQ=255)
                                                     aaGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAg                                       <  1:2196683/43‑1 (MQ=255)
                                                       gAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaa                                          <  1:204717/38‑1 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa  >  1:2090412/1‑69 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa  >  1:1380940/1‑69 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa  >  1:2138087/1‑69 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa  >  1:2492186/1‑69 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa  >  1:2382558/1‑69 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa  >  1:1713068/1‑69 (MQ=255)
                                                                aGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaa   >  1:1894221/1‑68 (MQ=255)
                                                                   |                                                                 
GTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAA  >  minE/89098‑89230

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