Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,419,194 C→T 9.1% 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
*minE2,419,1940CT 9.1% 191.8 / 7.8 66E273E (GAG→GAAwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base T (4/2);  total (39/27)
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

AATCGCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCACAT  >  minE/2419124‑2419261
                                                                      |                                                                   
aaTCGCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATtc                                                                     <  1:3695645/71‑1 (MQ=255)
aaTCGCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATtc                                                                     <  1:2693787/71‑1 (MQ=255)
      ggACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCAc                                                               >  1:3839652/1‑71 (MQ=255)
       gACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACg                                                              >  1:547457/1‑71 (MQ=255)
        aCAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACg                                                              >  1:472112/1‑70 (MQ=255)
        aCAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACg                                                              >  1:2009103/1‑70 (MQ=255)
         cAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGc                                                            >  1:2556921/1‑71 (MQ=255)
          aaTTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATtct                                                                    <  1:334234/62‑1 (MQ=255)
              ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCgag                                                         >  1:3156675/1‑69 (MQ=255)
                 gCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGgt                                                    <  1:3172485/71‑1 (MQ=255)
                    aGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTc                                                 >  1:898086/1‑71 (MQ=255)
                     ggCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTccc                                                                 >  1:2553226/1‑54 (MQ=255)
                      gCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                               <  1:1580691/71‑1 (MQ=255)
                      gCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                               >  1:2009199/1‑71 (MQ=255)
                      gCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                               <  1:3449674/71‑1 (MQ=255)
                         cTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAg                                                      >  1:352242/1‑61 (MQ=255)
                         cTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAg                                            >  1:1163538/1‑71 (MQ=255)
                              aGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc                                               <  1:1986602/63‑1 (MQ=255)
                               gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCgag                                                         <  1:2595752/52‑1 (MQ=255)
                                tCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                     <  1:2858406/71‑1 (MQ=255)
                                tCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTgaga                                       >  1:1460024/1‑69 (MQ=255)
                                  gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCgaga                                                        >  1:3383449/1‑50 (MQ=255)
                                  gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCgaga                                                        <  1:3109236/50‑1 (MQ=255)
                                  gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                     >  1:1312776/1‑69 (MQ=255)
                                    tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATagat                                  >  1:1891640/1‑69 (MQ=255)
                                    tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATagag                                  >  1:2183249/1‑70 (MQ=255)
                                    tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATagag                                  >  1:3866120/1‑70 (MQ=255)
                                     gttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa                                     <  1:2916502/66‑1 (MQ=255)
                                      ttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaaa                               <  1:886322/71‑1 (MQ=255)
                                      ttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaa                                <  1:3838462/70‑1 (MQ=255)
                                      ttgtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaa                                <  1:2588538/70‑1 (MQ=255)
                                        gtCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACgg                                                             >  1:1376917/1‑39 (MQ=255)
                                         tCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg                             >  1:1052483/1‑70 (MQ=255)
                                           gTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTAGCGAGTGAGATaga                                   >  1:1488441/1‑62 (MQ=255)
                                             aGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATaga                                   >  1:1164571/1‑60 (MQ=255)
                                              gTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga                                 <  1:3084963/61‑1 (MQ=255)
                                              gTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAg                        <  1:3184896/70‑1 (MQ=255)
                                              gTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAg                        >  1:2254335/1‑70 (MQ=255)
                                               ttCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGgtgt                                                  <  1:2232859/43‑1 (MQ=255)
                                                tCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaaa                               >  1:474380/1‑61 (MQ=255)
                                                 cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGgtgt                                                  <  1:2543902/41‑1 (MQ=255)
                                                 cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAAc                     <  1:1118108/70‑1 (MQ=255)
                                                 cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                    >  1:2130602/1‑71 (MQ=255)
                                                 cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                    >  1:492235/1‑71 (MQ=255)
                                                 cTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                    >  1:94378/1‑71 (MQ=255)
                                                    cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGta                           <  1:3836153/61‑1 (MQ=255)
                                                    cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGta                           >  1:1050380/1‑61 (MQ=255)
                                                    cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACg                    <  1:2442703/68‑1 (MQ=255)
                                                    cAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCtt                 >  1:1038969/1‑71 (MQ=255)
                                                      aCAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGtat                          >  1:293822/1‑60 (MQ=255)
                                                       cAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCAt              <  1:1775859/71‑1 (MQ=255)
                                                       cAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCCCGAGTGAGATAGAGAAACGTATAGAACGCTTCAt              <  1:2011107/71‑1 (MQ=255)
                                                        aGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATc             <  1:2494682/71‑1 (MQ=255)
                                                          cAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgc           >  1:963340/1‑71 (MQ=255)
                                                           aaCGCCAGATTCTCCCACGGCGATAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgcc          <  1:2072814/71‑1 (MQ=255)
                                                            aCGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgccg         >  1:2603502/1‑71 (MQ=255)
                                                            aCGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATcgccg         <  1:591803/71‑1 (MQ=255)
                                                              gCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCa               >  1:2853132/1‑63 (MQ=255)
                                                               ccAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTg                                          >  1:975478/1‑35 (MQ=255)
                                                               ccAGATTTTCCCACGGCGAGAAGGTGTCGAGAGTg                                          >  1:2275016/1‑35 (MQ=37)
                                                               ccAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCa               <  1:3695965/62‑1 (MQ=255)
                                                                cAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTca     <  1:3377227/71‑1 (MQ=255)
                                                                 aGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTcac    <  1:866637/71‑1 (MQ=255)
                                                                  gATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTcaca   <  1:3951470/71‑1 (MQ=255)
                                                                  gATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTc                >  1:2137652/1‑58 (MQ=255)
                                                                  gATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTcaca   >  1:11021/1‑71 (MQ=255)
                                                                   aTTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCAcat  >  1:1772434/1‑71 (MQ=255)
                                                                    ttCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga                                 >  1:3920318/1‑39 (MQ=255)
                                                                    ttCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga                                 >  1:3372326/1‑39 (MQ=255)
                                                                      |                                                                   
AATCGCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCACAT  >  minE/2419124‑2419261

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