Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,639,730 T→C 41.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,7300TC41.7% 7.3 / 25.9 24K222R (AAA→AGA) xerCsite‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base C (2/8);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

ATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGC  >  W3110S.gb/3639671‑3639793
                                                           |                                                               
aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa                                                      >  1:1071168/1‑71 (MQ=255)
   tATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcg                                                    >  1:1624103/1‑70 (MQ=255)
          aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                        <  1:834414/59‑1 (MQ=255)
                           gCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc                           <  1:474966/71‑1 (MQ=255)
                             ccTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                         <  1:1078034/71‑1 (MQ=255)
                              cTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                         >  1:674155/1‑70 (MQ=255)
                               ttGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGga                                              <  1:1364490/48‑1 (MQ=255)
                                tGTTTTATGCCTCATTCGGCAAAGCGTTTCTGcaca                                                         >  1:1548061/1‑36 (MQ=255)
                                tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcaca                                                         >  1:1586761/1‑36 (MQ=255)
                                tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcaca                                                         >  1:246481/1‑36 (MQ=255)
                                 gTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACAt                                                        >  1:780555/1‑36 (MQ=255)
                                  ttttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGa                     <  1:1497732/70‑1 (MQ=255)
                                  ttttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGa                     <  1:1112088/70‑1 (MQ=255)
                                    ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaa                   >  1:1006995/1‑70 (MQ=255)
                                     tATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa                  <  1:692224/70‑1 (MQ=255)
                                      aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa                         >  1:357326/1‑62 (MQ=255)
                                          cccATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATTCGCTTGCCCAGTTTCGACAGAAATAGcgcg            <  1:954908/71‑1 (MQ=255)
                                          cccATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGcgcg            <  1:156572/71‑1 (MQ=255)
                                            cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc          >  1:534266/1‑71 (MQ=255)
                                                 ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCtt     <  1:8705/71‑1 (MQ=255)
                                                 ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCtt     <  1:1017439/71‑1 (MQ=255)
                                                 ggCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCtt     <  1:313740/71‑1 (MQ=255)
                                                    aaaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGc  <  1:635468/71‑1 (MQ=255)
                                                      aGCGTTTCTGCACATTACGCGCGGAGAAACGCTTGCCCAGTTTCGACAGAAAAAGCGcgt           >  1:1215344/1‑60 (MQ=255)
                                                           |                                                               
ATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGC  >  W3110S.gb/3639671‑3639793

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: