Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,639,730 T→C 35.7% K222R (AAA→AGA)  xerC ← site‑specific tyrosine recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,639,7300TC35.7% 66.1 / 73.4 70K222R (AAA→AGA) xerCsite‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base T (17/28);  new base C (11/14);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

AGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTA  >  W3110S.gb/3639662‑3639795
                                                                    |                                                                 
aGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtctc                                                                 <  1:862325/71‑1 (MQ=255)
  tGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg                                                               >  1:1824898/1‑71 (MQ=255)
   gACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg                                                               >  1:3446331/1‑70 (MQ=255)
    aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGca                                                             >  1:860678/1‑71 (MQ=255)
    aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGca                                                             >  1:4607190/1‑71 (MQ=255)
    aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGca                                                             >  1:3940503/1‑71 (MQ=255)
         atatATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTa                                                        >  1:2642485/4‑71 (MQ=255)
         aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                  >  1:4501243/1‑61 (MQ=255)
          tttATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAc                                                       >  1:4276813/1‑71 (MQ=255)
          tttATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAc                                                       >  1:4282446/1‑71 (MQ=255)
                cGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtctc                                                                 <  1:4379241/55‑1 (MQ=255)
                cGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                  <  1:4089609/70‑1 (MQ=255)
                   aTGAACGTTATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGgaga                                              <  1:3170685/71‑1 (MQ=255)
                   aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                          >  1:4931163/1‑59 (MQ=255)
                   aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGgaga                                              <  1:3056689/71‑1 (MQ=255)
                   aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGgaga                                              <  1:2469410/71‑1 (MQ=255)
                      aaCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                  >  1:1240561/1‑48 (MQ=255)
                      aaCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                  >  1:535208/1‑48 (MQ=255)
                      aaCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                  >  1:1842269/1‑48 (MQ=255)
                      aaCGTGATTATTCAGCCCTTGTTTTATGCCACATTCGGCAAAGCGTTTCTGCACATTAcg                                                      >  1:2918064/1‑60 (MQ=255)
                        cGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc                                         >  1:3205705/1‑71 (MQ=255)
                         gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                          <  1:2314979/53‑1 (MQ=255)
                         gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                          <  1:4561696/53‑1 (MQ=255)
                         gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                          <  1:4039090/53‑1 (MQ=255)
                         gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCt                                        <  1:603982/71‑1 (MQ=255)
                         gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCt                                        <  1:4061325/71‑1 (MQ=255)
                         gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGAAAAGCGTCTCTGCACATTACGCGCGGAGATACGCt                                        <  1:1459016/71‑1 (MQ=255)
                          tGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCtt                                       <  1:4118472/71‑1 (MQ=255)
                           gATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGcaca                                                           >  1:560437/1‑50 (MQ=255)
                            attattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAt                                                          <  1:1654517/50‑1 (MQ=255)
                            attattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAt                                                          <  1:3943326/50‑1 (MQ=255)
                             ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                        <  1:2307490/51‑1 (MQ=255)
                             ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                        <  1:4278263/51‑1 (MQ=255)
                             ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                        <  1:4993854/51‑1 (MQ=255)
                             ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                        <  1:3502550/51‑1 (MQ=255)
                             ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                        <  1:1242189/51‑1 (MQ=255)
                               attCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTctct                                                                >  1:3131971/1‑41 (MQ=255)
                               attCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTctct                                                                >  1:720799/1‑41 (MQ=255)
                                ttCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGAt                                             <  1:1203060/59‑1 (MQ=255)
                                ttCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGAt                                             <  1:4325503/59‑1 (MQ=255)
                                       cTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGga                                                <  1:525249/49‑1 (MQ=255)
                                       cTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGgaga                                              <  1:4270061/51‑1 (MQ=255)
                                       cTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                           >  1:3567515/1‑70 (MQ=255)
                                         tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgaga                                              >  1:3932372/1‑49 (MQ=255)
                                         tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGAGCGgaga                                              >  1:1774728/1‑49 (MQ=255)
                                          gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAg                        <  1:4397031/70‑1 (MQ=255)
                                          gTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGga                                                >  1:913376/1‑46 (MQ=255)
                                          gTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                             >  1:1749/1‑65 (MQ=255)
                                          gTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                             >  1:591907/1‑65 (MQ=255)
                                             ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaa                     <  1:3207965/70‑1 (MQ=255)
                                             ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaa                     <  1:1287199/70‑1 (MQ=255)
                                              tATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa                    <  1:273603/70‑1 (MQ=255)
                                              tATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa                    <  1:2603067/70‑1 (MQ=255)
                                              tATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa                    <  1:50710/70‑1 (MQ=255)
                                              tATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa                    <  1:141895/70‑1 (MQ=255)
                                               aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAg                  >  1:3692715/1‑71 (MQ=255)
                                                  ccccATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcg                <  1:527303/70‑1 (MQ=255)
                                                  ccccATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcg                <  1:1445720/70‑1 (MQ=255)
                                                   cccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGcgcg              <  1:3461454/71‑1 (MQ=255)
                                                     cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc            >  1:2096923/1‑71 (MQ=255)
                                                     cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc            >  1:2406119/1‑71 (MQ=255)
                                                       ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGAt                                             <  1:1860534/36‑1 (MQ=255)
                                                       ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGAt                                             <  1:237147/36‑1 (MQ=255)
                                                       ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGAt                                             <  1:2504801/36‑1 (MQ=255)
                                                       ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCg                            <  1:3267093/53‑1 (MQ=255)
                                                       ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcgc               <  1:4860807/66‑1 (MQ=255)
                                                       ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt          >  1:1557978/1‑71 (MQ=255)
                                                       ttCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATAc                                           <  1:2015541/38‑1 (MQ=255)
                                                         cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCg                            <  1:3069417/51‑1 (MQ=255)
                                                         cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt          >  1:2599512/1‑69 (MQ=255)
                                                         cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt          >  1:348638/1‑69 (MQ=255)
                                                         cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt        >  1:3593757/1‑71 (MQ=255)
                                                         cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt        >  1:1273151/1‑71 (MQ=255)
                                                         cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt        >  1:2493667/1‑71 (MQ=255)
                                                         cGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGc                                         >  1:2221238/1‑38 (MQ=255)
                                                         cGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCgtcgtc         <  1:5111296/70‑1 (MQ=255)
                                                          ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt        <  1:1346177/70‑1 (MQ=255)
                                                               aGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTa  <  1:2433105/71‑1 (MQ=255)
                                                                   tttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTc      <  1:3000389/63‑1 (MQ=255)
                                                                    ttCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc            <  1:4550742/56‑1 (MQ=255)
                                                                    |                                                                 
AGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTA  >  W3110S.gb/3639662‑3639795

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: