Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 711,127 A→G 29.0% D96D (GAT→GAC solA ← N‑methyltryptophan oxidase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE711,1270AG29.0% 58.0 / 26.0 62D96D (GAT→GACsolAN‑methyltryptophan oxidase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (27/17);  new base G (4/14);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

TCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGG  >  minE/711062‑711190
                                                                 |                                                               
tCGAGCTTTTTAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCa                                                              >  1:1222167/1‑69 (MQ=255)
tCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCa                                                              >  1:3361276/1‑69 (MQ=255)
tCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCa                                                              >  1:3114483/1‑69 (MQ=255)
tCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCa                                                              >  1:1128223/1‑69 (MQ=255)
tCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCa                                                              >  1:1145706/1‑69 (MQ=255)
tCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCa                                                              >  1:1828784/1‑69 (MQ=255)
                 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:2841336/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:517966/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1885946/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:185930/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1235327/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:249056/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:2565831/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:2706738/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1721927/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:736621/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:3045072/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:304587/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1434624/69‑1 (MQ=255)
                 gAGTTGCCACTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTAAt                                             <  1:1356903/69‑1 (MQ=255)
                    ttGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGAc                                          >  1:2910220/1‑69 (MQ=255)
                    ttGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGAc                                          >  1:1446687/1‑69 (MQ=255)
                            gTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCagag                                   >  1:2941748/1‑68 (MQ=255)
                            gTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCagag                                   >  1:1462301/1‑68 (MQ=255)
                            gTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCagag                                   >  1:2162732/1‑68 (MQ=255)
                                         ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:1059978/69‑1 (MQ=255)
                                         ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:487201/69‑1 (MQ=255)
                                         ggCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATggg                     <  1:2423460/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:651304/1‑43 (MQ=25)
                                         ggCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTa                                               >  1:1109369/1‑43 (MQ=25)
                                         ggCGACGTTGGCGAGAAATGCGGAGTCAGCCGGACCGAGGTTa                                               >  1:2748530/1‑43 (MQ=25)
                                            gACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGAtc                  >  1:2528474/1‑69 (MQ=255)
                                            gACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGAtc                  >  1:1695774/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:2255308/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:1280387/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:839331/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:1726768/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:1244584/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:2612548/1‑69 (MQ=255)
                                                  ggCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCttcttc            >  1:608916/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:1775995/69‑1 (MQ=255)
                                                            gtgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:1747448/69‑1 (MQ=255)
                                                             tgGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  <  1:1954179/68‑1 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:312132/1‑67 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:2823281/1‑67 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:2465538/1‑67 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:1476283/1‑67 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCgg  >  1:1279087/1‑67 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCg   >  1:3249827/1‑66 (MQ=255)
                                                              ggAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTACGTTGTGGCgg  >  1:166955/1‑67 (MQ=255)
                                                                 |                                                               
TCGAGCTTTTCAACGTTGAGTTGCCATTGTTCGGCGCTGTGGGCGACGTTGGCGAGAAATGTGGAATCAGCCGGGCCAAGGTTAATGACGCCAGAGCGTACAAAAATGGGATCTTCTTCGTTGTGGCGG  >  minE/711062‑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: