Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,657,715 A→G 8.8% G122G (GGT→GGC wzyE ← predicted Wzy protein involved in ECA polysaccharide chain elongation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,657,7150AG 8.8% 147.9 / 5.3 57G122G (GGT→GGCwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (2/3);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGACGGCGCGG  >  W3110S.gb/3657645‑3657780
                                                                      |                                                                 
tGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATa                                                                   >  1:2846856/1‑71 (MQ=255)
    ggAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCa                                                               <  1:824980/71‑1 (MQ=255)
    ggAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCa                                                               <  1:2148292/71‑1 (MQ=255)
      aGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATc                                                             <  1:272206/71‑1 (MQ=255)
      aGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATc                                                             <  1:1857327/71‑1 (MQ=255)
          aaGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAgg                                                          <  1:2608280/70‑1 (MQ=255)
            gCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCATCAGGAt                                                        >  1:1583749/1‑70 (MQ=255)
            gCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCATCAGGAt                                                        >  1:3258201/1‑70 (MQ=255)
                 aaCAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCAt                                                              <  1:1369515/59‑1 (MQ=255)
                 aaCAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCAt                                                              <  1:171220/59‑1 (MQ=255)
                 aaCAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCAt                                                              <  1:2404098/59‑1 (MQ=255)
                  aCAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCAc                                                 >  1:1032550/1‑71 (MQ=255)
                   cAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACg                                                >  1:822389/1‑71 (MQ=255)
                     gCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGt                                               >  1:601794/1‑70 (MQ=255)
                      cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:472810/71‑1 (MQ=255)
                      cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:588600/71‑1 (MQ=255)
                      cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:302761/71‑1 (MQ=255)
                       aaaaaGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:1928293/70‑1 (MQ=255)
                       aaaaaGCCGTTGTGCATGAAGAAGACGCCGACGCTTACCATCGCGATACCCATCAGGATCACCCa                                                  >  1:2044673/1‑65 (MQ=255)
                       aaaaaGCCGTTGTGCATGAAGAAGACGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCa                                                  >  1:460032/1‑65 (MQ=255)
                        aaaaGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAg                                            <  1:3849356/70‑1 (MQ=255)
                         aaaGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGGGCGATACCTATCAGGATCACCCACGTAAGAt                                          <  1:192721/71‑1 (MQ=255)
                         aaaGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGAt                                          <  1:820743/71‑1 (MQ=255)
                         aaaGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGAt                                          <  1:821356/71‑1 (MQ=255)
                         aaaGCCGTTGCGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGAt                                          <  1:3817639/71‑1 (MQ=255)
                             ccGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:830130/64‑1 (MQ=255)
                               gTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTg                                        >  1:1405113/1‑67 (MQ=255)
                               gTTGTGCATGAAGAAGATACCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTg                                        >  1:517022/1‑67 (MQ=255)
                                ttGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGa                                           <  1:3663230/63‑1 (MQ=255)
                                ttGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGa                                           <  1:3657795/63‑1 (MQ=255)
                                ttGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGtctc                                   <  1:3661338/71‑1 (MQ=255)
                                  gtgCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTgg                                       >  1:2918484/1‑65 (MQ=255)
                                  gtgCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTgg                                       >  1:3646389/1‑65 (MQ=255)
                                  gtgCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCa                                 >  1:480651/1‑71 (MQ=255)
                                     cATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATAACCCACGTAAGATTGGt                                      >  1:3495130/1‑63 (MQ=255)
                                      aTGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCAcc                                                    <  1:1813971/48‑1 (MQ=255)
                                        gaagaaGATGCCGACGCTTACCAGCGCGATACCCATCAg                                                           <  1:1028186/39‑1 (MQ=255)
                                        gaagaaGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGt                                      >  1:2266159/1‑60 (MQ=255)
                                             agaTGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATg                      >  1:767575/1‑71 (MQ=255)
                                              gaTGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATgg                     >  1:2916349/1‑71 (MQ=255)
                                               aTGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGt                                               <  1:2984825/44‑1 (MQ=255)
                                                tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                    >  1:3699530/1‑70 (MQ=255)
                                                tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                    >  1:3851919/1‑70 (MQ=255)
                                                tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                    >  1:1895402/1‑70 (MQ=255)
                                                tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                    >  1:2588458/1‑70 (MQ=255)
                                                tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                    >  1:1360564/1‑70 (MQ=255)
                                                 gCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTa                   >  1:2039836/1‑70 (MQ=255)
                                                    gACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCAt                       >  1:3070417/1‑63 (MQ=255)
                                                      cGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACa               >  1:3343828/1‑69 (MQ=255)
                                                      cGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACa               >  1:118143/1‑69 (MQ=255)
                                                       gCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:3371757/38‑1 (MQ=255)
                                                       gCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:1833349/38‑1 (MQ=255)
                                                       gCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                             <  1:181900/38‑1 (MQ=255)
                                                       gCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCg            >  1:3785811/1‑71 (MQ=255)
                                                           aCCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGAc         <  1:2163520/70‑1 (MQ=255)
                                                               gcgcGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGACGgcgc    <  1:2829513/71‑1 (MQ=255)
                                                                 gcgATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGACGGCGCgg  >  1:3469320/1‑71 (MQ=255)
                                                                     tACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCAt                       <  1:2270461/46‑1 (MQ=255)
                                                                      |                                                                 
TGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGACGGCGCGG  >  W3110S.gb/3657645‑3657780

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: