Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 83,976 G→A 53.8% G278G (GGG→GGA 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,9760GA53.8% 9.1 / 79.6 93G278G (GGG→GGAftsWintegral membrane protein involved in stabilising FstZ ring during cell division
Reads supporting (aligned to +/- strand):  ref base G (28/15);  new base A (39/11);  total (67/26)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGC  >  minE/83913‑84044
                                                               |                                                                    
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2742598/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1707862/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1976847/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1993843/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2022332/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2091244/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:213836/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2402882/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2449906/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2689384/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1032817/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2863094/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:2988065/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:3126020/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:3163902/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:3184084/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:645831/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:828528/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:888530/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:992039/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1621978/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1377736/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAgg                                                                 >  1:1041866/1‑69 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAg                                                                  >  1:2321157/1‑68 (MQ=255)
cTTTGGCAGCGGCTATCAGTTAACGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAg                                                                  >  1:3214470/1‑68 (MQ=255)
                 aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGt                                                  >  1:1387973/1‑67 (MQ=255)
                 aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                                                 >  1:3028158/1‑68 (MQ=255)
                 aGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTa                                                 >  1:2970635/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:1128627/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:1237706/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:2482609/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:588438/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:236951/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:2821219/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:626565/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:3190245/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:201992/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGGTGAACTGTGGGGACAAGGttt                                                              <  1:455061/51‑1 (MQ=37)
                     aaCGCAATCGTTGATGGCGTTTGGTCGCGATGAACTGTGGGGACAAGGttt                                                              <  1:2628487/51‑1 (MQ=255)
                            tCGCTGATGGGGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGa                                      >  1:1587064/1‑68 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:3099613/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:2610414/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:1098543/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:1348281/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:588596/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:114181/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:2131759/1‑69 (MQ=255)
                            tCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAg                                     >  1:69436/1‑69 (MQ=255)
                              gCTGATGGCGTTTTGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:340490/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:1346005/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:3038392/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:2397804/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:3119967/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:2452635/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:2436548/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:2238389/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:2143900/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:978898/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:3287423/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:3385140/69‑1 (MQ=255)
                              gCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTa                                   <  1:1864786/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:1079741/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1802276/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1807511/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1906635/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1958668/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:3348674/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1553271/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:2997836/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:2201737/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1352081/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:2625173/1‑69 (MQ=255)
                                         ttGGTCGCGGTGAACTGTGGGGACAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1298716/1‑69 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGc                              >  1:549283/1‑62 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:71295/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1016419/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:491518/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)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1371736/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:2093404/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1536738/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1974028/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1788525/1‑68 (MQ=255)
                                          tGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAg                        >  1:1663652/1‑68 (MQ=255)
                                               gcggcgAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGca                    >  1:2928866/1‑67 (MQ=255)
                                                               gCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGc  <  1:1300438/69‑1 (MQ=255)
                                                               gCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGc  <  1:956975/69‑1 (MQ=255)
                                                               |                                                                    
CTTTGGCAGCGGCTATCAGTTAACGCAATCGCTGATGGCGTTTGGTCGCGGCGAACTTTGGGGGCAAGGTTTAGGTAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGC  >  minE/83913‑84044

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: