Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,927,490 C→T 9.4% H113Y (CAT→TAT)  yobB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,927,4900CT 9.4% 148.6 / 5.2 53H113Y (CAT→TAT) yobBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (3/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAGGTATGGA  >  W3110S.gb/1927421‑1927557
                                                                     |                                                                   
gATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCa                                                                    >  1:922817/1‑71 (MQ=255)
 aTCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt                                                                   <  1:615304/71‑1 (MQ=255)
 aTCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt                                                                   <  1:528390/71‑1 (MQ=255)
    gCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATg                                                                  >  1:2232979/1‑69 (MQ=255)
       ttATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGc                                                              <  1:3365099/70‑1 (MQ=255)
        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGGTCTACCATCAGAGTCATGGTGc                                                              <  1:3674819/69‑1 (MQ=255)
        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGc                                                              <  1:2086083/69‑1 (MQ=255)
        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAt                                                            >  1:1136390/1‑71 (MQ=255)
        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAt                                                            >  1:847286/1‑71 (MQ=255)
        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAt                                                            >  1:1751288/1‑71 (MQ=255)
                    tGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCcgtc                                                >  1:2841839/1‑71 (MQ=255)
                    tGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCcgtc                                                >  1:2837209/1‑71 (MQ=255)
                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGaa                                       >  1:330852/1‑71 (MQ=255)
                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGaa                                       >  1:2127498/1‑71 (MQ=255)
                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGaa                                       >  1:1602468/1‑71 (MQ=255)
                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGa                                        <  1:3077828/70‑1 (MQ=255)
                                 ccGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTg                                                       <  1:1397984/51‑1 (MQ=255)
                                   gTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCa                                    >  1:1302973/1‑68 (MQ=255)
                                   gTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCa                                 <  1:501694/71‑1 (MQ=255)
                                    tGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCa                                 <  1:2961525/70‑1 (MQ=255)
                                    tGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCAc                                <  1:1917853/71‑1 (MQ=255)
                                      gCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACgg                              <  1:2412206/71‑1 (MQ=255)
                                      gCGAAAAACGCCTGGGATCTACAATCAGAGTCATGGTGCATGTCTGGGACGTCGTTCCAGAACCATCACgg                              >  1:874922/1‑71 (MQ=255)
                                        gAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACggtg                            >  1:50938/1‑71 (MQ=255)
                                        gAAAAACGCCTGGGATCTACCATCAGAGTCATGGAGCATGTCTGGGCCGTCGTTCCAGAACCATCACggtg                            >  1:283765/1‑71 (MQ=255)
                                        gAAAAACGCCAGGGATCTACCATCAGAGTTATGGt                                                                <  1:13170/35‑1 (MQ=37)
                                        gAAAAACGCCAGGGATCTACCATCAGAGTTATGGt                                                                <  1:851643/35‑1 (MQ=37)
                                          aaaaCGCCAGGGATCTACCATCAGAGTTATGGTGCATGTCTGGGCCGTCGTTcc                                           >  1:2270927/1‑54 (MQ=255)
                                          aaaaCGCCAGGGATCTACCATCAGAGTTATGGTGCATGTCTGGGCCGTCGTTcc                                           >  1:613548/1‑54 (MQ=255)
                                          aaaaCGCCAGGGATCTACCATCAGAGTTATGGTGCATGTCTGGGCCGTCGTTcc                                           >  1:836011/1‑54 (MQ=255)
                                           aaaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTg                                                       >  1:999990/1‑41 (MQ=255)
                                           aaaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTc                         <  1:1872771/71‑1 (MQ=255)
                                            aaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTc                         >  1:3579386/1‑70 (MQ=255)
                                            aaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCg                        <  1:1085933/71‑1 (MQ=255)
                                            aaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCg                        <  1:444366/71‑1 (MQ=255)
                                              cGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGAt                      <  1:2363272/71‑1 (MQ=255)
                                                 cTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaa                   <  1:3692656/71‑1 (MQ=255)
                                                  tGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaa                   >  1:1990153/1‑70 (MQ=255)
                                                   gggATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaac                  >  1:3628039/1‑70 (MQ=255)
                                                    ggATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaacaa                <  1:1567697/71‑1 (MQ=255)
                                                    ggATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaaca                 >  1:3107815/1‑70 (MQ=255)
                                                         tACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGc            >  1:3190952/1‑70 (MQ=255)
                                                         tACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGc            >  1:1288109/1‑70 (MQ=255)
                                                          aCCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGc            >  1:1069615/1‑69 (MQ=255)
                                                          aCCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGc            >  1:1897558/1‑69 (MQ=255)
                                                           ccATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAg         >  1:1084942/1‑71 (MQ=255)
                                                            cATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAgg        <  1:1687554/71‑1 (MQ=255)
                                                            cATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAgg        <  1:3398786/71‑1 (MQ=255)
                                                              tCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGa                                        <  1:64691/37‑1 (MQ=255)
                                                                agagTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACg                               >  1:1584738/1‑44 (MQ=255)
                                                                  agTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAGGtatgga  <  1:2032584/71‑1 (MQ=255)
                                                                  agTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAGGtatgga  <  1:2135096/71‑1 (MQ=255)
                                                                   gTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAGGtatgga  >  1:760664/1‑70 (MQ=255)
                                                                     |                                                                   
GATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAACAACCGCAAGGTATGGA  >  W3110S.gb/1927421‑1927557

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: