Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,657,703 T→C 7.8% V126V (GTA→GTG 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,7030TC 7.8% 170.0 / 4.9 64V126V (GTA→GTGwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base C (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AAAGATCTGACTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCG  >  W3110S.gb/3657634‑3657770
                                                                     |                                                                   
aaaGATCTGACTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTa                                                                    >  1:243079/1‑71 (MQ=255)
aaaGATCTGACTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTa                                                                    >  1:1781209/1‑71 (MQ=255)
         aCTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGcgcg                                                            >  1:3051102/1‑70 (MQ=255)
           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)
                  gTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTa                                                                    >  1:2019713/1‑53 (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)
                          ggAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCtt                                                                     <  1:322674/44‑1 (MQ=255)
                          ggAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCtt                                                                     <  1:2807104/44‑1 (MQ=255)
                           gAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTa                                                                    <  1:3826922/44‑1 (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)
                               aGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAgcgc                                                             <  1:3712474/47‑1 (MQ=255)
                                gCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGt                                     >  1:601794/1‑70 (MQ=255)
                                 cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGa                                                           >  1:1539470/1‑47 (MQ=255)
                                 cAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGa                                                           >  1:73910/1‑47 (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:821356/71‑1 (MQ=255)
                                    aaaGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGAt                                <  1:820743/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:3646389/1‑65 (MQ=255)
                                             gtgCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTgg                             >  1:2918484/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:1895402/1‑70 (MQ=255)
                                                           tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt          >  1:1360564/1‑70 (MQ=255)
                                                           tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt          >  1:3851919/1‑70 (MQ=255)
                                                           tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt          >  1:2588458/1‑70 (MQ=255)
                                                           tGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGt          >  1:3699530/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)
                                                                     |                                                                   
AAAGATCTGACTGTAGGAGTTAAGCCGGAACAGCAAAAAGCCGTTGTGCATGAAGAAGATGCCGACGCTTACCAGCGCGATACCCATCAGGATCACCCACGTAAGATTGGTCTCCACGCGGTTCATGGTAAACAGCG  >  W3110S.gb/3657634‑3657770

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: