Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,657,262 C→T 25.9% E273E (GAG→GAA 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,2620CT25.9% 156.7 / 56.1 86E273E (GAG→GAAwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base C (35/28);  new base T (12/10);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCACATT  >  W3110S.gb/3657195‑3657330
                                                                   |                                                                    
cgcgGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGAtcctcc                                                                   >  1:4237296/1‑71 (MQ=255)
 gcgGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTccc                                                                  <  1:4963938/71‑1 (MQ=255)
      cAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGc                                                             <  1:326196/71‑1 (MQ=255)
       aaTTGGAGCCAGGCCCTGTAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGc                                                             >  1:819244/1‑70 (MQ=255)
       aaTTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGc                                                             >  1:559137/1‑70 (MQ=255)
       aaTTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGc                                                             >  1:700431/1‑70 (MQ=255)
        aTTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCga                                                           <  1:1448995/71‑1 (MQ=255)
         ttGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCga                                                           <  1:3230834/70‑1 (MQ=255)
           ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCgag                                                          >  1:3630011/1‑69 (MQ=255)
            gAGCCAGGCCCTGGAAGTCGTTGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGaa                                                        <  1:4947025/70‑1 (MQ=255)
            gAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGaa                                                        <  1:3436229/70‑1 (MQ=255)
            gAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGaa                                                        <  1:4358206/70‑1 (MQ=255)
              gCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGgt                                                     <  1:2789630/71‑1 (MQ=255)
               ccAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGgtg                                                    <  1:2410135/71‑1 (MQ=255)
                   gCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCgag                                                          <  1:3483914/61‑1 (MQ=255)
                    cccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATtct                                                                     >  1:3275543/1‑49 (MQ=255)
                    cccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATtct                                                                     >  1:4930722/1‑49 (MQ=255)
                    cccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgcg                                               <  1:5134665/71‑1 (MQ=255)
                     ccTTGAAGTCGATGTTGTCGTAGTTCTGCAACATCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                                >  1:2229556/1‑69 (MQ=255)
                     ccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                                >  1:4116187/1‑69 (MQ=255)
                     ccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                                >  1:4820346/1‑69 (MQ=255)
                     ccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACAGCGAGAAGGTGTcgc                                                >  1:1196049/1‑69 (MQ=255)
                         gAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGt                                            <  1:1570024/69‑1 (MQ=255)
                          aaGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCa                                                                 >  1:505399/1‑47 (MQ=255)
                            gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACgg                                                              >  1:205035/1‑48 (MQ=255)
                            gTCGATGTTGTCGTAGCTCTGCAACAGCAACGCCAGATTCTCCCACgg                                                              >  1:2459675/1‑48 (MQ=255)
                              cGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTgaga                                        >  1:4840137/1‑68 (MQ=255)
                               gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                      >  1:3969625/1‑69 (MQ=255)
                               gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                      >  1:520406/1‑69 (MQ=255)
                                 tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATagag                                   >  1:2713406/1‑70 (MQ=255)
                                  gttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga                                  >  1:2431832/1‑70 (MQ=255)
                                  gttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga                                  >  1:4656409/1‑70 (MQ=255)
                                   ttgtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGaa                                                        >  1:751446/1‑47 (MQ=255)
                                   ttgtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAg                                                       <  1:3925429/48‑1 (MQ=255)
                                   ttgtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGAt                                       >  1:349385/1‑64 (MQ=255)
                                   ttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaa                                 <  1:4060153/70‑1 (MQ=255)
                                   ttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaa                                 <  1:857133/70‑1 (MQ=255)
                                     gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg                              >  1:5842/1‑71 (MQ=255)
                                     gtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAg                                                       <  1:407108/46‑1 (MQ=255)
                                     gtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAg                                                       <  1:1533089/46‑1 (MQ=255)
                                     gtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg                              >  1:1642710/1‑71 (MQ=255)
                                     gtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg                              >  1:4491743/1‑71 (MQ=255)
                                      tCGTAGCTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg                              <  1:971212/70‑1 (MQ=255)
                                        gTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAg                                             >  1:1455125/1‑53 (MQ=255)
                                        gTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAg                                             >  1:1942499/1‑53 (MQ=255)
                                         tAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAg                                             <  1:2718357/52‑1 (MQ=255)
                                         tAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAg                                             <  1:1585836/52‑1 (MQ=255)
                                           gTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAg                                                       >  1:215276/1‑40 (MQ=255)
                                            ttCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                                <  1:2399736/46‑1 (MQ=255)
                                              cTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     >  1:594085/1‑71 (MQ=255)
                                              cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaa                                 >  1:3589396/1‑59 (MQ=255)
                                              cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     >  1:1133549/1‑71 (MQ=255)
                                              cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     >  1:594067/1‑71 (MQ=255)
                                              cTGCAACAGCAACGCCAGATTCTCCCACGGAGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     >  1:3727403/1‑71 (MQ=255)
                                                 cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     <  1:4854949/68‑1 (MQ=255)
                                                 cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     <  1:2511174/68‑1 (MQ=255)
                                                 cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACATATAGAACGCt                   <  1:4777252/70‑1 (MQ=255)
                                                  aaCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                      >  1:3631769/1‑50 (MQ=255)
                                                  aaCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                      >  1:2471381/1‑50 (MQ=255)
                                                  aaCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                      >  1:453331/1‑50 (MQ=255)
                                                  aaCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCtt                  <  1:1294256/70‑1 (MQ=255)
                                                  aaCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTc                 <  1:4832375/71‑1 (MQ=255)
                                                   aCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                     <  1:2535641/66‑1 (MQ=255)
                                                   aCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCa                <  1:2080751/71‑1 (MQ=255)
                                                    cAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCAt               >  1:1524116/1‑71 (MQ=255)
                                                    cAGAAACGCCAGATTATCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTc                 >  1:1970497/1‑69 (MQ=255)
                                                     aGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                              >  1:3385674/1‑39 (MQ=255)
                                                     aGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATc              >  1:2606186/1‑71 (MQ=255)
                                                      gCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTgaga                                        >  1:3541634/1‑44 (MQ=255)
                                                      gCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTgaga                                        >  1:2038536/1‑44 (MQ=255)
                                                      gCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcg             >  1:3381764/1‑71 (MQ=255)
                                                      gCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATc              >  1:2168600/1‑70 (MQ=255)
                                                         aCGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgcag          >  1:4806466/1‑69 (MQ=255)
                                                           gCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgccctc        <  1:4529637/71‑4 (MQ=255)
                                                           gCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCt        >  1:3890424/1‑71 (MQ=255)
                                                            ccAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTc       <  1:4431622/71‑1 (MQ=255)
                                                            ccAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTc       <  1:2797209/71‑1 (MQ=255)
                                                            ccAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTgaga                                        >  1:3286395/1‑38 (MQ=255)
                                                            ccAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg                              <  1:2232831/48‑1 (MQ=255)
                                                             cAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgc            >  1:1580324/1‑65 (MQ=255)
                                                             cAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTca      >  1:1760305/1‑71 (MQ=255)
                                                             cAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTca      <  1:1484149/71‑1 (MQ=255)
                                                             cAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTc       <  1:2256105/70‑1 (MQ=255)
                                                              aGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTca      >  1:256936/1‑70 (MQ=255)
                                                               gATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTcaca    <  1:3865546/71‑1 (MQ=255)
                                                               gATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTcaca    <  1:2396312/71‑1 (MQ=255)
                                                                aTTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTTTAGAACGCTTCATCGCCGCTcaca    <  1:1929552/70‑1 (MQ=255)
                                                                aTTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAg                         <  1:2754788/49‑1 (MQ=255)
                                                                 ttCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCAcatt  >  1:4259273/1‑71 (MQ=255)
                                                                  tctcCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgc            <  1:3446915/60‑1 (MQ=255)
                                                                  tctcCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgc            <  1:3339340/60‑1 (MQ=255)
                                                                   |                                                                    
CGCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCACATT  >  W3110S.gb/3657195‑3657330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: