Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,028 G→A 54.9% A573V (GCA→GTA)  uvrA ← ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,811,0280GA54.9% 18.7 / 94.6 71A573V (GCA→GTA) uvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base G (25/7);  new base A (27/12);  total (52/19)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTCAGCGGC  >  minE/2810963‑2811089
                                                                 |                                                             
tGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGAc                                                           <  1:1840529/70‑1 (MQ=255)
   cGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCAc                                                        <  1:1049163/70‑1 (MQ=255)
      tCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCa                                                    >  1:1105284/1‑71 (MQ=255)
      tCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCa                                                    >  1:3482445/1‑71 (MQ=255)
        aTGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCAcc                                                  >  1:2684273/1‑71 (MQ=255)
         tGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCAccg                                                 <  1:2909782/71‑1 (MQ=255)
         tGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGAcc                                                          <  1:225661/62‑1 (MQ=255)
           tACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGc                                                              >  1:2651324/1‑56 (MQ=255)
           tACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgcc                                               <  1:3296725/71‑1 (MQ=255)
           tACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgcc                                               <  1:511513/71‑1 (MQ=255)
            aCTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgcc                                               <  1:783463/70‑1 (MQ=255)
             cTGCCCGGGCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCAccgccg                                              <  1:1714828/70‑1 (MQ=255)
                cccGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGaa                                          >  1:3091605/1‑71 (MQ=255)
                 ccGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGaa                                          >  1:408434/1‑70 (MQ=255)
                   ggTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAacac                                       <  1:167422/71‑1 (MQ=255)
                   ggTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAacac                                       <  1:1477054/71‑1 (MQ=255)
                     tCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCt                                     >  1:488330/1‑71 (MQ=255)
                     tCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCt                                     <  1:1095384/71‑1 (MQ=255)
                      cAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt                                                      >  1:1940305/1‑53 (MQ=255)
                      cAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt                                                      >  1:2488173/1‑53 (MQ=255)
                      cAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTg                                    >  1:2098093/1‑71 (MQ=255)
                      aaaCGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACAcc                                      >  1:852358/2‑69 (MQ=255)
                       aaCGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgc                                   >  1:1061951/1‑71 (MQ=255)
                        aCGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTgcg                                  >  1:205210/1‑71 (MQ=255)
                        aCGACTCCGGCACCGCCATAATCGCTTCAAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTgcg                                  >  1:3496016/1‑71 (MQ=255)
                         cGACTCCGGCACCGCCATAATCGCTTACAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTgcg                                  <  1:2119332/70‑1 (MQ=255)
                          gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAac                                         >  1:3070300/1‑62 (MQ=255)
                          gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGAc                                                           <  1:913364/44‑1 (MQ=255)
                          gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGAc                                                           <  1:732347/44‑1 (MQ=255)
                             tCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCg                              >  1:855989/1‑70 (MQ=255)
                             tCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCg                              >  1:3421776/1‑70 (MQ=255)
                             tCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGAc                                                           <  1:93966/41‑1 (MQ=255)
                             tCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGa                                           >  1:239750/1‑57 (MQ=255)
                             tCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCgg                             >  1:3101480/1‑71 (MQ=255)
                               cGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgc                                   >  1:285407/1‑63 (MQ=255)
                               cGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCgg                             <  1:2915442/69‑1 (MQ=255)
                               cGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTc                                                     >  1:2070475/1‑45 (MQ=255)
                               cGGCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGc                           >  1:234621/1‑71 (MQ=255)
                                ggCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCa                                                         >  1:3483597/1‑40 (MQ=255)
                                 gCACCGCCATAATCGCTTCCAGCGGACCTTGTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGcc                          >  1:3410679/1‑70 (MQ=255)
                                 gCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGcc                          >  1:3342375/1‑70 (MQ=255)
                                 gCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAac                                         >  1:2639996/1‑55 (MQ=255)
                                 gCACCGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTCACAGCCGTGAac                                         >  1:3698724/1‑55 (MQ=255)
                                     cGCCATAATCGCTTCCAGCGGACCTTCTACGACCACTTc                                                     <  1:2039100/39‑1 (MQ=255)
                                      gCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTc                    >  1:3680066/1‑71 (MQ=255)
                                       ccATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGaa                                          >  1:2569344/1‑48 (MQ=255)
                                       ccATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGaa                                          >  1:270076/1‑48 (MQ=255)
                                       ccATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTc                    >  1:1180594/1‑70 (MQ=255)
                                            aTCGCTTCCAGCGGACCTTCTGCGCCCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca               >  1:322035/1‑70 (MQ=255)
                                            aTCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGcc                                >  1:1165009/1‑53 (MQ=255)
                                            aTCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca               >  1:3116267/1‑70 (MQ=255)
                                            aTCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca               <  1:2831757/70‑1 (MQ=255)
                                            aTCGCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca               >  1:3487182/1‑70 (MQ=255)
                                              cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACAt            >  1:3073070/1‑71 (MQ=255)
                                              cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACAt            >  1:185811/1‑71 (MQ=255)
                                               gCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGaa                                          >  1:3524395/1‑40 (MQ=255)
                                               gCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGaa                                          >  1:1687958/1‑40 (MQ=255)
                                               gCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGaa                                          >  1:1691105/1‑40 (MQ=255)
                                               gCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGaa                                          >  1:3160824/1‑40 (MQ=255)
                                               gCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGaa                                          >  1:514976/1‑40 (MQ=255)
                                               gCTTCCAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATg           >  1:2705992/1‑71 (MQ=255)
                                                   ccAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGt         >  1:1783828/1‑69 (MQ=255)
                                                   ccAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGt         >  1:2762187/1‑69 (MQ=255)
                                                   ccAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGt         >  1:1708295/1‑69 (MQ=255)
                                                   ccAGCGGACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGt         >  1:3033727/1‑69 (MQ=255)
                                                        ggACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTCAgcggc  >  1:1883795/1‑71 (MQ=255)
                                                        ggACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTCAgcgg   >  1:1326201/1‑70 (MQ=255)
                                                        ggACCTTCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTCAgcgg   <  1:2939328/70‑1 (MQ=255)
                                                              tCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCg                   >  1:3526066/1‑48 (MQ=255)
                                                              tCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCg                   >  1:2974053/1‑48 (MQ=255)
                                                              tCTACGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTCAgcgg   >  1:1236172/1‑64 (MQ=255)
                                                                 |                                                             
TGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTCAGCGGC  >  minE/2810963‑2811089

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: