Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,316,188 A→G 37.1% E502E (GAA→GAG glpK → glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,316,1880AG37.1% 69.9 / 78.4 121E502E (GAA→GAGglpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base G (35/10);  total (70/51)
Fisher's exact test for biased strand distribution p-value = 6.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 122 total aligned reads displayed.

ACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGCAACCCCG  >  minE/2316127‑2316256
                                                             |                                                                    
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1851393/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               >  1:400322/1‑69 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1593552/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1436377/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:3320439/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:3192551/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:3069432/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:565520/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1829109/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1158985/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:852884/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               >  1:1148805/1‑69 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1855286/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               >  1:2058195/1‑69 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:2666088/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:2119803/69‑1 (MQ=255)
aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:231877/69‑1 (MQ=255)
      aCGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGc                                                          >  1:101841/1‑68 (MQ=255)
      aCGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGc                                                          >  1:14539/1‑68 (MQ=255)
        gCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCg                                                        <  1:511232/68‑1 (MQ=255)
        gCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCg                                                        <  1:2248352/68‑1 (MQ=255)
        gCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCg                                                        <  1:913707/68‑1 (MQ=255)
             cTGGAAAAAAGGGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:428772/56‑1 (MQ=255)
             cTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATgca                                                  >  1:2520634/1‑67 (MQ=255)
             cTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGaa                                                  >  1:1527475/1‑69 (MQ=255)
             cTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGaa                                                  >  1:521671/1‑69 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:2509065‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:2661354‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1039767‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:2754557‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1998055‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1111567‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:686114‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1733121‑M2/69‑20 (MQ=255)
                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1311068‑M2/69‑20 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:486029/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:3057286/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:342234/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:2308456/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2524643/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:3214827/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2530378/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2608624/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2645272/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:3208947/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2655493/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:3120281/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:853991/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2807395/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:3078246/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2934068/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2988868/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:551255/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:781933/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:2365848/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2336938/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1616157/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2233562/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:1509542/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:2196335/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1645191/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1506235/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:1472134/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:1406415/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2015855/1‑69 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1744783/69‑1 (MQ=255)
                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1954386/69‑1 (MQ=255)
                            tAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:73130/41‑1 (MQ=255)
                              aaCGCGCGATTGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:152771/39‑1 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGATGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:3431201‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcttcttat                                >  1:3744‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1197676‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1258210‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1789789‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1061513‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:172837‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:2034669‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:2087966‑M2/1‑35 (MQ=255)
                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:2114535‑M2/1‑35 (MQ=255)
                                cgcgcgATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:2099689/37‑1 (MQ=255)
                                cgcgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1015131‑M2/1‑34 (MQ=255)
                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:919467‑M2/1‑32 (MQ=255)
                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:130344‑M2/1‑32 (MQ=255)
                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1025624‑M2/1‑32 (MQ=255)
                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1696489‑M2/1‑32 (MQ=255)
                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3156531‑M2/1‑32 (MQ=255)
                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1737154‑M2/1‑32 (MQ=255)
                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3197549‑M2/1‑31 (MQ=255)
                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1476526‑M2/1‑31 (MQ=255)
                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:126346‑M2/1‑31 (MQ=255)
                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3036409‑M2/1‑31 (MQ=255)
                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:238090‑M2/1‑31 (MQ=255)
                                       tGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1152198‑M2/1‑27 (MQ=255)
                                        ggCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTccc                        >  1:2223753/1‑68 (MQ=255)
                                              ggAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3043948‑M2/1‑20 (MQ=255)
                                               gaagaaCACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1631303‑M2/1‑19 (MQ=255)
                                                   aaCACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:3311120‑M2/1‑15 (MQ=255)
                                                      acgacgAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1946650‑M2/1‑12 (MQ=255)
                                                      acgacgAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGc         <  1:2891201/69‑1 (MQ=255)
                                                             ataatGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGAAACCCCg  <  1:384538/69‑1 (MQ=255)
                                                             |                                                                    
ACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGCAACCCCG  >  minE/2316127‑2316256

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: