Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,331,307 A→G 68.2% I93I (ATT→ATC cpsG ← phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,331,3070AG68.2% 40.2 / 20.3 47I93I (ATT→ATCcpsGphosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (5/10);  new base G (14/18);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

GAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATA  >  minE/1331243‑1331373
                                                                |                                                                  
gAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgc                                                                <  1:1371157/69‑1 (MQ=255)
gAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgc                                                                <  1:1396818/69‑1 (MQ=255)
gAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgc                                                                <  1:499248/69‑1 (MQ=255)
gAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgc                                                                <  1:3334682/69‑1 (MQ=255)
gAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgc                                                                <  1:3005658/69‑1 (MQ=255)
gAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgc                                                                <  1:2549971/69‑1 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:982628/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:1222465/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:897224/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:805447/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:571594/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:2982465/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:2854249/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:2683491/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:2676858/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:2355276/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:1422855/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:2162857/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:1599911/1‑69 (MQ=255)
          gcgcACCAGTTTCATGCCGTTGTAATCCATCGGGTTATGGCTGGCGGTAACTTCGATGCCGCCATCCAc                                                      >  1:1997813/1‑69 (MQ=255)
               ccAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCga                                                 <  1:3267885/69‑1 (MQ=255)
               ccAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCga                                                 <  1:3208687/69‑1 (MQ=255)
               ccAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCga                                                 <  1:2889421/69‑1 (MQ=255)
               ccAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCga                                                 <  1:2085940/69‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGTAACGTGGCGAAATAg                             <  1:3327543/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:495885/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:2466786/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:910727/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1532704/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:874026/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:853066/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1572441/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:733661/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:2675095/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1613153/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:2834256/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:3420888/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1631283/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1685960/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1692993/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:2700651/59‑1 (MQ=255)
                                             tATGGCTGGCGGTAACTTCGATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                             <  1:1884593/59‑1 (MQ=255)
                                                          aaCTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATAcc      >  1:3056450/1‑69 (MQ=255)
                                                              tCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAata  >  1:2164512/1‑69 (MQ=255)
                                                              tCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAata  >  1:825232/1‑69 (MQ=255)
                                                              tCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAata  >  1:3036113/1‑69 (MQ=255)
                                                              tCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAata  >  1:975302/1‑69 (MQ=255)
                                                                |                                                                  
GAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATA  >  minE/1331243‑1331373

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