Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 711,136 G→A 32.0% G93G (GGC→GGT solA ← N‑methyltryptophan oxidase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE711,1360GA32.0% 46.3 / 26.1 56G93G (GGC→GGTsolAN‑methyltryptophan oxidase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base A (4/14);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGG  >  minE/711079‑711190
                                                         |                                                      
gAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCTGGCCAAGGTTAAt                                             <  1:1713897/69‑1 (MQ=255)
gAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAAt                                             <  1:1804041/69‑1 (MQ=255)
gAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAAt                                             <  1:1728315/69‑1 (MQ=255)
gAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAAt                                             <  1:1470783/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:736621/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:185930/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1885946/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:249056/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:2565831/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1721927/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:2706738/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:2841336/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1434624/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1356903/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:3045072/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:304587/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1235327/69‑1 (MQ=255)
gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:517966/69‑1 (MQ=255)
   ttGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGAc                                          >  1:2910220/1‑69 (MQ=255)
   ttGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGAc                                          >  1:1446687/1‑69 (MQ=255)
           gTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCagag                                   >  1:1462301/1‑68 (MQ=255)
           gTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCagag                                   >  1:2941748/1‑68 (MQ=255)
           gTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCagag                                   >  1:2162732/1‑68 (MQ=255)
                        ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:487201/69‑1 (MQ=255)
                        ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:2423460/69‑1 (MQ=255)
                        ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:1059978/69‑1 (MQ=255)
                        ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:1846877/69‑1 (MQ=255)
                        ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:1392887/69‑1 (MQ=255)
                        ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:1750019/69‑1 (MQ=255)
                        ggCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTa                                               >  1:3366557/1‑43 (MQ=25)
                        ggCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTa                                               >  1:1109369/1‑43 (MQ=25)
                        ggCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTa                                               >  1:2748530/1‑43 (MQ=25)
                        ggCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTa                                               >  1:651304/1‑43 (MQ=25)
                           gACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGAtc                  >  1:2528474/1‑69 (MQ=255)
                           gACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGAtc                  >  1:1695774/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:1726768/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:2612548/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:608916/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:2255308/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:1244584/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:839331/1‑69 (MQ=255)
                                 ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:1280387/1‑69 (MQ=255)
                                           gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:102184/69‑1 (MQ=255)
                                           gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:432860/69‑1 (MQ=255)
                                           gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:351698/69‑1 (MQ=255)
                                           gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:843732/69‑1 (MQ=255)
                                           gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:1747448/69‑1 (MQ=255)
                                           gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:1775995/69‑1 (MQ=255)
                                            tgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:1954179/68‑1 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:312132/1‑67 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:1279087/1‑67 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:1476283/1‑67 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:2823281/1‑67 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:2465538/1‑67 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCg   >  1:3249827/1‑66 (MQ=255)
                                             ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTACGTTGTGGCgg  >  1:166955/1‑67 (MQ=255)
                                                         |                                                      
GAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGG  >  minE/711079‑711190

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: