Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,554,371 G→A 53.8% L89L (CTC→CTT 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,3710GA53.8% 8.3 / 78.1 91L89L (CTC→CTTenvCprotease with a role in cell division
Reads supporting (aligned to +/- strand):  ref base G (26/16);  new base A (18/31);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCCAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGCTTCTTC  >  minE/2554315‑2554424
                                                        |                                                     
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2274620/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2308866/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1744151/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:234418/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:3147799/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:344898/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:416736/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:44421/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:1202171/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1186044/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2670687/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2220498/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:1167499/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:821932/1‑69 (MQ=255)
tgctCTAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGTTTTGc                                           >  1:2216472/1‑69 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2961031/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1743104/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1754142/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2282230/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1773414/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2565059/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1873859/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1888388/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2136500/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:274140/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:1304659/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:246307/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2468640/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2408550/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:609840/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:58978/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:575225/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:2804135/57‑1 (MQ=255)
         ttGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTCAGTTGATTAAGCGTGttt                                              <  1:2797881/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:2505200/69‑1 (MQ=255)
                     gcgTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTg                      <  1:685408/69‑1 (MQ=255)
                     gcgTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTg                      <  1:882413/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:1614813/54‑1 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTa                            >  1:1919689/1‑54 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:1498901/1‑64 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:2181536/1‑64 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:458189/1‑64 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:1384327/1‑64 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCtt                  >  1:872662/1‑64 (MQ=255)
                              tcatcaATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCt                   >  1:1019603/1‑63 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1478414/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:536048/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:3107578/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:632526/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1080727/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1068434/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1658310/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:824117/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2084783/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2671647/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:1037842/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2609846/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:898896/1‑48 (MQ=255)
                                  caATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCt                              >  1:2245262/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:823068/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:1610464/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:3414835/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:3371581/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2597466/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2940920/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:1638298/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2652252/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGcttct    <  1:2733377/65‑1 (MQ=255)
                                           tATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGCttcttc  <  1:2846301/67‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCTCAGCTTACg                          <  1:2830677/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:2494857/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:2386249/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:933879/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:1760053/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:1624157/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:3099257/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:1392415/35‑1 (MQ=255)
                                                   tgattgaGCGTGTTTTGCGTTTCGCGCAGCTTACg                          <  1:1561820/35‑1 (MQ=255)
                                                        |                                                     
TGCTCCAGTTTGGCAATCGACGCGTTCATCTCATCAATCTGTTTATTGAGTTGATTGAGCGTGTTTTGCGTTTCGCGCAGCTTACGGGTGGCTTCAGAGATCGCTTCTTC  >  minE/2554315‑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: