Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,675 C→T 22.9% intergenic (+57/‑318) gyrB → / → ibpA DNA gyrase, subunit B/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,493,6750CT22.9% 136.8 / 30.0 61intergenic (+57/‑318)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (18/29);  new base T (6/8);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGC  >  minE/2493606‑2493744
                                                                     |                                                                     
tatCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCt                                                                      >  1:2587836/1‑71 (MQ=255)
tatCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCt                                                                      <  1:3694531/71‑1 (MQ=255)
  tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                     <  1:3356194/70‑1 (MQ=255)
    gATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCa                                                                  <  1:3012516/71‑1 (MQ=255)
        tttAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATc                                                                       <  1:1464544/62‑1 (MQ=255)
        tttAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAAt                                                                <  1:2921712/69‑1 (MQ=255)
          tAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTa                                                             <  1:359626/70‑1 (MQ=255)
                 gcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                   <  1:3775165/57‑1 (MQ=255)
                 gcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAttt                                                     <  1:2201371/71‑1 (MQ=255)
                              gCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          >  1:177543/1‑69 (MQ=255)
                              gCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTa                                        <  1:928018/71‑1 (MQ=255)
                               aatGCCTGATGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAg                                       >  1:2701367/3‑71 (MQ=14)
                                atGCCTGATGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAgg                                      <  1:1678206/70‑1 (MQ=14)
                                 tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                        >  1:6608/1‑52 (MQ=255)
                                 tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                        >  1:2386982/1‑52 (MQ=255)
                                  gCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                     >  1:2278556/1‑38 (MQ=255)
                                   ccTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   >  1:741336/1‑71 (MQ=255)
                                    cTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAtttt                                                    <  1:1089963/53‑1 (MQ=255)
                                    cTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   >  1:229203/1‑70 (MQ=255)
                                      gATGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGa                                 >  1:277434/1‑70 (MQ=14)
                                       aTGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGAt                                <  1:1366568/70‑1 (MQ=14)
                                       aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                <  1:1127597/70‑1 (MQ=255)
                                       aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                >  1:1751399/1‑70 (MQ=255)
                                        tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAGTGTATTGAATTTGCACGATTTTGTAGGCCGGAt                                >  1:3531496/1‑69 (MQ=14)
                                        tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                >  1:3830590/1‑69 (MQ=255)
                                          cgcTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAg                             <  1:1473834/70‑1 (MQ=14)
                                          cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAttt                                                     <  1:3251488/46‑1 (MQ=255)
                                          cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAttt                                                     <  1:318566/46‑1 (MQ=255)
                                          cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAgg                            <  1:1981069/71‑1 (MQ=255)
                                           gcTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAAtttg                                                    <  1:2458354/46‑2 (MQ=14)
                                           gcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                           <  1:2128235/55‑1 (MQ=255)
                                           gcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          >  1:2761066/1‑56 (MQ=255)
                                            cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                              <  1:1559875/51‑1 (MQ=255)
                                            cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                              <  1:1709508/51‑1 (MQ=255)
                                            cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                              <  1:2565870/51‑1 (MQ=255)
                                            cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCg                                   <  1:3240778/62‑1 (MQ=255)
                                             tACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGt                         >  1:1674431/1‑71 (MQ=255)
                                               cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGa                                                         <  1:632998/37‑1 (MQ=17)
                                               cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGtt                        >  1:981727/1‑70 (MQ=14)
                                                 cTTATCAGGCCTACGAAAATTCTGCAATGTATTGa                                                         >  1:3149691/1‑35 (MQ=14)
                                                 cTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAc                     <  1:3193602/71‑1 (MQ=255)
                                                   tATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:354434/48‑1 (MQ=255)
                                                   tATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                   >  1:3533810/1‑71 (MQ=255)
                                                     tCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTa                                        >  1:1000507/1‑48 (MQ=255)
                                                        ggCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCAt              <  1:1782256/71‑1 (MQ=14)
                                                        ggCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt              <  1:3135301/71‑1 (MQ=255)
                                                        ggCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt              >  1:3576027/1‑71 (MQ=255)
                                                         gCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATc             <  1:3501910/71‑1 (MQ=11)
                                                         gCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc             <  1:1685497/71‑1 (MQ=255)
                                                         gCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATc             <  1:893433/71‑1 (MQ=255)
                                                          ccTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATcc            <  1:1653805/71‑1 (MQ=11)
                                                           cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTa                                        >  1:3602198/1‑42 (MQ=14)
                                                            tACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                           <  1:258896/38‑1 (MQ=255)
                                                             aCGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg          >  1:3912605/1‑70 (MQ=255)
                                                               gAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                          <  1:989866/36‑1 (MQ=255)
                                                               gAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATaa                              <  1:3297536/48‑1 (MQ=255)
                                                               gAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa       <  1:1248346/71‑1 (MQ=255)
                                                                 aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                        >  1:1134761/1‑36 (MQ=255)
                                                                 aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaag      <  1:1423851/70‑1 (MQ=255)
                                                                  aaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATa                               <  1:3883058/44‑1 (MQ=255)
                                                                   aTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaag   >  1:1047723/1‑71 (MQ=255)
                                                                    tCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGc  >  1:1307273/1‑71 (MQ=255)
                                                                     |                                                                     
TATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGC  >  minE/2493606‑2493744

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: