Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,401,943 A→G 31.2% N128N (AAT→AAC xerC ← site‑specific tyrosine recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,401,9430AG31.2% 68.2 / 39.5 80N128N (AAT→AACxerCsite‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base A (7/48);  new base G (2/23);  total (9/71)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AAAGACGCAGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCA  >  minE/2401879‑2401998
                                                                |                                                       
aaaGACGCAGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGa                                                     <  1:2727970/69‑1 (MQ=255)
aaaGACGCAGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGa                                                     <  1:2860462/69‑1 (MQ=255)
 aaGACGCAGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGa                                                     <  1:2899252/68‑1 (MQ=255)
        aGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAata                                             <  1:2066586/69‑1 (MQ=255)
        aGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAata                                             <  1:575009/69‑1 (MQ=255)
        aGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAata                                             <  1:716583/69‑1 (MQ=255)
         gACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGa                                                     <  1:264465/60‑1 (MQ=255)
           cccGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATcc                                          <  1:873463/69‑1 (MQ=255)
           cccGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATcc                                          <  1:3397887/69‑1 (MQ=255)
           cccGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAAGGGATCATTGATATCAATATcc                                          <  1:2704612/69‑1 (MQ=255)
                    tACATCACTTCCAGCATTGCACGGTCGCGTAGAGCGAGGGGATCATTGATATc                                                 >  1:735111/1‑53 (MQ=255)
                    tACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATc                                                 >  1:1389035/1‑53 (MQ=255)
                     aCATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCa                                                         >  1:860320/1‑44 (MQ=255)
                         cACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:558451/69‑1 (MQ=255)
                         cACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:3319117/69‑1 (MQ=255)
                         cACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:3330411/69‑1 (MQ=255)
                         cACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:117692/69‑1 (MQ=255)
                         cACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:1196829/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCATCAGCCGATTCata                           <  1:3113142/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1171565/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:255689/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:263712/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1031034/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1132720/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1141571/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1161094/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1986483/69‑1 (MQ=255)
                          aCTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1430639/69‑1 (MQ=255)
                          aCTTCCAGCAATGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1335783/69‑1 (MQ=255)
                           cTTCCAGCATTGCACGGTCGCGTACAGCGAGGTGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:116398/67‑1 (MQ=255)
                           cTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTCatat                          >  1:2045813/1‑69 (MQ=255)
                           cTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTCatat                          >  1:2587915/1‑69 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:3080564/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:421105/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:3162830/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:2676499/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:233790/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:519839/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:3435139/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:2400859/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:2363378/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:16206/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:2232678/68‑1 (MQ=255)
                                aGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATAtcgtc                      <  1:164955/68‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2041331/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:1064346/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:524361/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:1464444/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:517273/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:1471010/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:1575708/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:3283104/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:3239048/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:3238917/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:1739326/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:669855/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:625571/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:3016737/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2976418/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2970596/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2082523/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2846833/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2258056/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2292233/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:58655/54‑1 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCATTGATATc                                                 >  1:22897/1‑35 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGAtt                               >  1:1692787/1‑53 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGAtt                               >  1:166192/1‑53 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGAtt                               >  1:2589829/1‑53 (MQ=255)
                                      gCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATcgtcgacgtc                >  1:718295/1‑68 (MQ=255)
                                      gCACGGTCGCGCACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:2435765/54‑1 (MQ=39)
                                       cACGGTCGCGTACAGCGAGGGGATCGTTGATATCAATATCCAGCAGCCGATTc                              <  1:3037960/53‑1 (MQ=255)
                                            tCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:1425852/50‑1 (MQ=255)
                                                gTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCat                            <  1:201315/46‑1 (MQ=255)
                                                   cAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCa  <  1:721251/69‑1 (MQ=255)
                                                   cAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCa  <  1:90262/69‑1 (MQ=255)
                                                   cAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCa  <  1:2570851/69‑1 (MQ=255)
                                                   cAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCa  <  1:2045645/69‑1 (MQ=255)
                                                   cAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCa  <  1:2033420/69‑1 (MQ=255)
                                                   cAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCa  <  1:1435830/69‑1 (MQ=255)
                                                            gATCATTGATATCAATATCCAGCAGCCGATTCata                           <  1:1104525/35‑1 (MQ=255)
                                                                |                                                       
AAAGACGCAGACCCGCGCCGTACATCACTTCCAGCATTGCACGGTCGCGTACAGCGAGGGGATCATTGATATCAATATCCAGCAGCCGATTCATATCGTCGACGTCGATGTTTTTCGGCA  >  minE/2401879‑2401998

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: