Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,560,998 A→G 29.8% 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,9980AG29.8% 60.6 / 36.6 47F44S (TTT→TCT) yihXpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (19/14);  new base G (7/7);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TAGAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCA  >  W3110S.gb/3560928‑3561062
                                                                      |                                                                
tAGAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaa                                                                  <  1:3732193/71‑1 (MQ=255)
 agagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaaa                                                                 >  1:3401180/1‑71 (MQ=255)
  gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCTTCGCTAATTTCCCCACGCTCATGCTGATGAGAc                                                                <  1:1877925/71‑1 (MQ=255)
    gCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGc                                                              >  1:1585068/1‑71 (MQ=255)
        tctcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTcc                                                          >  1:2451014/1‑71 (MQ=255)
          tcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGcc                                                             <  1:884666/66‑1 (MQ=255)
            aTGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaaa                                                                 >  1:753050/1‑60 (MQ=255)
            aTGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTc                                                           >  1:4629261/1‑66 (MQ=255)
                cacaGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATg                                                  <  1:2154390/71‑1 (MQ=255)
                 acaGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATg                                                  <  1:1847470/70‑1 (MQ=255)
                  caGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGc                                                              <  1:2479069/57‑1 (MQ=255)
                   aGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGaa                                                >  1:1375062/1‑70 (MQ=255)
                     cgcCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGaaa                                               >  1:545137/1‑69 (MQ=255)
                           tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaaa                                                                 >  1:910965/1‑45 (MQ=255)
                           tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaaaa                                              >  1:3332346/1‑64 (MQ=255)
                           tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:1694401/62‑1 (MQ=255)
                            gCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                >  1:3705694/1‑45 (MQ=255)
                            gCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                >  1:1514316/1‑45 (MQ=255)
                             cGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTcttctt                                      >  1:3849061/1‑70 (MQ=255)
                                 cGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                  >  1:1265027/1‑70 (MQ=255)
                                 cGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                  >  1:1090320/1‑70 (MQ=255)
                                  gCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGAc                                                                <  1:4440242/39‑1 (MQ=255)
                                  gCTTCGTCGCTAATTTCCCCACCCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:257919/55‑1 (MQ=255)
                                   cTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGAt                               <  1:1406274/71‑1 (MQ=255)
                                     tcgtcgCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGc                             >  1:3400036/1‑71 (MQ=255)
                                      cgtcgCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGc                             >  1:1628048/1‑70 (MQ=255)
                                       gtcgCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGc                             >  1:1768642/1‑69 (MQ=255)
                                         cgCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGc                         >  1:3851937/1‑71 (MQ=255)
                                          gCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGc                         <  1:3048111/70‑1 (MQ=255)
                                           cTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:4147032/71‑1 (MQ=255)
                                           cTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:884990/71‑1 (MQ=255)
                                              aTTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                     <  1:916061/70‑1 (MQ=255)
                                              aTTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGa                      <  1:2016665/69‑1 (MQ=255)
                                               tttCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                  <  1:3639909/56‑1 (MQ=255)
                                               tttCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:38409/67‑1 (MQ=255)
                                               tttCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:2424160/67‑1 (MQ=255)
                                                ttCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACtcttct                                       >  1:2849044/1‑50 (MQ=255)
                                                 tCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAg                                   >  1:4059446/1‑53 (MQ=255)
                                                   cccACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                     <  1:574140/65‑1 (MQ=255)
                                                      aCGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTa            <  1:4514702/71‑1 (MQ=255)
                                                      aCGCTCATGCTGATGAAACGCCTCCCCCATATGAAAAc                                             >  1:3967362/1‑38 (MQ=255)
                                                      aCGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                     >  1:1784545/1‑62 (MQ=255)
                                                        gCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTaa           <  1:3573876/70‑1 (MQ=255)
                                                         cTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATg                              >  1:1352104/1‑50 (MQ=255)
                                                          tCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAAt         >  1:1184854/1‑70 (MQ=255)
                                                          tCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATc        >  1:1559610/1‑71 (MQ=255)
                                                           cATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                     >  1:114740/1‑57 (MQ=255)
                                                             tGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCt     <  1:3554236/71‑1 (MQ=255)
                                                                    gaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATa                   >  1:1687182/1‑50 (MQ=255)
                                                                    gaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATa                   >  1:2115446/1‑50 (MQ=255)
                                                                    gaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCa  >  1:3034138/1‑67 (MQ=255)
                                                                     aaaCGCCTCCCCCATATGAAAACTCTTCTTAAGCGa                                >  1:457135/1‑36 (MQ=255)
                                                                     aaaCGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATa                   >  1:3105830/1‑49 (MQ=255)
                                                                      |                                                                
TAGAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCA  >  W3110S.gb/3560928‑3561062

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: