Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,677,415 C→T 15.3% R172H (CGT→CAT) ‡ yhgA ← predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,677,4150CT15.3% 201.3 / 23.9 85R172H (CGT→CAT) ‡yhgApredicted transposase
Reads supporting (aligned to +/- strand):  ref base C (39/33);  new base T (2/11);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAG  >  minE/2677345‑2677485
                                                                      |                                                                      
cATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATAc                                                                        <  1:4373978/71‑1 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:5485432/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:1580988/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:1750307/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:4488858/1‑71 (MQ=255)
    tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGAGCCATACGTCg                                                                    >  1:4893237/1‑71 (MQ=255)
      ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:2703192/70‑1 (MQ=255)
      ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:422764/70‑1 (MQ=255)
      ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:1545122/70‑1 (MQ=255)
       tctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:140882/1‑71 (MQ=255)
           ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:2235050/1‑70 (MQ=255)
           ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:2848944/1‑70 (MQ=255)
           ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:15778/1‑70 (MQ=255)
             aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:2862998/1‑65 (MQ=255)
             aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:4615940/1‑65 (MQ=255)
              atcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                                           <  1:2602926/68‑1 (MQ=255)
               tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa                                                         >  1:5350984/2‑69 (MQ=255)
                 ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                       <  1:4987377/71‑1 (MQ=255)
                 ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                       <  1:1087263/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:3514677/1‑70 (MQ=255)
                      aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:3072811/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:5350784/50‑2 (MQ=255)
                                   gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:3023090/50‑2 (MQ=255)
                                   gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:3992153/50‑2 (MQ=255)
                                    aaTATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 <  1:2917339/42‑1 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:1891867/1‑71 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:1906093/1‑71 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   <  1:2240814/71‑1 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:665995/1‑71 (MQ=255)
                                     atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:5192344/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:391025/1‑71 (MQ=255)
                                                         cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:3098297/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:4718488/71‑1 (MQ=255)
                                                             agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:3949063/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)
                                                                      cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg  >  1:4786541/1‑71 (MQ=255)
                                                                      cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg  >  1:1007875/1‑71 (MQ=255)
                                                                      |                                                                      
CATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAG  >  minE/2677345‑2677485

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