Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,401,662 T→C 8.5% K222R (AAA→AGA)  xerC ← site‑specific tyrosine recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,401,6620TC 8.5% 258.6 / 12.6 94K222R (AAA→AGA) xerCsite‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base T (38/48);  new base C (4/4);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 103 total aligned reads displayed.

CGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCAAA  >  minE/2401592‑2401732
                                                                      |                                                                      
cGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtt                                                                        <  1:1181895/71‑1 (MQ=255)
 gAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                       <  1:590890/71‑1 (MQ=255)
 gAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                       <  1:123110/71‑1 (MQ=255)
   gTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCt                                                                     <  1:1215433/71‑1 (MQ=255)
   gTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCt                                                                     <  1:952014/71‑1 (MQ=255)
    tGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTc                                                                      <  1:758453/69‑1 (MQ=255)
    tGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTc                                                                      <  1:3646005/69‑1 (MQ=255)
    tGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg                                                                    >  1:2754664/1‑71 (MQ=255)
      aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca                                                                  <  1:3336985/71‑1 (MQ=255)
      aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca                                                                  >  1:834341/1‑71 (MQ=255)
       cGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca                                                                  >  1:30647/1‑70 (MQ=255)
       cGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca                                                                  >  1:3406706/1‑70 (MQ=255)
       cGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca                                                                  >  1:3709446/1‑70 (MQ=255)
         tAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGTCAAAGCGTTTCTGc                                                                   <  1:141802/67‑1 (MQ=255)
         tAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGc                                                                   <  1:120435/67‑1 (MQ=255)
         tAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAt                                                               <  1:2250100/71‑1 (MQ=255)
         tAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAt                                                               <  1:1271621/71‑1 (MQ=255)
         tAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAt                                                               >  1:1263412/1‑71 (MQ=255)
           aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                             <  1:583256/71‑1 (MQ=255)
           aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                             >  1:2274618/1‑71 (MQ=255)
           aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTCCACATTa                                                             <  1:3944849/71‑1 (MQ=255)
            tttATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCt                                                                     <  1:962201/62‑1 (MQ=255)
            tttATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAc                                                            >  1:1699001/1‑71 (MQ=255)
             ttATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcg                                                           >  1:478246/1‑71 (MQ=255)
             ttATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcg                                                           >  1:3314535/1‑71 (MQ=255)
              tATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcg                                                           >  1:2548818/1‑70 (MQ=255)
              tATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcg                                                           <  1:2511685/70‑1 (MQ=255)
                   ggATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                       >  1:3243339/1‑69 (MQ=255)
                     aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTc                                                                      <  1:431360/52‑1 (MQ=255)
                      tGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGAt                                                  <  1:1230371/71‑1 (MQ=255)
                        aaCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATAc                                                >  1:2802609/1‑71 (MQ=255)
                         aCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACg                                               <  1:2008093/71‑1 (MQ=255)
                         aCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACg                                               <  1:1778964/71‑1 (MQ=255)
                          cGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc                                              <  1:3615256/71‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTc                                                                      >  1:2416798/1‑46 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             <  1:1416217/71‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             >  1:2081537/1‑71 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             <  1:1411625/71‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             <  1:1415218/71‑1 (MQ=255)
                            tGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCtt                                            >  1:3938043/1‑71 (MQ=255)
                              attattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATAc                                                <  1:652505/65‑1 (MQ=255)
                               ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                       >  1:79570/1‑41 (MQ=255)
                               ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                       >  1:382704/1‑57 (MQ=255)
                               ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                       >  1:3622360/1‑57 (MQ=255)
                               ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                       >  1:1954101/1‑57 (MQ=255)
                                tattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGcc                                         >  1:2507119/1‑70 (MQ=255)
                                 attCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGccc                                        <  1:1387979/70‑1 (MQ=255)
                                  ttCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAg                                      <  1:3535730/71‑1 (MQ=255)
                                  ttCAGACCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCa                                       <  1:1104031/70‑1 (MQ=255)
                                      gCCCTTGTTTTATGCCCGATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                                  <  1:3698493/71‑1 (MQ=255)
                                      gCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                                  <  1:2674888/71‑1 (MQ=255)
                                      gCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                                  <  1:1361550/71‑1 (MQ=255)
                                        ccTTGTTTTATGCCCCATTCTGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                <  1:3494043/71‑1 (MQ=255)
                                         cTTGTTTTTTGCCCCATTCGGCAAAGCGTCTCTGCTCATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                <  1:3498622/70‑1 (MQ=255)
                                          ttGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                >  1:72254/1‑69 (MQ=255)
                                           tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                >  1:1460411/1‑68 (MQ=255)
                                           tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                >  1:548173/1‑68 (MQ=255)
                                           tGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                >  1:2409410/1‑68 (MQ=255)
                                            gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAg                             <  1:2933775/70‑1 (MQ=255)
                                            gTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                               <  1:2776764/36‑1 (MQ=255)
                                            gTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAg                             >  1:63198/1‑70 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             >  1:3216884/1‑51 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             >  1:2470430/1‑51 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                             >  1:2052770/1‑51 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaa                          <  1:2817506/70‑1 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaa                          >  1:1980214/1‑70 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaa                          <  1:229332/70‑1 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGAt                                                  >  1:957688/1‑46 (MQ=255)
                                                 aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaaa                        <  1:3952686/70‑1 (MQ=255)
                                                  tGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAg                       <  1:425150/70‑1 (MQ=255)
                                                    ccccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcgc                    <  1:1984254/71‑1 (MQ=255)
                                                     cccATTCGGCAAAGCGTTTCTGCACATTACGCGCGga                                                     >  1:3298146/1‑37 (MQ=255)
                                                     cccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGcgcg                   <  1:1587078/71‑1 (MQ=255)
                                                     cccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGcgcg                   <  1:1922289/71‑1 (MQ=255)
                                                       cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc                 >  1:3939226/1‑71 (MQ=255)
                                                       cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgt                  >  1:3517588/1‑70 (MQ=255)
                                                         ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcgc                    <  1:83120/66‑1 (MQ=255)
                                                         ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt               >  1:311938/1‑71 (MQ=255)
                                                         ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcg                >  1:3438681/1‑70 (MQ=255)
                                                           cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATa                                                 >  1:780626/1‑35 (MQ=255)
                                                           cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt             >  1:1956801/1‑71 (MQ=255)
                                                            ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGa                            <  1:1429110/55‑1 (MQ=255)
                                                            ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCtt            >  1:1558097/1‑71 (MQ=255)
                                                            ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt             <  1:2199791/70‑1 (MQ=255)
                                                             gCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTc           <  1:2816940/71‑1 (MQ=255)
                                                             gCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTc           <  1:2664988/71‑1 (MQ=255)
                                                              cAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTc           <  1:1402445/70‑1 (MQ=255)
                                                              cAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCg          <  1:188144/71‑1 (MQ=255)
                                                              cAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCg          <  1:1481547/71‑1 (MQ=255)
                                                                aaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCgtcgtc              <  1:487685/65‑1 (MQ=255)
                                                                aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                                  >  1:1201771/1‑45 (MQ=255)
                                                                 aGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTa       >  1:1124252/1‑71 (MQ=255)
                                                                  gCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTa       <  1:880803/70‑1 (MQ=255)
                                                                   cGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTAcc     <  1:2796240/71‑1 (MQ=255)
                                                                    gTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                                >  1:1351684/1‑43 (MQ=255)
                                                                     tttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGa                            <  1:3130797/46‑1 (MQ=255)
                                                                     tttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgc                      >  1:379727/1‑52 (MQ=255)
                                                                     tttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCaa   <  1:493553/71‑1 (MQ=255)
                                                                      ttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCaaa  >  1:2985391/1‑71 (MQ=255)
                                                                      ttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGAGTCGTCTTCGCTACCaa   >  1:760864/1‑70 (MQ=255)
                                                                      |                                                                      
CGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCAAA  >  minE/2401592‑2401732

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: