Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 83,964 C→T 56.4% G274G (GGC→GGT ftsW → integral membrane protein involved in stabilising FstZ ring during cell division

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE83,9640CT56.4% 29.0 / 93.2 117G274G (GGC→GGTftsWintegral membrane protein involved in stabilising FstZ ring during cell division
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base T (55/11);  total (82/35)
Fisher's exact test for biased strand distribution p-value = 5.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

CCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCA  >  minE/83899‑84026
                                                                 |                                                              
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:862474/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:1618889/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:1538443/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:1835293/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:1958767/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:1002597/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                                             <  1:2276695/69‑1 (MQ=255)
ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGGCGCGGCGaa                                                             <  1:460399/69‑1 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACtg                                                          >  1:3081069/1‑57 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTggg                                                      >  1:912152/1‑62 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2091244/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2022332/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:3163902/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2988065/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:1976847/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2402882/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:213836/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:3184084/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:1707862/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:645831/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2689384/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:1032817/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:1041866/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2863094/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:888530/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:828528/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:1377736/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:2742598/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                               >  1:1621978/1‑69 (MQ=255)
              cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAg                                                >  1:3214470/1‑68 (MQ=255)
                               aGTTAACGCAATCGTTGATGGCGTTTGTTCGCGGTGa                                                              >  1:2506548/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:285133/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:3363617/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:1576975/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:2670980/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:612853/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:1998255/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:2797751/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:308580/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:1946045/1‑37 (MQ=37)
                               aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                                              >  1:1916765/1‑37 (MQ=37)
                               aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGt                                >  1:1387973/1‑67 (MQ=255)
                               aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                               >  1:2970635/1‑68 (MQ=255)
                               aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                               >  1:3028158/1‑68 (MQ=255)
                               aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                               >  1:393344/1‑68 (MQ=255)
                               aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                               >  1:2172894/1‑68 (MQ=255)
                                  tAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACaaa                           <  1:1304616/69‑1 (MQ=255)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:3190245/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:236951/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:201992/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:2482609/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:1237706/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:455061/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                            <  1:626565/51‑1 (MQ=37)
                                   aaCGCAATCGTTGATGGCGTTTGGTCGCGATGAACTGTGGGGACAAGGttt                                            <  1:2628487/51‑1 (MQ=255)
                                          tCGCTGATGGGGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGa                    >  1:1587064/1‑68 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:588596/1‑69 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:1098543/1‑69 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:114181/1‑69 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:1348281/1‑69 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:69436/1‑69 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:2610414/1‑69 (MQ=255)
                                          tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                   >  1:3099613/1‑69 (MQ=255)
                                            gCTGATGGCGTTTTGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:340490/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:3119967/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:3385140/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:978898/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:2452635/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:3038392/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:1864786/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:1346005/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:2143900/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:2238389/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:2397804/69‑1 (MQ=255)
                                            gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:2436548/69‑1 (MQ=255)
                                             cTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                 <  1:2488172/68‑1 (MQ=255)
                                                       ttGGTCTCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:5261/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGTTACAAAAACTGGAGTATCTGCCGGAAg      >  1:2322288/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1298716/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:3348674/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:2201737/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1079741/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:2625173/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1802276/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1958668/1‑69 (MQ=255)
                                                       ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1807511/1‑69 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGc            >  1:549283/1‑62 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1663652/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1536738/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1016419/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:71295/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1371736/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1788525/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:491518/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:1974028/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:2093404/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:339005/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:3126400/1‑68 (MQ=255)
                                                        tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg      >  1:2604186/1‑68 (MQ=255)
                                                             gcggcgAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGca  >  1:2928866/1‑67 (MQ=255)
                                                                 |                                                              
CCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCA  >  minE/83899‑84026

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