Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,677,204 A→T 8.5% S242R (AGT→AGA yhgA ← predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,677,2040AT 8.5% 294.1 / 15.7 106S242R (AGT→AGAyhgApredicted transposase
Reads supporting (aligned to +/- strand):  ref base A (48/49);  new base T (2/7);  total (50/56)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 120 total aligned reads displayed.

ATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGC  >  minE/2677133‑2677271
                                                                       |                                                                   
attattTTGGCTATATGCAGGGCTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                                     >  1:1993270/1‑71 (MQ=255)
 ttattTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                                     <  1:3674607/71‑1 (MQ=255)
 ttattTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                                     <  1:2999804/71‑1 (MQ=255)
   attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAAct                                                                   >  1:1442851/1‑71 (MQ=255)
   attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAAc                                                                    >  1:1614431/1‑70 (MQ=255)
      ttGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAActct                                                                 <  1:4030237/70‑1 (MQ=255)
         gCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCtt                                                             <  1:2238824/71‑1 (MQ=255)
            atatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTtgtg                                                          >  1:3018182/1‑71 (MQ=255)
             tatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTgtgt                                                         >  1:2640550/1‑71 (MQ=255)
               tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAActc                                                                  <  1:2322354/60‑1 (MQ=255)
               tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                                       >  1:5674497/1‑71 (MQ=255)
               tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                                       >  1:458701/1‑71 (MQ=255)
               tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                                       >  1:2039009/1‑71 (MQ=255)
                gCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTc                                                      >  1:2855571/1‑71 (MQ=255)
                 cAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCg                                                     >  1:2114348/1‑71 (MQ=255)
                  aGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCgg                                                    <  1:3123443/71‑1 (MQ=255)
                   gggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                                       <  1:4605188/67‑1 (MQ=255)
                    ggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGt                                                   >  1:2284974/1‑70 (MQ=255)
                      cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGTAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                                                >  1:665471/1‑71 (MQ=255)
                      cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTg                                                  <  1:2654664/69‑1 (MQ=255)
                      cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                                                >  1:5566757/1‑71 (MQ=255)
                        ttGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                                                >  1:1695325/1‑69 (MQ=255)
                         tGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                                             >  1:2339278/1‑71 (MQ=255)
                         tGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                                             >  1:3608937/1‑71 (MQ=255)
                         tGGATTGTTCCACCTCCTGCCGCAATATTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                                             >  1:3721111/1‑71 (MQ=255)
                           gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGtt                                            <  1:2158191/70‑1 (MQ=255)
                           gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTc                                           >  1:1466136/1‑71 (MQ=255)
                             ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCgg                                         <  1:2933328/71‑1 (MQ=255)
                             ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCgg                                         <  1:3575475/71‑1 (MQ=255)
                              tGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                                     >  1:1461287/1‑42 (MQ=255)
                               gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                                        <  1:1473989/70‑1 (MQ=255)
                               gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                                        <  1:2483895/70‑1 (MQ=255)
                                ttCCCACTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCa                                       <  1:1604583/70‑1 (MQ=255)
                                 tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                                     >  1:2550433/1‑39 (MQ=255)
                                 tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCg                                          >  1:783954/1‑66 (MQ=255)
                                 tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCg                                          >  1:4082537/1‑66 (MQ=255)
                                 tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      >  1:5549790/1‑70 (MQ=255)
                                 tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                                     >  1:2568676/1‑71 (MQ=255)
                                 tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                                     >  1:3991758/1‑71 (MQ=255)
                                   ccccTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACg                                    <  1:1969183/70‑1 (MQ=255)
                                      ttcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgta                                                        >  1:3008479/2‑46 (MQ=12)
                                      ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      <  1:3778954/65‑1 (MQ=255)
                                      ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      <  1:4730658/65‑1 (MQ=255)
                                       tcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgt                                                         <  1:5020024/45‑1 (MQ=12)
                                       tcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgt                                                         <  1:1771849/45‑1 (MQ=12)
                                       tcctGCCGCAGCCGTTCAGCAATAGTCATTAATCTCTCCt                                                              >  1:2142348/1‑40 (MQ=16)
                                           gCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGAt                            >  1:5014097/1‑70 (MQ=255)
                                            ccGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTgtgt                                                         <  1:1110434/40‑1 (MQ=255)
                                              gCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaaa                        >  1:1918759/1‑71 (MQ=255)
                                                aaTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATc                      <  1:1136212/71‑1 (MQ=255)
                                                 aTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCg                     >  1:1451415/1‑71 (MQ=255)
                                                 aTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCg                     >  1:4838314/1‑71 (MQ=255)
                                                  tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTg                                                  <  1:1744585/41‑1 (MQ=255)
                                                  tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaaa                        <  1:4637351/67‑1 (MQ=255)
                                                  tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGt                    >  1:799715/1‑71 (MQ=255)
                                                      tCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      >  1:1179428/1‑49 (MQ=255)
                                                      tCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      >  1:4531809/1‑49 (MQ=255)
                                                       ccGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTaaa                <  1:181053/70‑1 (MQ=255)
                                                       ccGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTaaa                <  1:2082963/70‑1 (MQ=255)
                                                       ccGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTaaa                <  1:2252254/70‑1 (MQ=255)
                                                        cGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGTCAACGCCGTCGATAAAATCGTTa                  <  1:4852011/67‑1 (MQ=255)
                                                        cGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTa                  <  1:2641821/67‑1 (MQ=255)
                                                         gCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAAcc             <  1:3684715/71‑2 (MQ=255)
                                                         gCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACg             <  1:2738195/71‑1 (MQ=255)
                                                           aaTAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                                             <  1:2020595/37‑1 (MQ=255)
                                                            aTAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTa           >  1:3205731/1‑70 (MQ=255)
                                                            aTAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTAc          <  1:4614381/71‑1 (MQ=255)
                                                            aTAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTAc          <  1:3259337/71‑1 (MQ=255)
                                                             tAGTCATTAATCTCTCCTTTTGTAGTGGTGAACGTTCGGCa                                       <  1:1756549/41‑1 (MQ=19)
                                                             tAGTCATTAATCTCTCCTTTTGTAGTGGTGAACGTTCGGCa                                       <  1:4071968/41‑1 (MQ=19)
                                                             tAGTCATTAATCTCTCCTTTTGTAGTGGTGAACGTTCGGCa                                       <  1:475597/41‑1 (MQ=19)
                                                             tAGTCATTAATCTCTCCTTTTGTAGTGGTGAACGTTCGGCa                                       <  1:4997389/41‑1 (MQ=19)
                                                             tAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                                        <  1:4194965/40‑1 (MQ=255)
                                                             tAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                                        <  1:378502/40‑1 (MQ=255)
                                                             tAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                                        <  1:2912889/40‑1 (MQ=255)
                                                             tAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTa                  <  1:5292242/62‑1 (MQ=255)
                                                             tAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTa                  <  1:5620757/62‑1 (MQ=255)
                                                              aGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      >  1:1310955/1‑41 (MQ=255)
                                                              aGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                                      >  1:1619937/1‑41 (MQ=255)
                                                               gTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTc                                           >  1:4938327/1‑35 (MQ=255)
                                                               gTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTc                                           >  1:368193/1‑35 (MQ=255)
                                                               gTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTaaa                >  1:2873567/1‑62 (MQ=255)
                                                                  aTTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTAc          >  1:1352341/1‑65 (MQ=255)
                                                                  aTTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGt     <  1:568789/70‑1 (MQ=255)
                                                                   ttAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                                     <  1:2176697/37‑1 (MQ=255)
                                                                   ttAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaa                         <  1:4230020/49‑1 (MQ=255)
                                                                   ttAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaa                         <  1:1524692/49‑1 (MQ=255)
                                                                   ttAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaa                          <  1:4219297/48‑1 (MQ=255)
                                                                   ttAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTaa                 <  1:3595419/57‑1 (MQ=255)
                                                                    tAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGt                    <  1:4097980/53‑1 (MQ=255)
                                                                    tAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGc  >  1:4361196/1‑71 (MQ=255)
                                                                    tAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGc  >  1:311675/1‑71 (MQ=255)
                                                                    tAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGc  >  1:3939126/1‑71 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAAc              >  1:3317998/1‑58 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACa         >  1:5430790/1‑63 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAg        >  1:3799282/1‑64 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAg        >  1:3792578/1‑64 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGc  >  1:4483713/1‑70 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGc  >  1:718666/1‑70 (MQ=255)
                                                                     aaaCTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGc  >  1:775808/1‑70 (MQ=255)
                                                                       |                                                                   
ATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGTTAAAACGTACAGCGTCGC  >  minE/2677133‑2677271

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