Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,842,330 T→C 37.3% Q52Q (CAA→CAG)  ftsB ← cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,842,3300TC37.3% 63.2 / 70.8 111Q52Q (CAA→CAG) ftsBcell division protein
Reads supporting (aligned to +/- strand):  ref base T (41/28);  new base C (24/17);  total (66/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 112 total aligned reads displayed.

TTCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCGGG  >  minE/1842264‑1842397
                                                                  |                                                                   
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGtt                                                                    >  1:913590/1‑68 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:2594825/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:3147786/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   <  1:3171911/69‑1 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:1066858/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:386806/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:1591389/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:1488763/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   <  1:624750/69‑1 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGATg                                                                   >  1:912429/1‑69 (MQ=255)
ttCATTACGCGCACGCTCTTCGAGCGCCCCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   >  1:5102/1‑69 (MQ=255)
 tCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   <  1:2343801/68‑1 (MQ=255)
 tCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTg                                                                   <  1:503614/68‑1 (MQ=255)
  cATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGAt                                                                 >  1:1040916/1‑69 (MQ=255)
  cATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGAt                                                                 >  1:116624/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGTCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:1383549/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGtt                                                             <  1:960011/69‑1 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGtt                                                             >  1:3010968/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGtt                                                             <  1:1702724/69‑1 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGtt                                                             >  1:785213/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:3169048/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:1244930/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:2594671/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:2257472/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:563883/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:690629/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:63826/1‑69 (MQ=255)
      aCGCGCACGCTCTTCGAGCGCCTACTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGtt                                                             >  1:2587400/1‑69 (MQ=255)
           cACGCTCTTCGCGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTc                                                           >  1:699503/1‑66 (MQ=255)
           cACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTcgcg                                                        >  1:42099/1‑69 (MQ=255)
               ctctTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCt                                                      >  1:1730947/1‑67 (MQ=255)
               ctctTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCt                                                      >  1:950265/1‑67 (MQ=255)
               ctctTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCt                                                      >  1:2676006/1‑67 (MQ=255)
                    cGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGtt                                               >  1:297299/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTcg                                                          >  1:224824/1‑52 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTcgcg                                          >  1:2190254/1‑68 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTcgcg                                          >  1:755296/1‑68 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2280222/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2527238/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2815748/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2514245/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2741729/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2581812/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:2023733/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:349557/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:411029/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:149718/1‑69 (MQ=255)
                          cctcctGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGt                                         >  1:1495782/1‑69 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:2577073/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:3166133/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:3170072/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:943561/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:3316458/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:556573/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:563809/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:1038599/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:1232025/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:1375414/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:1496268/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:2392275/68‑1 (MQ=255)
                                   cgccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:1924929/68‑1 (MQ=255)
                                    gccATTGAGATCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:2183778/67‑1 (MQ=255)
                                        ttGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCtgt                           <  1:1540107/69‑1 (MQ=255)
                                        ttGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCtgt                           <  1:2417544/69‑1 (MQ=255)
                                          gataTCGTCAATTTCGGAAAAAAGTTGATCGTTTCGAGCTTTAAGTTTCGCGTTTGTAGGTTGCtgt                           <  1:2699076/64‑1 (MQ=255)
                                          gagaTCGTCAATTTCGTCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:666971/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:989179/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:1473730/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:1663764/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:1817403/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:3317038/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTcgcg                                          <  1:2679435/52‑1 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgc                         >  1:204884/1‑69 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgc                         >  1:267985/1‑69 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgc                         >  1:2692674/1‑69 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgc                         >  1:1399471/1‑69 (MQ=255)
                                          gagaTCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgc                         >  1:1495181/1‑69 (MQ=255)
                                           agaTCGTCAATTTCGGCAAAAAGCTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCt                                 <  1:2856976/60‑1 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:2674778/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:740806/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:1320654/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:245594/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:1531047/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:2142214/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:1734422/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:1816217/1‑67 (MQ=255)
                                               cGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTgccgc                      >  1:2242425/1‑67 (MQ=255)
                                                gTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCtgt                           >  1:817774/1‑61 (MQ=255)
                                                gTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCcac                   <  1:503050/69‑1 (MQ=255)
                                                     tttCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGt                                     <  1:2247244/46‑1 (MQ=255)
                                                     tttCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGt                                     <  1:1066081/46‑1 (MQ=255)
                                                     tttCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGt                                     <  1:891679/46‑1 (MQ=255)
                                                     tttCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGt                                     <  1:254727/46‑1 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGGCGCCACATCATCATTGACGCggg  >  1:1251430/1‑69 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATGATCATTGACGCggg  >  1:1150741/1‑69 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCggg  >  1:2208890/1‑69 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCggg  >  1:1148138/1‑69 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCggg  >  1:105113/1‑69 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCggg  <  1:1045524/69‑1 (MQ=255)
                                                                 gTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCgg   >  1:416935/1‑68 (MQ=255)
                                                                  |                                                                   
TTCATTACGCGCACGCTCTTCGAGCGCCTCCTGGCCGCCATTGAGATCGTCAATTTCGGCAAAAAGTTGATCGTTTCGCGCTTTAAGTTTCGCGTTTGTAGCTTGCTGTGCCGCCACATCATCATTGACGCGGG  >  minE/1842264‑1842397

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: ‑