Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,224,919 G→A 34.7% intergenic (‑423/‑123) ymgG ← / → ymgH hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,224,9190GA34.7% 87.6 / 72.4 75intergenic (‑423/‑123)ymgG/ymgHhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/20);  new base A (9/17);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 5.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAAAAAT  >  W3110S.gb/1224855‑1224988
                                                                |                                                                     
ttCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGa                                                                  <  1:4930478/70‑1 (MQ=255)
ttCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGa                                                                  <  1:1263101/70‑1 (MQ=255)
   atcaGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAg                                                               >  1:3175784/1‑70 (MQ=255)
   atcaGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGt                                                              <  1:1715376/71‑1 (MQ=255)
     caGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                     >  1:2250454/1‑62 (MQ=255)
     caGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                     >  1:1506315/1‑62 (MQ=255)
      aGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAg                                                               >  1:1544515/1‑67 (MQ=255)
      aGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtat                                                            <  1:218442/70‑1 (MQ=255)
      aGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtat                                                            >  1:4568988/1‑70 (MQ=255)
      aGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtat                                                            <  1:2638084/70‑1 (MQ=255)
      aGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtat                                                            >  1:3176074/1‑70 (MQ=255)
      aGCAACGATTACATTAGTCATTTTATTTTGCCGACAGCCTCATTGTCGAAAGATAAGCGTACGACAGtat                                                            <  1:316071/70‑1 (MQ=255)
       gCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtatta                                                          <  1:3049876/71‑1 (MQ=255)
       gCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtatta                                                          <  1:2454715/71‑1 (MQ=255)
          aCGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                     >  1:719837/1‑57 (MQ=255)
          aCGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                     >  1:3482500/1‑57 (MQ=255)
           cGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtattat                                                         >  1:3784145/1‑68 (MQ=255)
           cGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtattat                                                         >  1:3144627/1‑68 (MQ=255)
              ttACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGtat                                                            <  1:4181094/62‑1 (MQ=255)
              ttACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAg                                                      >  1:2196734/1‑68 (MQ=255)
              ttACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaaa                                                   <  1:3491279/71‑1 (MQ=255)
              ttACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    >  1:4147203/1‑70 (MQ=255)
              ttACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    >  1:4284174/1‑70 (MQ=255)
               tACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                    >  1:4296597/1‑69 (MQ=255)
               tACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    >  1:2196/1‑69 (MQ=255)
               tACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    >  1:1255447/1‑69 (MQ=255)
                    tAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGt                                              >  1:831736/1‑70 (MQ=255)
                       tCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGta                                                             <  1:2244317/52‑1 (MQ=255)
                       tCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAtt                                          <  1:523411/71‑1 (MQ=255)
                         attttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaaaa                                                  >  1:1101383/1‑61 (MQ=255)
                         attttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaaaa                                                  >  1:3265486/1‑61 (MQ=255)
                         attttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAtttt                                        >  1:2331483/1‑71 (MQ=255)
                         attttattttGCCGACGCCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAttt                                         <  1:2154732/70‑1 (MQ=255)
                          ttttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                    <  1:1327830/58‑1 (MQ=255)
                           tttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAttttt                                       >  1:2084839/1‑70 (MQ=255)
                            ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTa                                     >  1:703546/1‑71 (MQ=255)
                            ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaaaa                                                  >  1:3704006/1‑58 (MQ=255)
                             tattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTa                                     >  1:2600716/1‑70 (MQ=255)
                             tattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTAt                                    >  1:904910/1‑71 (MQ=255)
                             tattttCCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTAt                                    >  1:1862/1‑71 (MQ=255)
                                    ccGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTTATTTTTTATcc                                  <  1:1222095/66‑1 (MQ=38)
                                    ccGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATcc                                  <  1:3157617/66‑1 (MQ=255)
                                      gACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCaa                                >  1:1180776/1‑66 (MQ=255)
                                       aCGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCaa                                >  1:1448035/1‑65 (MQ=255)
                                          gCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTc                       <  1:187100/71‑1 (MQ=255)
                                             tCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTaca                           <  1:1149571/64‑1 (MQ=255)
                                             tCATTGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    <  1:2517534/39‑1 (MQ=255)
                                             tCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAtt                                          <  1:1765857/49‑1 (MQ=255)
                                              cATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAgcg                   >  1:666585/1‑71 (MQ=255)
                                              cATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACtt                        <  1:465344/66‑1 (MQ=255)
                                              cATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACtt                        <  1:3729575/66‑1 (MQ=255)
                                                ttGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTAt                                    >  1:2977909/1‑52 (MQ=255)
                                                ttGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    <  1:869308/36‑1 (MQ=255)
                                                ttGTCGAAAGATAAGCATACGACAGTATTATCAGaa                                                    <  1:1122086/36‑1 (MQ=255)
                                                  gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTc                       >  1:4392336/1‑63 (MQ=255)
                                                  gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACt               >  1:2800913/1‑71 (MQ=255)
                                                  gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACt               >  1:4750343/1‑71 (MQ=255)
                                                  gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACt               >  1:4441965/1‑71 (MQ=255)
                                                   tCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTg              <  1:4912178/71‑1 (MQ=255)
                                                     gAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAtt                                          <  1:4211685/41‑1 (MQ=255)
                                                     gAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAtt                                          <  1:4036897/41‑1 (MQ=255)
                                                     gAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCg            <  1:2199970/71‑1 (MQ=255)
                                                       aaGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCgtg          >  1:1795396/1‑71 (MQ=255)
                                                        aGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCgtg          >  1:3227011/1‑70 (MQ=255)
                                                        aGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTa                                     >  1:914674/1‑43 (MQ=255)
                                                            aaGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaa      >  1:3151692/1‑70 (MQ=255)
                                                             aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaaa    <  1:4121540/71‑1 (MQ=255)
                                                             aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaaa    <  1:1933263/71‑1 (MQ=255)
                                                             aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaa     <  1:4946418/70‑1 (MQ=255)
                                                             aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaa     <  1:3340380/70‑1 (MQ=255)
                                                             aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaa     <  1:1728038/70‑1 (MQ=255)
                                                             aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaa     <  1:105215/70‑1 (MQ=255)
                                                             aGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaa     <  1:4772033/70‑1 (MQ=255)
                                                               cGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAAAAAt  <  1:4679632/71‑1 (MQ=255)
                                                               cGTACGACAGTATTATAAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAAAAAt  <  1:3152453/71‑1 (MQ=255)
                                                                |                                                                     
TTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAAAAAT  >  W3110S.gb/1224855‑1224988

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: