Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,039,558 C→T 38.7% P115P (CCC→CCT nth → DNA glycosylase and apyrimidinic (AP) lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,039,5580CT38.7% 33.5 / 42.5 67P115P (CCC→CCTnthDNA glycosylase and apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base C (13/28);  new base T (14/12);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 8.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCT  >  minE/1039491‑1039621
                                                                   |                                                               
gCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTAAGCCCTGCCCgg                                                               <  1:558403/69‑1 (MQ=255)
   gTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg                                                               <  1:847243/67‑1 (MQ=255)
   gTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg                                                               <  1:2297208/67‑1 (MQ=255)
   gTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg                                                               <  1:296488/67‑1 (MQ=255)
      tCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAg                                                          >  1:2877045/1‑69 (MQ=255)
      tCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAg                                                          >  1:3272381/1‑69 (MQ=255)
          gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                       >  1:1004949/1‑68 (MQ=255)
          gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                       >  1:95574/1‑68 (MQ=255)
          gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                       >  1:1113102/1‑68 (MQ=255)
          gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                       >  1:2731509/1‑68 (MQ=255)
          gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                       >  1:3243870/1‑68 (MQ=255)
             ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa                                                   <  1:992508/69‑1 (MQ=255)
             ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa                                                   <  1:2365172/69‑1 (MQ=255)
             ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa                                                   <  1:2361242/69‑1 (MQ=255)
             ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa                                                   <  1:22935/69‑1 (MQ=255)
             ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa                                                   <  1:2244559/69‑1 (MQ=255)
                      tAATGGCGAGGTTCCTGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCa                                           >  1:791961/1‑68 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACa                                               >  1:1672572/1‑64 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2598168/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:1804106/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2653480/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2498331/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2469014/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:913401/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:1537942/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2955234/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:1608366/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2132180/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:3093530/1‑69 (MQ=255)
                      tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa                                          >  1:2013940/1‑69 (MQ=255)
                        aTGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg                                        >  1:1443645/1‑69 (MQ=255)
                        aTGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg                                        >  1:417172/1‑69 (MQ=255)
                        aTGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg                                        >  1:2537885/1‑69 (MQ=255)
                          ggCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc                                                <  1:312947/59‑1 (MQ=255)
                          ggCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc                                                <  1:3426769/59‑1 (MQ=255)
                          ggCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc                                                <  1:1811971/59‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:513989/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:86150/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:2457228/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:2525071/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:237546/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:2319572/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:22482/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:1264407/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:1687090/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:1131766/69‑1 (MQ=255)
                             gAGGTTCCGGAAGATCGTGCTGCCCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:981471/69‑1 (MQ=255)
                              aGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa                                   <  1:1917544/68‑1 (MQ=255)
                                gTTCCGGGAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:2453568/69‑1 (MQ=255)
                                gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:806535/69‑1 (MQ=255)
                                gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:1518378/69‑1 (MQ=255)
                                gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:1095858/69‑1 (MQ=255)
                                gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:18220/69‑1 (MQ=255)
                                gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:3014795/69‑1 (MQ=255)
                                 ttCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:871516/68‑1 (MQ=255)
                                 ttCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                >  1:316670/1‑68 (MQ=255)
                                 ttCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                <  1:2025755/68‑1 (MQ=255)
                                   ccGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAac                             >  1:3419881/1‑69 (MQ=255)
                                                ctgcGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCt                >  1:3201873/1‑69 (MQ=255)
                                                  gcgcTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTg               <  1:1090486/68‑1 (MQ=255)
                                                       tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc         <  1:1544071/69‑1 (MQ=255)
                                                       tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc         <  1:159355/69‑1 (MQ=255)
                                                       tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc         <  1:353748/69‑1 (MQ=255)
                                                       tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc         <  1:176119/69‑1 (MQ=255)
                                                       tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc         <  1:2652711/69‑1 (MQ=255)
                                                              cTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCt  <  1:1757263/69‑1 (MQ=255)
                                                               tGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCt  <  1:2030818/68‑1 (MQ=255)
                                                                   |                                                               
GCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCT  >  minE/1039491‑1039621

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: