Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,639,730 T→C 10.6% 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,7300TC10.6% 171.1 / 11.7 66K222R (AAA→AGA) xerCsite‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (4/3);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCA  >  W3110S.gb/3639660‑3639798
                                                                      |                                                                    
cGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtt                                                                      <  1:1601252/71‑1 (MQ=255)
cGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtt                                                                      >  1:2048812/1‑71 (MQ=255)
 gAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                     <  1:988673/71‑1 (MQ=255)
 gAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt                                                                     <  1:2071013/71‑1 (MQ=255)
  aGTGACGTAATTTATGAGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTc                                                                    <  1:279887/71‑1 (MQ=255)
   gTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCt                                                                   >  1:2844512/1‑71 (MQ=255)
    tGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg                                                                  >  1:1594197/1‑71 (MQ=255)
      aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca                                                                >  1:1987175/1‑71 (MQ=255)
      aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGca                                                                >  1:4142793/1‑71 (MQ=255)
        gTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcaca                                                              >  1:19870/1‑71 (MQ=255)
                  cGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                     <  1:767229/70‑1 (MQ=255)
                  cGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg                                                     <  1:2194307/70‑1 (MQ=255)
                  cGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                                     <  1:2088540/70‑1 (MQ=255)
                   ggATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGga                                                   <  1:3304459/71‑1 (MQ=255)
                    gATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgag                                                  <  1:3238763/71‑1 (MQ=255)
                    gATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgag                                                  <  1:1127600/71‑1 (MQ=255)
                      tGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgaga                                                 >  1:3822456/1‑70 (MQ=255)
                        aaCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATAc                                              >  1:1401723/1‑71 (MQ=255)
                         aCGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAtt                                                            >  1:804923/1‑56 (MQ=255)
                          cGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc                                            <  1:3773521/71‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgc                                                      <  1:2099829/60‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgc                                                      <  1:4237074/60‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                           <  1:4100942/71‑1 (MQ=255)
                           gTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCt                                           <  1:3828569/71‑1 (MQ=255)
                             gATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTg                                         >  1:3027935/1‑71 (MQ=255)
                              attattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtt                                                                      <  1:2954809/41‑1 (MQ=255)
                              attattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGtt                                                                      <  1:2819549/41‑1 (MQ=255)
                              attattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGc                                        <  1:3276968/71‑1 (MQ=255)
                               ttattaAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGcc                                       >  1:1836417/1‑71 (MQ=255)
                               ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTc                                                                    <  1:1873199/42‑1 (MQ=255)
                                 attCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGga                                                   >  1:297447/1‑57 (MQ=255)
                                 attCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGc                                            <  1:349157/64‑1 (MQ=255)
                                   tCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgc                                                      <  1:4106970/52‑1 (MQ=255)
                                   tCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc                                            <  1:1120880/62‑1 (MQ=255)
                                   tCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCt                                           <  1:3268873/63‑1 (MQ=255)
                                    cAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                           >  1:4016263/1‑46 (MQ=255)
                                    cAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGtt                                  >  1:3728810/1‑71 (MQ=255)
                                     cagccTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGttt                                 >  1:1758278/4‑71 (MQ=255)
                                     aGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgaga                                                 <  1:916748/55‑1 (MQ=255)
                                     aGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATAc                                              >  1:190319/1‑58 (MQ=255)
                                      gCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                                <  1:12241/71‑1 (MQ=255)
                                       cccTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCg                               <  1:3583909/71‑1 (MQ=255)
                                        ccTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgaga                                                 <  1:4092035/52‑1 (MQ=255)
                                         cTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACg                                             >  1:3863714/1‑55 (MQ=255)
                                         cTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                              >  1:1020706/1‑70 (MQ=255)
                                          ttGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAc                                                          <  1:2975579/41‑1 (MQ=255)
                                           tGTTTTATGCCCCATTCGGCAAAGAGTTTCTGCACATTAcgcgcg                                                     >  1:934447/1‑45 (MQ=255)
                                            gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcg                                                         >  1:3348034/1‑40 (MQ=255)
                                            gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATa                                               >  1:1693968/1‑50 (MQ=255)
                                            gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACa                            <  1:3311935/69‑1 (MQ=255)
                                              tttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaa                         >  1:881009/1‑70 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCa                                     >  1:2093842/1‑57 (MQ=255)
                                               ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                              >  1:2551833/1‑64 (MQ=255)
                                                 aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgaga                                                 <  1:2909961/43‑1 (MQ=255)
                                                 aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaaa                      <  1:1102562/70‑1 (MQ=255)
                                                      ccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcgc                  >  1:1317240/1‑69 (MQ=255)
                                                       cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc               <  1:4099446/71‑1 (MQ=255)
                                                         ttCGGCAAAGCGTTTCTGCACATTACGCGCGgaga                                                 <  1:3596202/35‑1 (MQ=255)
                                                         ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt             >  1:3676668/1‑71 (MQ=255)
                                                           cGGCAAAGCGTTTCTGTACATTACGCGCGGAGATACGCATGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt           >  1:1392189/1‑71 (MQ=255)
                                                           cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCgtcgtc            <  1:1710886/70‑1 (MQ=255)
                                                           cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt           >  1:3905818/1‑71 (MQ=255)
                                                           cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt           >  1:3623330/1‑71 (MQ=255)
                                                           cGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt           >  1:863541/1‑71 (MQ=255)
                                                           cGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcg                   >  1:961794/1‑63 (MQ=255)
                                                            ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGGCt           <  1:2608646/70‑1 (MQ=255)
                                                            ggCAAAGCGTTTCAGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc               <  1:3886799/66‑1 (MQ=255)
                                                               aaaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGc       >  1:460071/1‑71 (MQ=255)
                                                               aaaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGc       <  1:4280940/71‑1 (MQ=255)
                                                               aaaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGc       >  1:2885359/1‑71 (MQ=255)
                                                                aaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCt      <  1:3169891/71‑1 (MQ=255)
                                                                    gTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCa  <  1:3668421/71‑1 (MQ=255)
                                                                      |                                                                    
CGAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTACCA  >  W3110S.gb/3639660‑3639798

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: