Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 83,944 C→T 49.5% L268L (CTG→TTG)  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,9440CT49.5% 2.3 / 98.4 107L268L (CTG→TTG) ftsWintegral membrane protein involved in stabilising FstZ ring during cell division
Reads supporting (aligned to +/- strand):  ref base C (22/32);  new base T (42/11);  total (64/43)
Fisher's exact test for biased strand distribution p-value = 6.70e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 108 total aligned reads displayed.

CCGTGTTACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTA  >  minE/83877‑84011
                                                                   |                                                                   
ccgTGTTACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    <  1:1120252/69‑1 (MQ=255)
ccgTGTTACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    <  1:3183570/69‑1 (MQ=255)
     ttACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    <  1:9705/64‑1 (MQ=255)
     ttACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    <  1:2757739/64‑1 (MQ=255)
     ttACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    >  1:1299951/1‑64 (MQ=255)
     ttACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    >  1:1770627/1‑64 (MQ=255)
     ttACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    <  1:2423947/64‑1 (MQ=255)
     ttACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACCCAATCGCt                                                                    <  1:2473811/64‑1 (MQ=255)
       aCCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCg                                                             >  1:3080845/1‑69 (MQ=255)
       aCCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCg                                                             >  1:2988890/1‑69 (MQ=255)
       aCCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCg                                                             >  1:2913062/1‑69 (MQ=255)
       aCCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCg                                                             >  1:3330965/1‑69 (MQ=255)
       aCCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCg                                                             >  1:910910/1‑69 (MQ=255)
           cATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTg                                                         <  1:478/69‑1 (MQ=255)
               cTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCt                                                                    <  1:3230130/54‑1 (MQ=255)
               cTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTcg                                                     >  1:1535600/1‑69 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:1538443/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:1618889/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:2276695/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:1002597/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:1835293/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:862474/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:797777/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGaa                                              <  1:1958767/69‑1 (MQ=255)
                      ccGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGGCGCGGCGaa                                              <  1:460399/69‑1 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCgcgg                                                  >  1:2046222/1‑51 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACtg                                           >  1:3081069/1‑57 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTggg                                       >  1:912152/1‑62 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:2449906/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:2742598/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:2689384/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:888530/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:1377736/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:3126020/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:645831/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:1621978/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:2402882/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:3163902/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:1707862/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:1976847/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:3184084/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:1993843/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:2022332/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:213836/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:2091244/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                >  1:828528/1‑69 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAg                                 >  1:2321157/1‑68 (MQ=255)
                                    cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAg                                 >  1:3214470/1‑68 (MQ=255)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:2621106/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:3363617/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:2670980/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:612853/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:1576975/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:2797751/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:308580/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:285133/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:1916765/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:428686/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:1998255/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:1946045/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGa                                               >  1:170877/1‑37 (MQ=37)
                                                     aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                >  1:393344/1‑68 (MQ=255)
                                                     aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                >  1:3028158/1‑68 (MQ=255)
                                                     aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                >  1:2970635/1‑68 (MQ=255)
                                                     aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                >  1:2172894/1‑68 (MQ=255)
                                                        tAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACaaa            <  1:1304616/69‑1 (MQ=255)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:236951/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:1128627/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:1237706/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:455061/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:201992/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:588438/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:626565/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:2482609/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                             <  1:2821219/51‑1 (MQ=37)
                                                         aaCGCAATCGTTGATGGCGTTTGGTCGCGATGAACTGTGGGGACAAGGttt                             <  1:2628487/51‑1 (MQ=255)
                                                                tCGCTGATGGGGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGa     >  1:1587064/1‑68 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:69436/1‑69 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:2610414/1‑69 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:1098543/1‑69 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:114181/1‑69 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:3099613/1‑69 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:1348281/1‑69 (MQ=255)
                                                                tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg    >  1:2131759/1‑69 (MQ=255)
                                                                  gCTGATGGCGTTTTGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:340490/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:2397804/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:978898/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:1346005/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:1864786/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:2143900/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:2436548/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:2452635/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:3038392/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:3119967/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:3287423/69‑1 (MQ=255)
                                                                  gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:3385140/69‑1 (MQ=255)
                                                                   cTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa  <  1:2488172/68‑1 (MQ=255)
                                                                   |                                                                   
CCGTGTTACCGCATTCTGGAACCCGTGGGAAGATCCCTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTA  >  minE/83877‑84011

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: