Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,249,224 C→T 45.2% D17N (GAT→AAT)  ruvB ← ATP‑dependent DNA helicase, component of RuvABC resolvasome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,249,2240CT45.2% 24.9 / 99.3 73D17N (GAT→AAT) ruvBATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base C (26/14);  new base T (12/21);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

AACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGC  >  minE/1249158‑1249289
                                                                  |                                                                 
aacgaacCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTcc                                                               >  1:1146627/1‑71 (MQ=255)
    aacCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCTGCa                                                           <  1:703556/71‑1 (MQ=255)
    aacCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCa                                                           <  1:830925/71‑1 (MQ=255)
     acCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCa                                                           <  1:4163913/70‑1 (MQ=255)
       cTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAg                                                        <  1:3009676/71‑1 (MQ=255)
       cTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAg                                                        <  1:1289046/71‑1 (MQ=255)
       cTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAg                                                        <  1:1592501/71‑1 (MQ=255)
          cGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGt                                                       <  1:4225050/69‑1 (MQ=255)
          cGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGt                                                       <  1:3170994/69‑1 (MQ=255)
          cGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGt                                                       <  1:2105487/69‑1 (MQ=255)
               gACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCg                                                >  1:2915738/1‑71 (MQ=255)
                  cAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAg                                                        >  1:563035/1‑60 (MQ=255)
                  cAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCgg                                               >  1:475323/1‑69 (MQ=255)
                   aaCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGt                                                       <  1:1204116/60‑1 (MQ=255)
                   aaCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCgg                                               >  1:834757/1‑68 (MQ=255)
                   aaCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCgg                                               >  1:395449/1‑68 (MQ=255)
                   aaCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCgg                                               >  1:1249954/1‑68 (MQ=255)
                    aCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAg                                            >  1:3181956/1‑70 (MQ=255)
                    aCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAg                                            >  1:2299415/1‑70 (MQ=255)
                     cATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTgg                                                     >  1:4031286/1‑60 (MQ=255)
                     cATACTATTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCaaa                                                         >  1:2387989/1‑56 (MQ=255)
                      aTACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGaaa                                         >  1:3056356/1‑71 (MQ=255)
                      aTACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGaaa                                         <  1:787148/71‑1 (MQ=255)
                       tACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGaaa                                         >  1:981451/1‑70 (MQ=255)
                          tctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGaaa                                         >  1:3361186/1‑67 (MQ=255)
                          tctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAAt                                        >  1:1776411/1‑68 (MQ=255)
                          tctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAAt                                        >  1:2614867/1‑68 (MQ=255)
                          tctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGG‑‑GAAATCAGAc                                   >  1:3568926/1‑71 (MQ=255)
                          tctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGCCAGAAAt                                        >  1:377011/1‑68 (MQ=255)
                           ctTCCAGTAATTTGGGGCGAATGGCGAGATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGa                                    >  1:1558141/1‑71 (MQ=255)
                            ttCCAGTAATTTGGTGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGt                                                    <  1:1266732/54‑1 (MQ=255)
                            ttCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGt                                                    <  1:3579695/54‑1 (MQ=255)
                            ttCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGt                                                    <  1:3635833/54‑1 (MQ=255)
                              ccAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:20616/1‑42 (MQ=255)
                              ccAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACgg                                 >  1:1246490/1‑71 (MQ=255)
                              ccAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACgg                                 >  1:1256308/1‑71 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:578143/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:3448699/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:25568/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:401346/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:1669431/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:2436797/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:2307278/41‑1 (MQ=255)
                                  tAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCa                                                           <  1:4113650/41‑1 (MQ=255)
                                    aTTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:321380/1‑36 (MQ=255)
                                    aTTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:1639508/1‑36 (MQ=255)
                                    aTTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:190319/1‑36 (MQ=255)
                                    aTTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:1920083/1‑36 (MQ=255)
                                    aTTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:2481764/1‑36 (MQ=255)
                                    aTTTGGGGCGAATGGCGCGATCTGCTACATCTTCCg                                                              >  1:2774765/1‑36 (MQ=255)
                                     tttGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGa                                    >  1:3520759/1‑61 (MQ=255)
                                     tttGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGa                                    >  1:1372634/1‑61 (MQ=255)
                                     tttGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGa                                    >  1:2135945/1‑61 (MQ=255)
                                     tttGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGa                                    >  1:2501970/1‑61 (MQ=255)
                                      ttGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCtt                          >  1:3395351/1‑70 (MQ=255)
                                            cGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAatcat                   >  1:3579350/1‑71 (MQ=255)
                                            cGAATGGCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAatca                    >  1:1560675/1‑70 (MQ=255)
                                               aTGGCGCGATCTGCTACATTTTCCGGCAAAGTGGt                                                    <  1:1040488/35‑1 (MQ=255)
                                               aTGGCGCGATCTGCTACATTTTCCGGCAAAGTGGt                                                    <  1:954882/35‑1 (MQ=255)
                                               aTGGCGCGATCTGCTACATTTTCCGGCAAAGTGGt                                                    <  1:959857/35‑1 (MQ=255)
                                               aTGGCGCGATCTGCTACATTTTCCGGCAAAGTGGt                                                    <  1:1565779/35‑1 (MQ=255)
                                                tGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGaa                                          <  1:1396016/44‑1 (MQ=255)
                                                tGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCtt               <  1:116473/71‑1 (MQ=255)
                                                 ggCGCGATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCtt                          >  1:4186057/1‑59 (MQ=255)
                                                     cgATCTGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCt          >  1:954661/1‑71 (MQ=255)
                                                          tGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCa        <  1:1552159/68‑1 (MQ=255)
                                                          tGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCa        <  1:3931981/68‑1 (MQ=255)
                                                          tGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCa        <  1:3616069/68‑1 (MQ=255)
                                                          tGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCa        <  1:3104118/68‑1 (MQ=255)
                                                          tGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCa        <  1:152259/68‑1 (MQ=255)
                                                          tGCTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCa        <  1:3497568/68‑1 (MQ=255)
                                                            cTACATTTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAAcg   >  1:2800280/1‑71 (MQ=255)
                                                               cATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAAcgc  <  1:1568004/69‑1 (MQ=255)
                                                                  |                                                                 
AACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGC  >  minE/1249158‑1249289

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: