Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 90,611 G→A 41.6% L189L (CTG→CTA ftsZ → GTP‑binding tubulin‑like cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE90,6110GA41.6% 43.9 / 78.1 108L189L (CTG→CTAftsZGTP‑binding tubulin‑like cell division protein
Reads supporting (aligned to +/- strand):  ref base G (40/23);  new base A (31/14);  total (71/37)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 110 total aligned reads displayed.

AACTGCTGAAAGTTCTGGGCCGCGGTATCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGA  >  minE/90544‑90679
                                                                   |                                                                    
aaCTGCTGAAAGTTCTGGGCAGCGGTATCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGa                                                                     <  1:2897266/69‑1 (MQ=255)
             tCTGGGCCTCGGTATCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCa                                                        <  1:345788/69‑1 (MQ=255)
             tCTGGGCCGCGGTATCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCa                                                        <  1:2445858/69‑1 (MQ=255)
                           tctcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGaa                                           >  1:3420211/1‑68 (MQ=255)
                             tcactgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAAc                                          >  1:1252950/4‑67 (MQ=255)
                             tcCCTGCTGGATGCTTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2914477/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAGCGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:1968314/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGaa                                           >  1:1172040/1‑66 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACt                                         >  1:72648/1‑68 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACt                                         >  1:1713344/1‑68 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACt                                         >  1:1947268/1‑68 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACt                                         >  1:1318696/1‑68 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:3195691/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2037005/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2401500/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2135346/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2114648/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2046722/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2012996/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2621638/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:3039121/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2627451/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2853112/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:2895311/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1862982/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1829650/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:678187/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1324651/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:551994/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1292041/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:911289/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1201638/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1402729/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:95353/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1623049/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                        >  1:1851884/1‑69 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:2951743/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:127734/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:308942/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:3144389/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:3183360/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:971688/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:3205864/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:2303115/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:2064840/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:3241054/68‑1 (MQ=255)
                             tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:523370/68‑1 (MQ=255)
                             tcCCAGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                         <  1:3204253/68‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:1903064/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2380702/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:3395721/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:922328/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:372729/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2674328/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:3290774/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:257119/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2269012/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2474928/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:732512/67‑1 (MQ=255)
                                 tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2380847/67‑1 (MQ=255)
                                  gctgTATGCGTTTGGCGCATCGAACGATTTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2472157/66‑1 (MQ=255)
                                         gCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCg                            >  1:256327/1‑69 (MQ=255)
                                         gCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCg                            >  1:1320325/1‑69 (MQ=255)
                                         gCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCg                            >  1:2348390/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatga                  <  1:2797491/69‑1 (MQ=255)
                                                   cAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatga                  <  1:1147882/69‑1 (MQ=255)
                                                   cAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatga                  <  1:1355734/69‑1 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTatgatg                   >  1:1368363/1‑68 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:559510/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1591915/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1589846/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:146815/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1400184/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:379873/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:410737/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3304389/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:568805/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1353335/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1319896/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1287913/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1191638/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:1135229/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:2896474/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:2913644/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:2622488/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3052854/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3054147/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:2617934/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3140890/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3155773/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:24696/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:235902/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:2198077/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3263010/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:2094147/1‑69 (MQ=255)
                                                   cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                  >  1:3316830/1‑69 (MQ=255)
                                                      cGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:1069289/46‑1 (MQ=255)
                                                          cGATGTACTGAAATGCGCTGTGCAAGGTATCGCTGAACTGAt                                      <  1:2533245/42‑1 (MQ=255)
                                                                   gAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGa  <  1:2963860/69‑1 (MQ=255)
                                                                   gAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGa  <  1:2746582/69‑1 (MQ=255)
                                                                   |                                                                    
AACTGCTGAAAGTTCTGGGCCGCGGTATCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGA  >  minE/90544‑90679

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: