Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,652 A→C 13.5% intergenic (+34/‑341) 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,6520AC13.5% 116.9 / 15.3 52intergenic (+34/‑341)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (19/26);  new base C (3/4);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

AGAGAACGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGT  >  minE/2493582‑2493721
                                                                      |                                                                     
agagAACGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTa                                                                       >  1:960644/1‑71 (MQ=255)
     aCGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCtt                                                                  <  1:3768355/71‑1 (MQ=255)
     aCGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCtt                                                                  <  1:3451787/71‑1 (MQ=255)
     aCGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCtt                                                                  <  1:1938221/71‑1 (MQ=255)
      cGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTa                                                                 <  1:3753495/71‑1 (MQ=255)
       gCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTAt                                                                <  1:1334290/71‑1 (MQ=255)
        cccTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCtt                                                                  >  1:188820/1‑68 (MQ=255)
                gcggcgAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGcc                                                          >  1:1506991/1‑68 (MQ=255)
                  ggcgAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTa                                                                 <  1:3703127/59‑1 (MQ=255)
                    taaaTATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATc                                                               >  1:3011991/3‑59 (MQ=255)
                    taaaTATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCAGATGCGCTACGCTTATc                                                               >  1:2572793/3‑59 (MQ=255)
                      aaTATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGaaaa                                                  >  1:372067/1‑70 (MQ=255)
                       atatCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCa                                                              >  1:1409141/1‑57 (MQ=255)
                       atatCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAg                                                             <  1:150729/58‑1 (MQ=255)
                        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)
                                  tAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTAt                                                                >  1:416423/1‑44 (MQ=255)
                                  tAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTAt                                                                >  1:2190821/1‑44 (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: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:2565870/51‑1 (MQ=255)
                                                                    cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                       <  1:1559875/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)
                                                                      |                                                                     
AGAGAACGCCCTGAAAGCGGCGAATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGT  >  minE/2493582‑2493721

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