Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 89,548 A→G 13.6% E276G (GAA→GGA)  ftsA → ATP‑binding cell division protein involved in recruitment of FtsK to Z ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE89,5480AG13.6% 218.1 / 23.4 81E276G (GAA→GGA) ftsAATP‑binding cell division protein involved in recruitment of FtsK to Z ring
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base G (4/7);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGCC  >  minE/89480‑89617
                                                                    |                                                                     
cGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAg                                                                     <  1:620400/71‑1 (MQ=255)
 gCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGt                                                                    >  1:3732611/1‑71 (MQ=255)
 gCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGt                                                                    >  1:165785/1‑71 (MQ=255)
   cGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGc                                                                  >  1:286375/1‑71 (MQ=255)
   cGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGc                                                                  >  1:2053555/1‑71 (MQ=255)
   cGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGc                                                                  >  1:3391404/1‑71 (MQ=255)
   cGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGc                                                                  >  1:2939742/1‑71 (MQ=255)
    gAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGcc                                                                 <  1:1684126/71‑1 (MQ=255)
    gAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGcc                                                                 <  1:3604008/71‑1 (MQ=255)
     aaGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGcc                                                                 <  1:1609403/70‑1 (MQ=255)
            tAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAg                                                         >  1:3389343/1‑71 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATTAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     <  1:3554565/71‑1 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     >  1:3428191/1‑71 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     <  1:2739196/71‑1 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     >  1:2203804/1‑71 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     <  1:1874898/71‑1 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     <  1:1200900/71‑1 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     <  1:3645248/71‑1 (MQ=255)
                gTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAggtgg                                                     <  1:778049/71‑1 (MQ=255)
                 ttCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAggtggt                                                    >  1:977371/1‑71 (MQ=255)
                 ttCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAggtgg                                                     >  1:1951577/1‑70 (MQ=255)
                  tCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGa                                                                       <  1:1376868/51‑1 (MQ=255)
                   cGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGgtc                                                   <  1:1127294/70‑1 (MQ=255)
                       aCGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCg                                              >  1:3778896/1‑71 (MQ=255)
                         ggTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGcc                                            <  1:918048/71‑1 (MQ=255)
                          gTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCa                                           <  1:1426665/71‑1 (MQ=255)
                             gtgCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:2923219/70‑1 (MQ=255)
                             gtgCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:443590/70‑1 (MQ=255)
                              tgCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:2062765/69‑1 (MQ=255)
                              tgCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:1382498/69‑1 (MQ=255)
                               gcgcTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAg                                      >  1:1290186/1‑71 (MQ=255)
                               gcgcTGGGTTCCATCGTTGGAAAAGATGAGAGCGAGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAg                                      >  1:55722/1‑71 (MQ=255)
                                 gcTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTc                                    >  1:2611066/1‑71 (MQ=255)
                                  cTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCt                                   >  1:1185074/1‑71 (MQ=255)
                                   tGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTg                                  >  1:1621455/1‑71 (MQ=255)
                                     ggTTCCATCGTTGGAAAAGATGAGAGCGTGGGAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGc                                 <  1:2539353/70‑1 (MQ=255)
                                     ggTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTg                                  <  1:2406200/69‑1 (MQ=255)
                                     ggTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGc                                 >  1:473422/1‑70 (MQ=255)
                                     ggTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGc                                 <  1:572187/70‑1 (MQ=255)
                                       ttCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAAc                              >  1:2465676/1‑71 (MQ=255)
                                         ccATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACg                             >  1:226678/1‑70 (MQ=255)
                                         ccATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGt                            >  1:378295/1‑71 (MQ=255)
                                         ccATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGt                            >  1:3020789/1‑71 (MQ=255)
                                         ccATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGt                            >  1:565190/1‑71 (MQ=255)
                                         ccATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGt                            >  1:1606269/1‑71 (MQ=255)
                                           aTCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCa                          <  1:630886/71‑1 (MQ=255)
                                              gTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:958333/53‑1 (MQ=255)
                                              gTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:352544/53‑1 (MQ=255)
                                                tGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGaca                      >  1:424915/1‑70 (MQ=255)
                                                 ggAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACt                    >  1:740274/1‑71 (MQ=255)
                                                   aaaaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:2244737/48‑1 (MQ=255)
                                                   aaaaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:1200342/48‑1 (MQ=255)
                                                   aaaaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:839989/48‑1 (MQ=255)
                                                   aaaaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:1281775/48‑1 (MQ=255)
                                                   aaaaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:3193690/48‑1 (MQ=255)
                                                     aaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTgg                  >  1:1280136/1‑69 (MQ=255)
                                                     aaGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTgg                  >  1:3339364/1‑69 (MQ=255)
                                                        aTGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACt                    <  1:3550031/64‑1 (MQ=255)
                                                        aTGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACt                    <  1:3578826/64‑1 (MQ=255)
                                                         tGAGAGCGTGGAAGTTCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAgg            <  1:3623093/71‑1 (MQ=255)
                                                          gagaGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACg                             >  1:1570650/1‑53 (MQ=255)
                                                          gagaGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACg                             >  1:1359800/1‑53 (MQ=255)
                                                          gagaGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACg                             >  1:758341/1‑53 (MQ=255)
                                                          gagaGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACg                             >  1:2361395/1‑53 (MQ=255)
                                                           agagCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGt           <  1:1630852/70‑1 (MQ=255)
                                                           agagCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTg          <  1:2341679/71‑1 (MQ=255)
                                                            gagCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:2299682/39‑1 (MQ=255)
                                                            gagCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:2585489/39‑1 (MQ=255)
                                                             agCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACg                                         <  1:2090343/38‑1 (MQ=255)
                                                              gCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATc       >  1:2734022/1‑71 (MQ=255)
                                                               cGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCg      <  1:1057839/71‑1 (MQ=255)
                                                                 tGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGt           <  1:640310/64‑1 (MQ=255)
                                                                 tGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGt           <  1:1234247/64‑1 (MQ=255)
                                                                 tGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGt           <  1:3559654/64‑1 (MQ=255)
                                                                 tGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGt           <  1:2203790/64‑1 (MQ=255)
                                                                 tGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGt           <  1:1871913/64‑1 (MQ=255)
                                                                   gAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGa         >  1:3503731/1‑64 (MQ=255)
                                                                   gAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGcc  >  1:3818075/1‑71 (MQ=255)
                                                                   gAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGcc  >  1:1109212/1‑71 (MQ=255)
                                                                   gAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGcc  >  1:1309727/1‑71 (MQ=255)
                                                                   gAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGc   <  1:482324/70‑1 (MQ=255)
                                                                   gAAGTGACGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGcc  >  1:525980/1‑71 (MQ=255)
                                                                    |                                                                     
CGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGCAACGTCAGACACTGGCAGAGGTGATCGAGCC  >  minE/89480‑89617

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: