Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,554,362 G→C 52.7% L92L (CTC→CTG envC ← protease with a role in cell division

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,554,3620GC52.7% 7.8 / 111.7 93L92L (CTC→CTGenvCprotease with a role in cell division
Reads supporting (aligned to +/- strand):  ref base G (27/17);  new base C (18/31);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CCTGGGCGGCTTTTTGCTGCTCCAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGCTTCTTC  >  minE/2554298‑2554424
                                                                |                                                              
ccTGGGCGGCTTTTTGCTGCTCCAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTGAGtt                                                            >  1:2719460/1‑69 (MQ=255)
   gggCGGCTTTTTGCTGCTCCAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTGAGttg                                                           <  1:2279538/67‑1 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:821932/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1744151/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:19461/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1983639/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2706489/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2670687/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:3147799/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2216472/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2220498/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:344898/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:416736/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2274620/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1202171/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1186044/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1167499/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:44421/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2308866/1‑69 (MQ=255)
                 tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:234418/1‑69 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2468640/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2565059/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2282230/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2136500/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2408550/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1888388/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:246307/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1873859/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1754142/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:609840/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:58978/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:575225/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1169026/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1212643/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1300762/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1304659/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1450366/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1456388/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:3413859/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2961031/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2804135/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2797881/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:274140/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1743104/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1773414/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTCttt                                              <  1:203316/57‑1 (MQ=255)
                          ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATCAAGCGTGttt                                              <  1:1681677/57‑1 (MQ=255)
                              cAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1689299/53‑1 (MQ=255)
                                   gACGCGTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACgg                         >  1:972738/1‑69 (MQ=255)
                                   gACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:3120414/48‑1 (MQ=255)
                                    aCGCTTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2743124/47‑1 (MQ=255)
                                      gcgTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTg                      <  1:882413/69‑1 (MQ=255)
                                      gcgTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTg                      <  1:2505200/69‑1 (MQ=255)
                                      gcgTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTg                      <  1:685408/69‑1 (MQ=255)
                                            atctcatcAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2553423/39‑1 (MQ=38)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTa                            <  1:1031950/54‑1 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTa                            >  1:1415883/1‑54 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTa                            >  1:1919689/1‑54 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTa                            <  1:1614813/54‑1 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:2181536/1‑64 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:872662/1‑64 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:458189/1‑64 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:1384327/1‑64 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:1498901/1‑64 (MQ=255)
                                               tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCt                   >  1:1019603/1‑63 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:632526/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1037842/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1068434/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1080727/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:536048/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2671647/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2084783/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1478414/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2609846/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:824117/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1658310/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:3107578/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2245262/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:898896/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:950805/1‑48 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACg                          >  1:2110981/1‑52 (MQ=255)
                                                   caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCg         >  1:112969/1‑69 (MQ=255)
                                                      tCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGct       <  1:2450087/68‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:3414835/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2652252/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:3371581/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:823068/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2733377/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2597466/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:1610464/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2940920/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:1638298/65‑1 (MQ=255)
                                                            tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGCttcttc  <  1:2846301/67‑1 (MQ=255)
                                                                |                                                              
CCTGGGCGGCTTTTTGCTGCTCCAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGCTTCTTC  >  minE/2554298‑2554424

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: