Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,836,489 G→C 67.5% A277A (GCC→GCG yjeK ← predicted lysine aminomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,836,4890GC67.5% 139.2 / 97.4 129A277A (GCC→GCGyjeKpredicted lysine aminomutase
Reads supporting (aligned to +/- strand):  ref base G (37/5);  new base C (33/54);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 5.11e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 137 total aligned reads displayed.

CACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGT  >  minE/2836422‑2836549
                                                                   |                                                            
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:31688/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:2958768/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:287414/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:283496/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:1226007/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:2430639/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:1671870/69‑1 (MQ=255)
cactcacCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCgg                                                             <  1:1570523/69‑1 (MQ=255)
        aTAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGaac                                                     <  1:3265966/69‑1 (MQ=255)
         tAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGaac                                                     <  1:3002456/68‑1 (MQ=255)
           aaaTGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGaacaa                                                   >  1:599726/1‑68 (MQ=255)
           aaaTGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGaacaa                                                   >  1:2038114/1‑68 (MQ=255)
           aaaTGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGaacaa                                                   >  1:749700/1‑68 (MQ=255)
           aaaTGCGCCGCGCCCTGTACTTTATCGAGAACATGCAGGTAATACGGCATTACGCCGGCATCGaac                                                     >  1:1019044/1‑66 (MQ=255)
                            tACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAgg                                      >  1:620/1‑64 (MQ=255)
                            tACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGttt                                   >  1:154598/1‑67 (MQ=255)
                            tACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGttt                                   >  1:48062/1‑67 (MQ=255)
                            tACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGttt                                   >  1:1597181/1‑67 (MQ=255)
                              cTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCa                               >  1:602290/1‑69 (MQ=255)
                              cTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCa                               >  1:798947/1‑69 (MQ=255)
                              cTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCa                               >  1:1819382/1‑69 (MQ=255)
                              cTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCa                               >  1:1706561/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:2264688/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:2235359/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:3365470/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:3087186/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:205650/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:1835633/1‑69 (MQ=255)
                                    cGAGCACATGCAGGTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGttt                         >  1:3330262/1‑69 (MQ=255)
                                           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGtt                  >  1:2439921/1‑69 (MQ=255)
                                           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGtt                  >  1:8713/1‑69 (MQ=255)
                                           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGtt                  >  1:1311767/1‑69 (MQ=255)
                                           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGtt                  >  1:2969048/1‑69 (MQ=255)
                                           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGtt                  >  1:1690016/1‑69 (MQ=255)
                                           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAAAGCATTACTCAGGTTTGCCAGCGTTTGTGCGtt                  >  1:1223375/1‑69 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGTTTTGCCAGCGTTTGCGCGTTATCg              <  1:1664241/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGTCAGCGTTTGCGCGTTATCg              <  1:1225780/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:1181463/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:1054738/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:273908/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:3306582/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:496475/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:2665164/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:2525703/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:1669179/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:1408956/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:1610757/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:2083445/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:2077750/68‑1 (MQ=255)
                                                ggTAATACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:1405045/68‑1 (MQ=255)
                                                 gTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGtt            <  1:1620800/69‑1 (MQ=255)
                                                 gTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGtt            <  1:2112080/69‑1 (MQ=255)
                                                 gTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGtt            <  1:3363283/69‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:1538619/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:3252704/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:3149861/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:3310230/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:1576160/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:922207/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:2861450/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:48579/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:2282212/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:2667918/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:885463/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:520911/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:2357477/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACt                                          >  1:899429/1‑37 (MQ=38)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:630055/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:937490/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:862962/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:3285503/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:968321/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:2410142/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:1035051/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:1135026/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:1286156/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:1439721/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:1716713/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:2169912/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:253241/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:2569465/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:3192648/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:2968190/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:2886750/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:2774197/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:281635/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTAGCCAGCGt                           <  1:1552424/52‑1 (MQ=255)
                                                   aaTACGGCATCACGCCCGCATCGAACAACGAATTACt                                          >  1:1939076/1‑37 (MQ=25)
                                                    aTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:230221/51‑1 (MQ=255)
                                                    aTACGGCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:579266/51‑1 (MQ=255)
                                                      aCGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACAt       >  1:3026657/1‑69 (MQ=255)
                                                      aCGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACAt       >  1:963285/1‑69 (MQ=255)
                                                       cGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATc      >  1:2011956/1‑69 (MQ=255)
                                                       cGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATc      >  1:2729465/1‑69 (MQ=255)
                                                       cGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATc      >  1:3220456/1‑69 (MQ=255)
                                                        ggCATCACGCCCGCATCGAACAACGCATTACTCAGGTTTGCCAGCGt                           <  1:918620/47‑1 (MQ=39)
                                                             caCGCCCTCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGCGCGTTATCg              <  1:484139/54‑1 (MQ=255)
                                                              aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTca          >  1:3399781/1‑58 (MQ=255)
                                                              aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGt  >  1:1170055/1‑66 (MQ=255)
                                                              aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGt  >  1:2706764/1‑66 (MQ=255)
                                                                   |                                                            
CACTCACCATAAAATGCGCCGCGCCCTGTACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGT  >  minE/2836422‑2836549

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: