Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,494,512 A→G 7.4% intergenic (+106/‑6) ibpA → / → ibpB heat shock chaperone/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,494,5120AG 7.4% 226.3 / 7.5 81intergenic (+106/‑6)ibpA/ibpBheat shock chaperone/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (42/33);  new base G (3/3);  total (45/36)
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.52e-01

TGCTCGCCGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCAACGC  >  minE/2494443‑2494582
                                                                     |                                                                      
tGCTCGCCGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAAt                                                                       >  1:2119663/1‑71 (MQ=255)
     gCCGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACt                                                                   <  1:1707597/70‑1 (MQ=255)
     gCCGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTa                                                                  >  1:2063649/1‑71 (MQ=255)
     gCCGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTa                                                                  <  1:3193395/71‑1 (MQ=255)
      ccGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAGTGACTAt                                                                 >  1:1305788/1‑71 (MQ=255)
         tCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCTCTTCTTAGAAGGAGAAATGACTATGCg                                                              >  1:1325762/1‑71 (MQ=255)
         tCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTATAAGGAGAAATGACTATGCg                                                              >  1:3478958/1‑71 (MQ=255)
         tCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCg                                                              >  1:2362154/1‑71 (MQ=255)
         tCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCg                                                              >  1:2584348/1‑71 (MQ=255)
         tCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCg                                                              >  1:760328/1‑71 (MQ=255)
           aGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGa                                                                     >  1:3666472/1‑62 (MQ=255)
           aGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGt                                                             <  1:3917379/70‑1 (MQ=255)
           aGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGt                                                             <  1:2820033/70‑1 (MQ=255)
           aGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGt                                                             <  1:1832134/70‑1 (MQ=255)
             ggAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAAc                                                          <  1:1173884/71‑1 (MQ=255)
              gAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAAt                                                                       <  1:12845/57‑1 (MQ=255)
                gCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTc                                                       <  1:1512655/71‑1 (MQ=255)
                 cATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCg                                                      <  1:976580/71‑1 (MQ=255)
                    atGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGa                                                     >  1:990624/1‑69 (MQ=255)
                    atGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGAtt                                                   >  1:891573/1‑71 (MQ=255)
                     tGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGAtt                                                   >  1:2881303/1‑70 (MQ=255)
                      gCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTa                                                 >  1:3238165/1‑71 (MQ=255)
                       cGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTa                                                 <  1:358241/70‑1 (MQ=255)
                       cGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTAt                                                >  1:1652131/1‑71 (MQ=255)
                        gAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTa                                                                  >  1:1888849/1‑52 (MQ=255)
                         aaTCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGa                                                     <  1:2040807/64‑1 (MQ=255)
                         aaTCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATcc                                              >  1:2125858/1‑71 (MQ=255)
                          aTCTTCGGATTTGCAGGTGCTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTAt                                                >  1:3420534/1‑68 (MQ=255)
                          aTCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTAt                                                >  1:3712549/1‑68 (MQ=255)
                          aTCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATccc                                             <  1:1372530/71‑1 (MQ=255)
                           tCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTa                                                                  >  1:1279856/1‑49 (MQ=255)
                             ttCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGa                                                                     >  1:3327167/1‑44 (MQ=255)
                               cGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGaaa                                                                        >  1:2678744/1‑39 (MQ=255)
                               cGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGaaa                                                                        >  1:2736987/1‑39 (MQ=255)
                               cGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATcccc                                            <  1:2811281/67‑1 (MQ=255)
                               cGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACt                                         >  1:3261431/1‑70 (MQ=255)
                               cGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACt                                         <  1:2828619/70‑1 (MQ=255)
                                  aTTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACt                                                         >  1:627460/1‑51 (MQ=255)
                                    ttGCAGGTACTTACTCGCTTCTTAGAAGGAGAAGTGACTATGCGTAACTTCGATTTATCCCCACTGATGc                                    >  1:16717/1‑70 (MQ=255)
                                        aGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCCCTGATGCGTCaa                               <  1:72917/71‑1 (MQ=255)
                                          gTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATg                             >  1:223037/1‑71 (MQ=255)
                                           tacttactCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATgg                            <  1:3045851/71‑1 (MQ=255)
                                           tacttactCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATgg                            <  1:2052169/71‑1 (MQ=255)
                                            acttactCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGa                           >  1:3537271/1‑71 (MQ=255)
                                             cttactCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATg                             >  1:2940976/1‑68 (MQ=255)
                                              ttactCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGAttt                                                  <  1:2278467/46‑1 (MQ=255)
                                              ttactCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGAt                                      <  1:2751381/58‑1 (MQ=255)
                                                actCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGa                                                     >  1:3568137/1‑41 (MQ=255)
                                                actCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGa                                                     >  1:368564/1‑41 (MQ=255)
                                                actCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTAt                                                >  1:1583712/1‑46 (MQ=255)
                                                actCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCgg                       >  1:2810794/1‑71 (MQ=255)
                                                  tCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGAttt                                                  <  1:3109130/42‑1 (MQ=255)
                                                  tCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACt                                         <  1:2836264/51‑1 (MQ=255)
                                                    gCTTCTTAGAAGGAGAAGTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGtt                     <  1:1540883/69‑1 (MQ=255)
                                                    gCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGtttt                   >  1:2696143/1‑71 (MQ=255)
                                                    gCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGtttt                   <  1:1161369/71‑1 (MQ=255)
                                                    gCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGtttt                   <  1:1130899/71‑1 (MQ=255)
                                                     cttcttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCaa                               >  1:1226756/1‑58 (MQ=255)
                                                      ttcttAGAAGGAGAAGTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATgg                            <  1:1761256/60‑1 (MQ=255)
                                                      ttcttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCAc                                          >  1:1793955/1‑46 (MQ=255)
                                                      ttcttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGa                 >  1:3387132/1‑71 (MQ=255)
                                                      ttcttAGAAGGAGAAATGACTATGAGTAACTTCGATTTATCCCCa                                           >  1:167155/1‑45 (MQ=255)
                                                       tcttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATCCGTCAATGGATCGGTTTTGa                 >  1:2312946/1‑70 (MQ=255)
                                                        cttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCa                                >  1:806855/1‑54 (MQ=255)
                                                         ttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGt                      <  1:234656/63‑1 (MQ=255)
                                                         ttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGt                      <  1:2976043/63‑1 (MQ=255)
                                                         ttAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGt                      <  1:3693833/63‑1 (MQ=255)
                                                           aGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGa                           >  1:1607513/1‑56 (MQ=255)
                                                            gAAGGAGAAGTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACa               <  1:2493123/67‑1 (MQ=255)
                                                            gAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACaa              <  1:3146274/68‑1 (MQ=255)
                                                            gAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAAc            <  1:3061622/70‑1 (MQ=255)
                                                             aaGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTg          <  1:3911412/71‑1 (MQ=255)
                                                              aGGAGAAGTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTgg         >  1:672578/1‑71 (MQ=255)
                                                              aGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTgg         <  1:2506644/71‑1 (MQ=255)
                                                               ggAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACa               <  1:1512739/64‑1 (MQ=255)
                                                               ggAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTgg         >  1:299768/1‑70 (MQ=255)
                                                               ggAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGc        <  1:669752/71‑1 (MQ=255)
                                                                 agaAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGcc       >  1:3793462/1‑70 (MQ=255)
                                                                  gaAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGtttt                   >  1:265440/1‑57 (MQ=255)
                                                                  gaAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACa               <  1:1022949/61‑1 (MQ=255)
                                                                  gaAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGcc       >  1:1724057/1‑69 (MQ=255)
                                                                  gaAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCaa     >  1:1940839/1‑71 (MQ=255)
                                                                  gaAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCaa     >  1:2231287/1‑71 (MQ=255)
                                                                     aTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACaaa             >  1:14686/1‑60 (MQ=255)
                                                                     aTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTgg         >  1:773289/1‑64 (MQ=255)
                                                                     aTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCAACg   >  1:578867/1‑70 (MQ=255)
                                                                     aTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCAACGc  <  1:1694447/71‑1 (MQ=255)
                                                                     aTGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCAACGc  <  1:3125373/71‑1 (MQ=255)
                                                                     |                                                                      
TGCTCGCCGTCAGGGAGCATATGCGAATCTTCGGATTTGCAGGTACTTACTCGCTTCTTAGAAGGAGAAATGACTATGCGTAACTTCGATTTATCCCCACTGATGCGTCAATGGATCGGTTTTGACAAACTGGCCAACGC  >  minE/2494443‑2494582

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: