Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,657,709 C→A 8.9% A124A (GCG→GCT 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,7090CA 8.9% 169.9 / 2.4 57A124A (GCG→GCTwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base A (2/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

ACTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGACGGCGC  >  W3110S.gb/3657643‑3657778
                                                                  |                                                                     
aCTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGcgcg                                                                    >  1:3051102/1‑70 (MQ=255)
  tGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATa                                                                 >  1:2846856/1‑71 (MQ=255)
      ggAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCa                                                             <  1:2148292/71‑1 (MQ=255)
      ggAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCa                                                             <  1:824980/71‑1 (MQ=255)
        aGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATc                                                           <  1:272206/71‑1 (MQ=255)
        aGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATc                                                           <  1:1857327/71‑1 (MQ=255)
          ttAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAgcgc                                                                     <  1:1025127/59‑1 (MQ=255)
            aaGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAgg                                                        <  1:2608280/70‑1 (MQ=255)
              gCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAgcgc                                                                     <  1:3613248/55‑1 (MQ=255)
              gCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCATCAGGAt                                                      >  1:3258201/1‑70 (MQ=255)
              gCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCATCAGGAt                                                      >  1:1583749/1‑70 (MQ=255)
                   aaCAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCAt                                                            <  1:171220/59‑1 (MQ=255)
                   aaCAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTCACCAGAGCGATGCCCAt                                                            <  1:1369515/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)
                      aGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAgcgc                                                                     <  1:3712474/47‑1 (MQ=255)
                       gCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGt                                             >  1:601794/1‑70 (MQ=255)
                        cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGa                                                                   >  1:73910/1‑47 (MQ=255)
                        cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGa                                                                   >  1:1539470/1‑47 (MQ=255)
                        cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                           <  1:302761/71‑1 (MQ=255)
                        cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                           <  1:472810/71‑1 (MQ=255)
                        cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                           <  1:588600/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:3657795/63‑1 (MQ=255)
                                  ttGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGa                                         <  1:3663230/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:1360564/1‑70 (MQ=255)
                                                  tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                  >  1:2588458/1‑70 (MQ=255)
                                                  tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                  >  1:3851919/1‑70 (MQ=255)
                                                  tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                  >  1:3699530/1‑70 (MQ=255)
                                                  tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt                  >  1:1895402/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:1833349/38‑1 (MQ=255)
                                                         gCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                           <  1:181900/38‑1 (MQ=255)
                                                         gCTTACCAGCGCGATACCCATCAGGATCACCCACGTaa                                           <  1:3371757/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)
                                                                  |                                                                     
ACTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCGGACGGCGC  >  W3110S.gb/3657643‑3657778

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: