Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,980,802 C→A 50.6% A103A (GCC→GCA yqgF → predicted Holliday junction resolvase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,980,8020CA50.6% 5.5 / 90.4 71A103A (GCC→GCAyqgFpredicted Holliday junction resolvase
Reads supporting (aligned to +/- strand):  ref base C (29/6);  new base A (16/20);  total (45/26)
Fisher's exact test for biased strand distribution p-value = 1.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCT  >  minE/1980750‑1980862
                                                    |                                                            
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:1070048/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:1821940/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:204000/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:356246/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:2352614/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:2368813/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:2733771/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:2801677/69‑1 (MQ=255)
gTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCCTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:1211546/69‑1 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAAc                                           >  1:2043663/1‑69 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAAc                                           >  1:448685/1‑69 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTcc                                                        >  1:507221/1‑56 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:3313173/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:252705/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:3234851/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:2264464/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:215060/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:3128992/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:2034083/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:30019/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:40150/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:1861287/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:1757821/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:114192/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:1144864/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:1170746/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCg                                                       >  1:1961081/1‑57 (MQ=255)
   tCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTg                                              <  1:3118231/66‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATc                               >  1:1436135/1‑68 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              >  1:1999645/1‑69 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              >  1:1772531/1‑69 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              >  1:2697350/1‑69 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:817833/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:1258076/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:2607908/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:1380186/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:2300372/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:2192745/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:1538093/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:2123833/69‑1 (MQ=255)
                aaaGCTCCATGACGAGCGTCTTAGCACCGTGGAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:607001/69‑1 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:833434/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:3351242/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:876159/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:683999/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:1071063/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:2843291/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:1801789/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:2176445/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:2618930/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCgg                             >  1:2623164/1‑69 (MQ=255)
                 aaGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              >  1:3399457/1‑68 (MQ=255)
                      ccATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGgcgc                         >  1:2030533/1‑68 (MQ=255)
                      ccATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGgcgc                         >  1:43384/1‑68 (MQ=255)
                                    ttAGCACTGTGGAAGCCCGTTGCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAg          >  1:3346685/1‑69 (MQ=255)
                                    ttAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAg          >  1:3215472/1‑69 (MQ=255)
                                    ttAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAg          >  1:1944392/1‑69 (MQ=255)
                                    ttAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAg          >  1:796412/1‑69 (MQ=255)
                                    ttAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAg          >  1:853701/1‑69 (MQ=255)
                                         aCTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGa      <  1:3293470/68‑1 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  >  1:2981063/1‑69 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  >  1:1205892/1‑69 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  >  1:93313/1‑69 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  <  1:874713/69‑1 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  <  1:866591/69‑1 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  <  1:1370335/69‑1 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  >  1:140959/1‑69 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  <  1:1775332/69‑1 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  >  1:44527/1‑69 (MQ=255)
                                            gtgGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGActct  <  1:2619819/69‑1 (MQ=255)
                                              ggAGGCACGTTCCGGTCTGTTTGAACAGGGCGGCTATCg                              <  1:1439409/39‑1 (MQ=38)
                                                    |                                                            
GTTTCGGTGTTGAAGTAAAGCTCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCT  >  minE/1980750‑1980862

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