Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,560,998 A→G 21.0% F44S (TTT→TCT)  yihX ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,560,9980AG21.0% 75.1 / 17.4 38F44S (TTT→TCT) yihXpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base G (2/6);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCAGGCTC  >  W3110S.gb/3560930‑3561067
                                                                    |                                                                     
gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                     >  1:3690143/1‑71 (MQ=255)
gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                     >  1:2270938/1‑71 (MQ=255)
gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                     >  1:2131865/1‑71 (MQ=255)
        tcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTcccc                                                             >  1:1746512/1‑71 (MQ=255)
        tcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTcccc                                                             >  1:343207/1‑71 (MQ=255)
          aTGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGcc                                                                  <  1:2138419/64‑1 (MQ=255)
          aTGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTccc                                                              <  1:988202/68‑1 (MQ=255)
             acacAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCatat                                                        >  1:397339/1‑71 (MQ=255)
                caGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGaa                                                     <  1:3528573/71‑1 (MQ=255)
                         tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACtctt                                              <  1:2233075/69‑1 (MQ=255)
                         tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACtcttct                                            >  1:578925/1‑71 (MQ=255)
                         tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCACCATATGAAAACtcttct                                            >  1:3933745/1‑71 (MQ=255)
                           cGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTcttctt                                           >  1:94290/1‑70 (MQ=255)
                             aaCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTc                                                                >  1:1155330/1‑47 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCTCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:1856553/71‑1 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:3638724/71‑1 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:810590/71‑1 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:777718/71‑1 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:1341526/71‑1 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       >  1:3496410/1‑71 (MQ=255)
                              aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:1884201/71‑1 (MQ=255)
                                    cgtcgCTAACTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCata                                                         >  1:3141158/1‑47 (MQ=255)
                                     gtcgCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCAAATGAAAACTCTTCTTAAGCGATGc                                  >  1:3970517/1‑69 (MQ=255)
                                     gtcgCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCAAATGAAAACTCTTCTTAAGCGATGc                                  >  1:1572347/1‑69 (MQ=255)
                                          tAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                     <  1:3270131/45‑1 (MQ=255)
                                            aTTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaaaa                                                   >  1:3860881/1‑45 (MQ=255)
                                            aTTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                          <  1:455638/70‑1 (MQ=255)
                                            aTTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                          <  1:176908/70‑1 (MQ=255)
                                            aTTTCCCCACGCTCATGCTAATGAAACGCCTCCCCCATATGaaaa                                                   >  1:1876443/1‑45 (MQ=255)
                                             tttCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:3765549/56‑1 (MQ=255)
                                             tttCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:1432496/56‑1 (MQ=255)
                                             tttCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                     <  1:1274275/42‑1 (MQ=255)
                                              ttCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATAc                       >  1:1708018/1‑71 (MQ=255)
                                                ccccACGCGCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                       <  1:2956816/53‑1 (MQ=255)
                                                         cATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTaaa               >  1:2012823/1‑68 (MQ=255)
                                                            gCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTc         >  1:504663/1‑71 (MQ=255)
                                                                  gaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCAGGCt   >  1:1953598/1‑71 (MQ=255)
                                                                  gaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCAGGCt   >  1:1703438/1‑71 (MQ=255)
                                                                   aaaCGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCAGGCTc  >  1:1049103/1‑71 (MQ=255)
                                                                    |                                                                     
GAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCAGGCTC  >  W3110S.gb/3560930‑3561067

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: