Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,250,792 A→C 80.4% L167L (CTT→CTG ruvC ← component of RuvABC resolvasome, endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,250,7920AC80.4% 217.1 / 27.8 86L167L (CTT→CTGruvCcomponent of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base A (11/6);  new base C (38/31);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCG  >  minE/1250727‑1250859
                                                                 |                                                                   
gggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt                                                                >  1:1656491/1‑71 (MQ=255)
 ggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg                                                                   >  1:4034332/1‑67 (MQ=255)
     tGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCa                                                                     <  1:926739/61‑1 (MQ=255)
      ggAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcg                                                          >  1:2847345/1‑71 (MQ=255)
       gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc                                                               <  1:2310956/65‑1 (MQ=255)
       gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc                                                               <  1:3858193/65‑1 (MQ=255)
       gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc                                                               <  1:3292639/65‑1 (MQ=255)
       gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc                                                               <  1:5232315/65‑1 (MQ=255)
       gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcgc                                                         >  1:776429/1‑71 (MQ=255)
        aaaCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcgc                                                         <  1:3940179/70‑1 (MQ=255)
          aCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCa                                                                     <  1:5542831/56‑1 (MQ=255)
            gCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTc                                                    >  1:2783068/1‑71 (MQ=255)
             ccTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCg                                                   >  1:3366392/1‑71 (MQ=255)
              cTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGc                                                  <  1:2606734/71‑1 (MQ=255)
              cTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGc                                                  >  1:1763296/1‑71 (MQ=255)
              cTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGc                                                  >  1:1584978/1‑71 (MQ=255)
               tCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCa                                                                     <  1:2921614/51‑1 (MQ=255)
                 aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCa                                               >  1:5379213/1‑71 (MQ=255)
                 aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCa                                               >  1:1588906/1‑71 (MQ=255)
                 aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCa                                                              >  1:1139307/1‑56 (MQ=255)
                  gCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGcc                                                           >  1:593392/1‑58 (MQ=255)
                  gCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCTCTCAt                                              >  1:2521923/1‑71 (MQ=255)
                   ccGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCAt                                              <  1:2621678/70‑1 (MQ=255)
                      gAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg                                                                  >  1:3654471/1‑47 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCTGCCGCGATTCGCTCATCTGCATCGCa                                    <  1:3995490/71‑1 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                      <  1:734720/69‑1 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                    <  1:3620735/71‑1 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                    <  1:5096119/71‑1 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                    >  1:3978096/1‑71 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                    <  1:3062776/71‑1 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                    <  1:1837034/71‑1 (MQ=255)
                            aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                    >  1:4057136/1‑71 (MQ=255)
                             ccGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGAtt                                                     >  1:5537674/1‑53 (MQ=255)
                              cGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAtt                                  >  1:3844930/1‑71 (MQ=255)
                                aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGt                                                                 >  1:3054566/1‑38 (MQ=255)
                                    ggTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGAc                             >  1:3307743/1‑70 (MQ=255)
                                    ggTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGAc                             >  1:2327451/1‑70 (MQ=255)
                                      tataACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcgcg                                                        <  1:4703544/41‑1 (MQ=255)
                                      tataACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcgcg                                                        <  1:2878748/41‑1 (MQ=255)
                                       ataACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCAt                                              >  1:3117652/1‑50 (MQ=255)
                                       ataACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCAt                                              >  1:2490159/1‑50 (MQ=255)
                                       ataACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCAt                                              >  1:2192877/1‑50 (MQ=255)
                                       ataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                              >  1:5252333/1‑66 (MQ=255)
                                       ataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                              >  1:1742769/1‑66 (MQ=255)
                                       ataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAAc                         >  1:4943233/1‑71 (MQ=255)
                                        taACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGcc                                                           <  1:4883637/36‑1 (MQ=255)
                                        taACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACg                        >  1:2442969/1‑71 (MQ=255)
                                         aaCTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGt                       >  1:2458016/1‑71 (MQ=255)
                                         aaCTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGt                       >  1:905482/1‑71 (MQ=255)
                                         aaCTTAACACAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGt                       >  1:3890810/1‑71 (MQ=255)
                                          aCTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGa                                                       >  1:477952/1‑38 (MQ=255)
                                          aCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGAt                                                      >  1:5252027/1‑39 (MQ=255)
                                          aCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCa                                               >  1:2012451/1‑46 (MQ=255)
                                           cTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     <  1:2120054/71‑1 (MQ=255)
                                           cTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     <  1:3208189/71‑1 (MQ=255)
                                           cTTAACGCAGTCGCCCTCTCGCCAGGTGCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     <  1:1684484/71‑1 (MQ=255)
                                           cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     <  1:4063503/71‑1 (MQ=255)
                                           cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     <  1:3337786/71‑1 (MQ=255)
                                           cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     <  1:4846725/71‑1 (MQ=255)
                                            ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATc                                             >  1:2089826/1‑46 (MQ=255)
                                            ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATc                                             >  1:760888/1‑46 (MQ=255)
                                            ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATc                                             >  1:5164677/1‑46 (MQ=255)
                                              aaCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGc                    <  1:2961143/69‑1 (MQ=255)
                                                cGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                              >  1:4125567/1‑57 (MQ=255)
                                                 gCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTgg               >  1:1321803/1‑71 (MQ=255)
                                                  cAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACg                        <  1:854131/61‑1 (MQ=255)
                                                  cAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACg                        <  1:2937428/61‑1 (MQ=255)
                                                  cAGTCGCCCTCTCGCCAGGTCCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACg                        <  1:4306517/61‑1 (MQ=255)
                                                    gTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                <  1:1768167/51‑1 (MQ=255)
                                                    gTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTg                      >  1:1920164/1‑61 (MQ=255)
                                                    gTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTg                      >  1:2940575/1‑61 (MQ=255)
                                                    gTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTg                      >  1:3316480/1‑61 (MQ=255)
                                                    gTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGt                 <  1:493833/66‑1 (MQ=255)
                                                         ccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                              >  1:4146381/1‑48 (MQ=255)
                                                         ccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                              >  1:2300581/1‑48 (MQ=255)
                                                         ccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                              >  1:5046921/1‑48 (MQ=255)
                                                         ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTACCGTGGCAGTGGGTGATAg        >  1:857045/1‑70 (MQ=255)
                                                         ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAg        <  1:2038286/70‑1 (MQ=255)
                                                          ctctCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCa      <  1:83791/71‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTc                                 <  1:3947265/43‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTc                                 <  1:4308801/43‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                <  1:879897/44‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                     >  1:2410675/1‑55 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTg            <  1:3427451/64‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTg            <  1:3884033/64‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTg            <  1:2008479/64‑1 (MQ=255)
                                                           tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCaa     >  1:676366/1‑71 (MQ=255)
                                                              cGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCg  >  1:3392074/1‑71 (MQ=255)
                                                               gCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCg  <  1:1763235/70‑1 (MQ=255)
                                                                 aaGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGt                       <  1:3947955/47‑1 (MQ=255)
                                                                 |                                                                   
GGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCG  >  minE/1250727‑1250859

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: