Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,250,792 A→C 10.7% 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,7920AC10.7% 199.6 / 21.6 87L167L (CTT→CTGruvCcomponent of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base A (43/34);  new base C (3/6);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCC  >  minE/1250728‑1250861
                                                                |                                                                     
ggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc                                                                 <  1:4107577/71‑1 (MQ=255)
ggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTc                                                                 <  1:111008/71‑1 (MQ=255)
   gTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCa                                                                       <  1:1768415/62‑1 (MQ=255)
   gTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCa                                                                       <  1:159332/62‑1 (MQ=255)
      gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAg                                                               >  1:1192865/1‑67 (MQ=255)
      gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCcgc                                                           >  1:3227745/1‑71 (MQ=255)
      gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCcgc                                                           >  1:741962/1‑71 (MQ=255)
      gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCcgc                                                           >  1:3736001/1‑71 (MQ=255)
      gAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCcgc                                                           >  1:1569508/1‑71 (MQ=255)
        aaCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCcgcg                                                          <  1:2672223/70‑1 (MQ=255)
        aaCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGa                                                         >  1:1997451/1‑71 (MQ=255)
         aCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGAt                                                        >  1:3257605/1‑71 (MQ=255)
           gCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGt                                                                   <  1:617147/58‑1 (MQ=255)
           gCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTc                                                      >  1:3966616/1‑71 (MQ=255)
           gCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTc                                                      >  1:4124444/1‑71 (MQ=255)
           gCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTc                                                      >  1:1709582/1‑71 (MQ=255)
              tCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGc                                                    >  1:2727059/1‑70 (MQ=255)
               cAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCa                                                                       <  1:2189709/50‑1 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGAtt                                                       >  1:2556503/1‑65 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTc                                                  <  1:3749612/70‑1 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTc                                                  <  1:3744711/70‑1 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTc                                                  <  1:588068/70‑1 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTc                                                  <  1:4046693/70‑1 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCa                                                 <  1:1077414/71‑1 (MQ=255)
                aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCCAGGTTCAGCCGCGATTCGCTCATc                                               >  1:1060061/1‑71 (MQ=255)
                   cGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCt                                              >  1:1598028/1‑71 (MQ=255)
                     gAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCt                                              >  1:1400935/1‑69 (MQ=255)
                     gAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTg                                             >  1:2557884/1‑70 (MQ=255)
                       aCTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCAt                                          >  1:787920/1‑71 (MQ=255)
                       aCTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCAt                                          >  1:1375457/1‑71 (MQ=255)
                        cTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATc                                         >  1:2381622/1‑71 (MQ=255)
                          gACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTc                                                                 >  1:3812072/1‑45 (MQ=255)
                           aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                        <  1:265575/69‑1 (MQ=255)
                           aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                        <  1:2403607/69‑1 (MQ=255)
                           aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                      >  1:3853015/1‑71 (MQ=255)
                            ccGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                        <  1:515561/68‑1 (MQ=255)
                            ccGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAt                                     <  1:1112308/71‑1 (MQ=255)
                               aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGa                                                         <  1:3170835/48‑1 (MQ=255)
                               aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                  <  1:1453051/71‑1 (MQ=255)
                               aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                  <  1:2329219/71‑1 (MQ=255)
                               aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                  <  1:2003373/71‑1 (MQ=255)
                               aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                  <  1:3141658/71‑1 (MQ=255)
                                    gTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCa                                                 <  1:669357/51‑1 (MQ=255)
                                     tataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCt                                                   <  1:1018588/48‑1 (MQ=255)
                                     tataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCt                                                   <  1:2951190/48‑1 (MQ=255)
                                     tataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                      <  1:984760/61‑1 (MQ=255)
                                       taACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGc                                                    >  1:585727/1‑45 (MQ=255)
                                       taACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACt                              <  1:2644459/67‑1 (MQ=255)
                                       taACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAAc                           >  1:703777/1‑70 (MQ=255)
                                        aaCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGc                                                    <  1:4081060/44‑1 (MQ=255)
                                        aaCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTc                                   >  1:71207/1‑61 (MQ=255)
                                         aCTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTg                        >  1:3720310/1‑71 (MQ=255)
                                         aCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCa                                           >  1:3411022/1‑52 (MQ=255)
                                         aCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCa                                           >  1:141833/1‑52 (MQ=255)
                                         aCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGt                         <  1:4041046/70‑1 (MQ=255)
                                         aCTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTg                        >  1:1620406/1‑71 (MQ=255)
                                          cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGc                                                    >  1:1892124/1‑42 (MQ=255)
                                          cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGc                                                    >  1:4033553/1‑42 (MQ=255)
                                          cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGc                                                    >  1:2488096/1‑42 (MQ=255)
                                          cTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCt                                                   >  1:3281834/1‑43 (MQ=255)
                                            tAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                       >  1:1194775/1‑69 (MQ=255)
                                             accGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGc                      <  1:3713282/67‑1 (MQ=255)
                                              aCGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAtt                                    <  1:4019884/54‑1 (MQ=255)
                                              aCGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAtt                                    <  1:1273242/54‑1 (MQ=255)
                                              aCGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGt                   <  1:2491451/71‑1 (MQ=255)
                                              aCGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGt                   <  1:635735/71‑1 (MQ=255)
                                               cGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGc                                       <  1:2363333/50‑1 (MQ=255)
                                                  aGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTc                                                  >  1:2384880/1‑36 (MQ=255)
                                                   gTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAgtgtgt               >  1:2199626/1‑70 (MQ=255)
                                                    tCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGa             <  1:796711/71‑1 (MQ=255)
                                                     cGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGa             >  1:2347600/1‑70 (MQ=255)
                                                     cGCCCTCTCGCAAGGTTCAGCCGAGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGc                      <  1:3522616/61‑1 (MQ=255)
                                                      gCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGa                                >  1:428815/1‑50 (MQ=255)
                                                        ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGc                                            <  1:2773746/36‑1 (MQ=255)
                                                        ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                        >  1:2984584/1‑40 (MQ=255)
                                                        ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                        >  1:2731396/1‑40 (MQ=255)
                                                        ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAtt                                    >  1:3673434/1‑44 (MQ=255)
                                                        ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGc         <  1:3424438/71‑1 (MQ=255)
                                                         ctctCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCa        <  1:1174511/71‑1 (MQ=255)
                                                          tctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAAc                           <  1:3787786/51‑1 (MQ=255)
                                                          tctcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                  <  1:1631965/44‑1 (MQ=255)
                                                           ctcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAtt                                    >  1:120664/1‑41 (MQ=255)
                                                            tcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATc     >  1:1682510/1‑71 (MQ=255)
                                                            tcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAAt      >  1:570804/1‑70 (MQ=255)
                                                            tcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAAt      >  1:269434/1‑70 (MQ=255)
                                                            tcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAAt      >  1:3059097/1‑70 (MQ=255)
                                                            tcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATc     >  1:185305/1‑71 (MQ=255)
                                                            tcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATc     >  1:2618886/1‑71 (MQ=255)
                                                              gCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGc   <  1:4038014/71‑1 (MQ=255)
                                                               cAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGcc  <  1:929415/71‑1 (MQ=255)
                                                                |                                                                     
GGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCC  >  minE/1250728‑1250861

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: