Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,488,502 C→T 41.3% L47L (CTG→TTG)  yijO → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,488,5020CT41.3% 24.7 / 53.3 46L47L (CTG→TTG) yijOpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (6/21);  new base T (9/10);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCAGCAGTAACGGCCCGGTGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCGATGTGCTGTACGTTCC  >  W3110S.gb/3488437‑3488569
                                                                 |                                                                   
ccAGCAGTAACGGCCCGTTGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGa                                                                 <  1:1595587/70‑1 (MQ=255)
ccAGCAGTAACGGCCCGGTGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGa                                                                 <  1:1692941/70‑1 (MQ=255)
         aCGGCCCGGTGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTCAGttt                                                       >  1:1095830/1‑71 (MQ=255)
            gCCCGGTGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTttgtt                                                    <  1:1441547/71‑1 (MQ=255)
                  tGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATAcc                                              <  1:1980146/71‑1 (MQ=255)
                  tGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATAcc                                              <  1:2609139/71‑1 (MQ=255)
                  tGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATAcc                                              <  1:874319/71‑1 (MQ=255)
                    cccGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCgg                                            <  1:2216586/71‑1 (MQ=255)
                        aTCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGa                                                  <  1:214912/61‑1 (MQ=255)
                        aTCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGa                                                  <  1:3503994/61‑1 (MQ=255)
                        aTCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGa                                                  <  1:1582518/61‑1 (MQ=255)
                          cTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTg                                       >  1:2539096/1‑70 (MQ=255)
                          cTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGCAAGGTGAGTTTGTTGATACCGGCGCTg                                       >  1:1470109/1‑70 (MQ=255)
                          cTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTg                                       >  1:863216/1‑70 (MQ=255)
                             gCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTttgtt                                                    <  1:1533269/54‑1 (MQ=255)
                             gCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTttgtt                                                    <  1:914299/54‑1 (MQ=255)
                             gCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTttgtt                                                    <  1:2901331/54‑1 (MQ=255)
                              cGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCa                                  <  1:2951246/71‑1 (MQ=255)
                              cGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCa                                  <  1:3762323/71‑1 (MQ=255)
                               gTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGa                                                  >  1:1392476/1‑54 (MQ=255)
                               gTATCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGa                                                  >  1:2910255/1‑54 (MQ=255)
                                tATCAGGTCGATTTTCCACGTCTGGAAATTGTGCTGGa                                                                 <  1:2324352/38‑1 (MQ=255)
                                  tCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGaa                                                                <  1:1732881/37‑1 (MQ=255)
                                  tCAGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGt                              >  1:1435568/1‑71 (MQ=255)
                                  tCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGttt                                                       <  1:3717328/46‑1 (MQ=255)
                                  tCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGttt                                                       <  1:1302706/46‑1 (MQ=255)
                                    aGGTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGt                              >  1:2318188/1‑69 (MQ=255)
                                     ggTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAg                           <  1:1945673/71‑1 (MQ=255)
                                      gTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGt                              <  1:217149/67‑1 (MQ=255)
                                      gTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGt                              <  1:3377181/67‑1 (MQ=255)
                                      gTCGATTTTCCACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGt                              <  1:3381655/67‑1 (MQ=255)
                                      gTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAg                                                          <  1:2084015/39‑1 (MQ=255)
                                      gTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAg                                                          <  1:3581862/39‑1 (MQ=255)
                                        cGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGa                                                  >  1:1236383/1‑45 (MQ=255)
                                        cGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGa                                                  >  1:3400740/1‑45 (MQ=255)
                                        cGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGa                                                  >  1:2384398/1‑45 (MQ=255)
                                        cGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGa                                                  >  1:14291/1‑45 (MQ=255)
                                        cGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTc                        >  1:94165/1‑71 (MQ=255)
                                              tCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGc                   <  1:2046589/70‑1 (MQ=255)
                                               ccACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCg                  <  1:3308628/70‑1 (MQ=255)
                                               ccACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCg                  <  1:2567763/70‑1 (MQ=255)
                                               ccACGGCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCg                  <  1:2140709/70‑1 (MQ=255)
                                                 aCGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGa                                                  <  1:3491534/36‑1 (MQ=255)
                                                    gCTGGAAATTGTGTTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCGAtgt              >  1:462962/1‑69 (MQ=255)
                                                    gCTGGAAATTGTGTCGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCGAtgtg             >  1:1700679/1‑70 (MQ=255)
                                                               tgCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCGATGTGCTGTACGTTcc  <  1:1809527/70‑1 (MQ=255)
                                                                 |                                                                   
CCAGCAGTAACGGCCCGGTGCCCGATCTTGCGTATCAGGTCGATTTTCCACGGCTGGAAATTGTGCTGGAAGGTGAGTTTGTTGATACCGGCGCTGGAGCAACGTTAGTTCCCGGCGATGTGCTGTACGTTCC  >  W3110S.gb/3488437‑3488569

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: