Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 90,635 T→C 36.3% A197A (GCT→GCC 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,6350TC36.3% 105.1 / 105.6 165A197A (GCT→GCCftsZGTP‑binding tubulin‑like cell division protein
Reads supporting (aligned to +/- strand):  ref base T (53/52);  new base C (46/14);  total (99/66)
Fisher's exact test for biased strand distribution p-value = 9.79e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 165 total aligned reads displayed.

TCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGA  >  minE/90571‑90700
                                                                |                                                                 
  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:1713344/1‑68 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACt                                                              >  1:72648/1‑68 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:1324651/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:1829650/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:1862982/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:1623049/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2012996/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:3195691/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2046722/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:1402729/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2114648/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:546840/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2135346/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:551994/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:1292227/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2895311/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:95353/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2853112/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTg                                                             >  1:2627451/1‑69 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                                              <  1:971688/68‑1 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                                              <  1:3183360/68‑1 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                                              <  1:2064840/68‑1 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                                              <  1:3144389/68‑1 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                                              <  1:3205864/68‑1 (MQ=255)
  tcCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACt                                                              <  1:308942/68‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:3290774/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:922328/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:2380702/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:2674328/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:257119/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:372729/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:1565018/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:3395721/67‑1 (MQ=255)
      tgctgGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:1903064/67‑1 (MQ=255)
       gctgTATGCGTTTGGCGCATCGAACGATTTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGAt                                                           <  1:2472157/66‑1 (MQ=255)
              gCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACt                                                       >  1:3066910/1‑63 (MQ=255)
              gCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCg                                                 >  1:256327/1‑69 (MQ=255)
                        cAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatga                                       <  1:2797491/69‑1 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTatgatg                                        >  1:1368363/1‑68 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2617934/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:24696/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:3263010/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1135229/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2622488/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2316404/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2896474/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1191638/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:3052854/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:568805/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1287913/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1319896/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1353335/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:379873/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:146815/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1589846/1‑69 (MQ=255)
                        cAGCGAACGATGTACTAAAAGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:3304389/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)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2195509/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2695980/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1202016/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:720212/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:716324/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1306073/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:509697/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:229800/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:1468495/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:159332/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGGTTtgatga                                       >  1:2227727/1‑50 (MQ=255)
                                           aGGCGCTGTGCAAGGTATCGCCGAACTGATTACTCGTCCGGCTTtgatga                                       >  1:2042523/1‑50 (MQ=255)
                                               gcTGTGCAAGTTATCGCTGAACTGATTACTCGTCCGGGTTtgatg                                        <  1:2189773/45‑1 (MQ=255)
                                               gcTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatg                                        <  1:1059188/45‑1 (MQ=255)
                                               gcTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatg                                        <  1:1986832/45‑1 (MQ=255)
                                               gcTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatg                                        <  1:3252981/45‑1 (MQ=255)
                                               gcTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTtgatg                                        <  1:2262891/45‑1 (MQ=255)
                                                cTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCa               <  1:1191097/69‑1 (MQ=255)
                                                cTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCa               <  1:836157/69‑1 (MQ=255)
                                                     cAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTa          >  1:236883/1‑69 (MQ=255)
                                                     cAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTa          >  1:2350321/1‑69 (MQ=255)
                                                           aTCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTg                           >  1:2092826/1‑46 (MQ=255)
                                                           aTCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTg                           <  1:1814658/46‑1 (MQ=255)
                                                           aTCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTg                           >  1:948565/1‑46 (MQ=255)
                                                            tCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGt      >  1:2394031/1‑66 (MQ=255)
                                                            tCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGt      >  1:1147577/1‑66 (MQ=255)
                                                            tCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGt      >  1:1252267/1‑66 (MQ=255)
                                                            tCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGt      >  1:1325434/1‑66 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:2318799/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:2789990/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:2741195/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:781639/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:1396535/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:2986345/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:1700916/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:1721909/68‑1 (MQ=255)
                                                              gCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTga  <  1:2242756/68‑1 (MQ=255)
                                                                |                                                                 
TCTCCCTGCTGGATGCGTTTGGCGCAGCGAACGATGTACTGAAAGGCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGA  >  minE/90571‑90700

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