Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,676 T→C 24.6% intergenic (+58/‑317) 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,6760TC24.6% 117.6 / 32.5 57intergenic (+58/‑317)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (16/27);  new base C (6/8);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

TATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTA  >  minE/2493606‑2493746
                                                                      |                                                                      
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)
        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:2386982/1‑52 (MQ=255)
                                 tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                          >  1:6608/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:1751399/1‑70 (MQ=255)
                                       aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                  <  1:1127597/70‑1 (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:318566/46‑1 (MQ=255)
                                          cgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAttt                                                       <  1:3251488/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:2565870/51‑1 (MQ=255)
                                            cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                                <  1:1709508/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)
                                                                      ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCta  >  1:1279665/1‑71 (MQ=255)
                                                                      ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCt   >  1:1898855/1‑70 (MQ=255)
                                                                      |                                                                      
TATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTA  >  minE/2493606‑2493746

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: