Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,677,417 T→G 6.5% R171R (CGA→CGC yhgA ← predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,677,4170TG 6.5% 323.4 / 18.3 124R171R (CGA→CGCyhgApredicted transposase
Reads supporting (aligned to +/- strand):  ref base T (50/66);  new base G (2/6);  total (52/72)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 125 total aligned reads displayed.

TTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGC  >  minE/2677348‑2677486
                                                                     |                                                                     
ttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                                                                      >  1:2458756/1‑71 (MQ=255)
 tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                                                                      <  1:780550/70‑1 (MQ=255)
 tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                                                                      <  1:785570/70‑1 (MQ=255)
 tGCTCTACAAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                     >  1:2220854/1‑71 (MQ=255)
  gCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                     <  1:1103915/70‑1 (MQ=255)
   ctctACCAGCCCCATCAGGTCGCGTTTACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGt                                                                   <  1:1936557/71‑1 (MQ=255)
          aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                                                            >  1:213696/1‑71 (MQ=255)
          aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                                                            <  1:2273773/71‑1 (MQ=255)
          aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                                                            <  1:3037099/71‑1 (MQ=255)
           gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                                                            >  1:686605/1‑70 (MQ=255)
           atcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                                            <  1:2748717/68‑1 (MQ=255)
             cccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAt                                                         <  1:2768629/71‑1 (MQ=255)
                aacAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtctc                                                      <  1:867482/69‑1 (MQ=255)
                aTCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                     >  1:718769/1‑56 (MQ=255)
                 tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                    >  1:1112838/1‑56 (MQ=255)
                 tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                     <  1:405496/55‑1 (MQ=255)
                 tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg                                                           <  1:2665079/65‑1 (MQ=255)
                   aGGTCGCGTTGATGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                        <  1:1152073/66‑1 (MQ=255)
                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                                                                       >  1:2140804/1‑51 (MQ=255)
                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                                                            <  1:2743089/62‑1 (MQ=255)
                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                    >  1:2628654/1‑70 (MQ=255)
                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                    >  1:3342403/1‑70 (MQ=255)
                     atCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                           <  1:1673469/60‑2 (MQ=255)
                      tCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcatc                                                <  1:2789595/71‑1 (MQ=255)
                       cgcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                 >  1:2633924/1‑69 (MQ=255)
                       cgcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                 >  1:1064608/1‑69 (MQ=255)
                       cgcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                 >  1:3066812/1‑69 (MQ=255)
                       cgcgCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                               >  1:3869491/1‑70 (MQ=255)
                        gcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcatca                                               <  1:643106/70‑1 (MQ=255)
                        gcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                              >  1:3474994/1‑71 (MQ=255)
                         cgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAgg                                             <  1:213280/71‑1 (MQ=255)
                           ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCa                                           >  1:579528/1‑71 (MQ=255)
                            tGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                          >  1:3965443/1‑71 (MQ=255)
                             gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                          >  1:2668514/1‑70 (MQ=255)
                             gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                          <  1:2432739/70‑1 (MQ=255)
                             gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                          <  1:2946684/70‑1 (MQ=255)
                             gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                          <  1:3796741/70‑1 (MQ=255)
                               cGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGAc                                       <  1:2133907/71‑1 (MQ=255)
                                gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCa                                           >  1:515905/1‑66 (MQ=255)
                                gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGAc                                       <  1:3300100/70‑1 (MQ=255)
                                gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGAc                                       <  1:2928519/70‑1 (MQ=255)
                                gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGAc                                       <  1:743027/70‑1 (MQ=255)
                                 aaTATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGAc                                       <  1:1756823/69‑1 (MQ=255)
                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCa                                                             >  1:1036151/1‑46 (MQ=255)
                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                    <  1:226263/71‑1 (MQ=255)
                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                    >  1:949823/1‑71 (MQ=255)
                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                    <  1:1742926/71‑1 (MQ=255)
                                  atatGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGt                                    <  1:972013/71‑1 (MQ=255)
                                   tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTg                                   <  1:3971231/71‑1 (MQ=255)
                                     tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGa                                        >  1:1157143/1‑64 (MQ=255)
                                     tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGa                                  >  1:4068641/1‑70 (MQ=255)
                                      gTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGa                                  >  1:2343386/1‑69 (MQ=255)
                                       tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg                                                           <  1:133888/43‑1 (MQ=255)
                                       tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcatca                                               >  1:96637/1‑55 (MQ=255)
                                       tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                    <  1:492058/66‑1 (MQ=255)
                                        ttctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatc                              >  1:2935972/4‑71 (MQ=255)
                                         tttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAgg                                             <  1:3571142/55‑1 (MQ=255)
                                         tttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatc                              <  1:3514089/70‑1 (MQ=255)
                                         tttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatc                              <  1:2776165/70‑1 (MQ=255)
                                          ttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCa                                           >  1:2054602/1‑56 (MQ=255)
                                          ttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatca                             <  1:1885155/70‑1 (MQ=255)
                                              aTTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAg                        >  1:3741100/1‑71 (MQ=255)
                                              aTTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAg                        <  1:2642815/71‑1 (MQ=255)
                                               ttAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAgg                                             <  1:305933/49‑1 (MQ=255)
                                               ttAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGc                       >  1:1876252/1‑71 (MQ=255)
                                                  aCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                    <  1:1013238/71‑1 (MQ=255)
                                                  aCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                    <  1:1950184/71‑1 (MQ=255)
                                                  aCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                    <  1:1965513/71‑1 (MQ=255)
                                                    tCCAGCAGCGCCATACGTCGGTGCTGCATGATTtc                                                      <  1:2435502/35‑1 (MQ=255)
                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                    <  1:92304/51‑1 (MQ=255)
                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                    <  1:642744/51‑1 (MQ=255)
                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                    <  1:886151/67‑1 (MQ=255)
                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                    <  1:940192/67‑1 (MQ=255)
                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAAt                 <  1:1429735/70‑1 (MQ=255)
                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAAt                 <  1:1845628/70‑1 (MQ=255)
                                                       agcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGc               <  1:13675/71‑1 (MQ=255)
                                                        gcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg              <  1:2711612/71‑1 (MQ=255)
                                                        gcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg              <  1:1155648/71‑1 (MQ=255)
                                                        gcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg              >  1:1354640/1‑71 (MQ=255)
                                                        gcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg              >  1:591913/1‑71 (MQ=255)
                                                        gcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg              >  1:2598025/1‑71 (MQ=255)
                                                           gcgcCATACGTCGGTGCTGCATGATTTAATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAg           >  1:3446834/1‑71 (MQ=255)
                                                            cgcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAgg          <  1:1393371/71‑1 (MQ=255)
                                                             gcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAgg          >  1:3236335/1‑70 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaatca                         <  1:1130800/52‑1 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaatca                         <  1:109437/52‑1 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCATTCAGCGGAAATGCGGAGGCATa      <  1:3300270/71‑1 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCAt       >  1:3073373/1‑70 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCAt       >  1:1469550/1‑70 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCAt       >  1:3647826/1‑70 (MQ=255)
                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa      <  1:873938/71‑1 (MQ=255)
                                                                 tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAg           <  1:757004/65‑1 (MQ=255)
                                                                  aCGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaaa                  >  1:794964/1‑57 (MQ=255)
                                                                  aCGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGc               <  1:3208772/60‑1 (MQ=255)
                                                                  aCGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg              <  1:465720/61‑1 (MQ=255)
                                                                  aCGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCa        <  1:23824/67‑1 (MQ=255)
                                                                   cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaat                           >  1:1810413/1‑47 (MQ=255)
                                                                   cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGc         >  1:424480/1‑65 (MQ=255)
                                                                   cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg   >  1:1005280/1‑71 (MQ=255)
                                                                    gTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGc  >  1:654495/1‑71 (MQ=255)
                                                                     |                                                                     
TTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGC  >  minE/2677348‑2677486

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: