Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,010,113 C→T 10.7% A100V (GCC→GTC)  yhiI → predicted HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,010,1130CT10.7% 166.2 / 7.3 56A100V (GCC→GTC) yhiIpredicted HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (3/3);  total (29/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

AAGATGGATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTGG  >  W3110S.gb/4010043‑4010178
                                                                      |                                                                 
aaGATGGATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAgcgc                                                                   <  1:3014224/71‑1 (MQ=255)
aaGATGGATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAgcgc                                                                   <  1:2399943/71‑1 (MQ=255)
  gATGGATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCTCCg                                                                 >  1:3530093/1‑71 (MQ=255)
      gATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAgcgc                                                                   <  1:3625366/65‑1 (MQ=255)
      gATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTtgc                                                             <  1:2993423/71‑1 (MQ=255)
      gATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTtgc                                                             <  1:3149676/71‑1 (MQ=255)
      gATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTtgc                                                             <  1:1542088/71‑1 (MQ=255)
           tCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGcc                                                                  <  1:4652704/61‑1 (MQ=255)
           tCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGcc                                                                  <  1:3735454/61‑1 (MQ=255)
              cgTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCgcgc                                                     <  1:2539177/71‑1 (MQ=255)
               gtgtTGCCGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCgcgc                                                     <  1:3795146/70‑1 (MQ=255)
                  ttGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGc                                                 >  1:1683902/1‑71 (MQ=255)
                    gCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCtt                                               >  1:2059731/1‑71 (MQ=255)
                    gCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCtt                                               >  1:2254388/1‑71 (MQ=255)
                    gCAGGAACAGCGACTGGAAGCCATAGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCa                                                    >  1:2228373/1‑66 (MQ=255)
                     cAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCttt                                              <  1:2588139/71‑1 (MQ=255)
                      aGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTg                                             <  1:2447788/71‑1 (MQ=255)
                      aGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTg                                             <  1:4558587/71‑1 (MQ=255)
                       ggAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGc                                            <  1:4252275/71‑1 (MQ=255)
                       ggAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGc                                            <  1:3107915/71‑1 (MQ=255)
                        gAACAGCGACTGGAATCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGc                                            <  1:4737075/70‑1 (MQ=255)
                        gAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCt                                           >  1:1815928/1‑71 (MQ=255)
                          aCAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCa                                                    >  1:1505972/1‑60 (MQ=255)
                          aCAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTg                                          >  1:1239272/1‑70 (MQ=255)
                            aGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAg                                       >  1:1825829/1‑71 (MQ=255)
                             gCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAg                                       >  1:2864963/1‑70 (MQ=255)
                              cGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCttt                                              <  1:2165142/62‑1 (MQ=255)
                                 cTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAAc                                   >  1:1101423/1‑70 (MQ=255)
                                  tGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCtt                                               <  1:299598/57‑1 (MQ=255)
                                    gAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGAc                                <  1:3981901/70‑1 (MQ=255)
                                    gAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGAc                                <  1:1110316/70‑1 (MQ=255)
                                      aGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACaaa                             <  1:2066160/71‑1 (MQ=255)
                                         cATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCg                          <  1:609781/71‑1 (MQ=255)
                                         cATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCg                          <  1:3558816/71‑1 (MQ=255)
                                          aTCGCGCAAATCAAAGAGGCACAAAGCGTCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGa                         <  1:4414450/71‑1 (MQ=255)
                                          aTCGCGCAAATCAAAGAGGCACAAAGCGTCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGa                         <  1:2892693/71‑1 (MQ=255)
                                          aTCGCGCAAATCAAAGAGGCACAAAGCGTCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGa                         <  1:1925375/71‑1 (MQ=255)
                                           tCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGa                         >  1:4577645/1‑70 (MQ=255)
                                            cgcgCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGaaa                       >  1:2856678/1‑71 (MQ=255)
                                            cgcgCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGaaa                       >  1:4071392/1‑71 (MQ=255)
                                            cgcgCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGaaa                       >  1:3660481/1‑71 (MQ=255)
                                            cgcgCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGaaa                       >  1:3740770/1‑71 (MQ=255)
                                            cgcgCAAATCAAAGAGGCAAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAAc                      >  1:1600950/1‑71 (MQ=255)
                                             gcgcAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAAc                      >  1:727206/1‑71 (MQ=255)
                                             gcgcAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAAc                      <  1:4171337/71‑1 (MQ=255)
                                                 aaaTCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGt                  <  1:3854653/71‑1 (MQ=255)
                                                 aaaTCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGt                  <  1:1358695/71‑1 (MQ=255)
                                                 aaaTCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGt                  <  1:424822/71‑1 (MQ=255)
                                                   aTCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGaa                        >  1:3731850/1‑63 (MQ=255)
                                                   aTCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGa                         <  1:425116/62‑1 (MQ=255)
                                                    tCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTg                                          >  1:1663671/1‑44 (MQ=255)
                                                      aaaGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCg              >  1:2785114/1‑70 (MQ=255)
                                                       aaGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGca            >  1:3559020/1‑71 (MQ=255)
                                                        agagGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGca            >  1:2657439/1‑70 (MQ=255)
                                                              acaAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCg      >  1:2629055/1‑70 (MQ=255)
                                                               caAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCt    <  1:3123983/71‑1 (MQ=255)
                                                                 aaGCGTCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTgg  >  1:4915363/1‑71 (MQ=255)
                                                                 aaGCGTCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTg   >  1:1483382/1‑70 (MQ=255)
                                                                 aaGCGTCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTg   >  1:1119359/1‑70 (MQ=255)
                                                                 aaGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTgg  >  1:441423/1‑71 (MQ=255)
                                                                 aaGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTgg  >  1:1291112/1‑71 (MQ=255)
                                                                      |                                                                 
AAGATGGATACTCGCGTGTTGCAGGAACAGCGACTGGAAGCCATCGCGCAAATCAAAGAGGCACAAAGCGCCGTTGCTGCCGCGCAGGCTTTGCTGGAGCAACGACAAAGCGAAACTCGTGCCGCACAGTCGCTGG  >  W3110S.gb/4010043‑4010178

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: