Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,677,414 A→T 14.9% R172R (CGT→CGA) ‡ yhgA ← predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,677,4140AT14.9% 221.5 / 28.5 87R172R (CGT→CGA) ‡yhgApredicted transposase
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base T (2/11);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAA  >  minE/2677345‑2677483
                                                                     |                                                                     
cATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATAc                                                                      <  1:4373978/71‑1 (MQ=255)
    tGCTCTGCCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATa                                                                       >  1:728882/1‑66 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                  >  1:4488858/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                  >  1:1750307/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                  >  1:1580988/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                  >  1:5485432/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGAGCCATACGTCg                                                                  >  1:4893237/1‑71 (MQ=255)
      ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                 <  1:422764/70‑1 (MQ=255)
      ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                 <  1:2703192/70‑1 (MQ=255)
      ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                 <  1:1545122/70‑1 (MQ=255)
       tctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                               >  1:140882/1‑71 (MQ=255)
           ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                            >  1:15778/1‑70 (MQ=255)
           ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                            >  1:2235050/1‑70 (MQ=255)
           ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                            >  1:2848944/1‑70 (MQ=255)
             aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                               >  1:4615940/1‑65 (MQ=255)
             aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                               >  1:2862998/1‑65 (MQ=255)
              atcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                                         <  1:2602926/68‑1 (MQ=255)
               tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa                                                       >  1:5350984/2‑69 (MQ=255)
                 ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                     <  1:1087263/71‑1 (MQ=255)
                 ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                     <  1:4987377/71‑1 (MQ=255)
                    tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATAc                                                                      <  1:55238/51‑1 (MQ=255)
                    tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCa                                                          >  1:4551637/1‑63 (MQ=255)
                    tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca                                                  >  1:2979993/1‑71 (MQ=255)
                     cAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca                                                  >  1:4320673/1‑70 (MQ=255)
                      aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                 >  1:3072811/1‑70 (MQ=255)
                      aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                 >  1:3514677/1‑70 (MQ=255)
                      aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                 >  1:1993096/1‑70 (MQ=255)
                        gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                                                                    <  1:4713278/49‑1 (MQ=255)
                        gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                              <  1:2660088/71‑1 (MQ=255)
                           gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                            <  1:4776253/70‑2 (MQ=255)
                           gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                            <  1:4977092/70‑2 (MQ=255)
                           gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                            <  1:1036944/70‑2 (MQ=255)
                           gcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                 >  1:5460968/1‑49 (MQ=255)
                            cgctGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCgg                                          <  1:5353350/68‑1 (MQ=255)
                              ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                                                                   <  1:886975/44‑1 (MQ=255)
                              ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCa                                        <  1:5202700/71‑1 (MQ=255)
                               tGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                       >  1:3925169/1‑71 (MQ=255)
                                gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                       >  1:1599291/1‑70 (MQ=255)
                                 aCGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                       <  1:921756/69‑1 (MQ=255)
                                 aCGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGa                                     <  1:5650889/71‑1 (MQ=255)
                                   gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACg                                   >  1:2702693/1‑71 (MQ=255)
                                   gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                        <  1:3992153/50‑2 (MQ=255)
                                   gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                        <  1:3023090/50‑2 (MQ=255)
                                   gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                        <  1:5350784/50‑2 (MQ=255)
                                    aaTATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                               <  1:2917339/42‑1 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                 <  1:2240814/71‑1 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                 >  1:5192344/1‑71 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                 >  1:665995/1‑71 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                 >  1:1906093/1‑71 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                 >  1:1891867/1‑71 (MQ=255)
                                      tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGa                                                       >  1:1334354/1‑48 (MQ=255)
                                      tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                 <  1:889489/70‑1 (MQ=255)
                                        tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                              <  1:389629/71‑1 (MQ=255)
                                        tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                              <  1:1823691/71‑1 (MQ=255)
                                        tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGAATGATTTCATCATCAGGCATGACGGTGAc                              <  1:2451865/71‑1 (MQ=255)
                                         gTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACa                             >  1:1055009/1‑71 (MQ=255)
                                          tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                           <  1:5547609/56‑1 (MQ=255)
                                          tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                           <  1:3642943/56‑1 (MQ=255)
                                           ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                            >  1:2630057/1‑70 (MQ=255)
                                           ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                            >  1:2246251/1‑70 (MQ=255)
                                            tttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatc                           <  1:3539959/70‑1 (MQ=255)
                                            tctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatc                                             >  1:336148/3‑52 (MQ=255)
                                            tctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCCGCACCACGGTAATatc                           <  1:5120915/68‑1 (MQ=255)
                                             ctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatcc                          <  1:1557292/69‑2 (MQ=255)
                                              tGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaat                        >  1:2548112/1‑71 (MQ=255)
                                               gAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatc                           <  1:2017244/67‑1 (MQ=255)
                                                aaTTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaatca                      <  1:5275029/71‑1 (MQ=255)
                                                   tAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                              >  1:3982781/1‑44 (MQ=255)
                                                    aaCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCg                   >  1:4759618/1‑70 (MQ=255)
                                                     aCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                 <  1:2503480/71‑1 (MQ=255)
                                                       tCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaa                >  1:3936979/1‑70 (MQ=255)
                                                        ccAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAAt              >  1:3397653/1‑71 (MQ=255)
                                                         cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaaa               >  1:4642543/1‑69 (MQ=255)
                                                         cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg             >  1:3098297/1‑71 (MQ=255)
                                                         cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg             >  1:391025/1‑71 (MQ=255)
                                                          agcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg             >  1:4353244/1‑70 (MQ=255)
                                                            cagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCgg          <  1:4105215/71‑1 (MQ=255)
                                                            cagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg           >  1:4473905/1‑70 (MQ=255)
                                                             agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa         <  1:494097/71‑1 (MQ=255)
                                                             agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa         <  1:3949063/71‑1 (MQ=255)
                                                             agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa         <  1:4718488/71‑1 (MQ=255)
                                                             agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa         <  1:3833914/71‑1 (MQ=255)
                                                               cgcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAgg       <  1:4749342/71‑1 (MQ=255)
                                                                gcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGc      <  1:4308832/71‑1 (MQ=255)
                                                                   aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa   <  1:5114754/71‑1 (MQ=255)
                                                                   aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa   <  1:1096705/71‑1 (MQ=255)
                                                                    tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa  <  1:5646981/71‑1 (MQ=255)
                                                                    tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa  <  1:1019130/71‑1 (MQ=255)
                                                                     |                                                                     
CATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAA  >  minE/2677345‑2677483

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