Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,637,539 A→G 13.6% V134A (GTT→GCT)  yfgM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,637,5390AG13.6% 107.8 / 10.1 44V134A (GTT→GCT) yfgMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (3/3);  total (23/21)
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.93e-01

CACCCTTCACCTTTGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGCCCCTGTTGTA  >  W3110S.gb/2637473‑2637608
                                                                  |                                                                     
cACCCTTCACCTTTGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAAcgcg                                                                   <  1:2380774/71‑1 (MQ=255)
   ccTTCACCTTTGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAGCGCGAGc                                                                <  1:511418/71‑1 (MQ=255)
       cACCTTTGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCaa                                                              <  1:3183910/69‑1 (MQ=255)
          cTTTGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACg                                                          >  1:3416495/1‑70 (MQ=255)
             tGATGTTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCa                                                        <  1:3169696/69‑1 (MQ=255)
             tGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAg                                                       <  1:788553/70‑1 (MQ=255)
               aTGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAtt                                                    <  1:4396575/71‑1 (MQ=255)
                tGGTATCAAGGGTTTTCAGCGCGGCATCAGCGTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAttt                                                   <  1:1215011/71‑1 (MQ=255)
                tGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAttt                                                   <  1:1452671/71‑1 (MQ=255)
                tGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAttt                                                   <  1:216254/71‑1 (MQ=255)
                  gTATCAAGGGTTTTCAGCGCGTCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTAt                                                 <  1:4220830/71‑1 (MQ=255)
                     tCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCaa                                                              <  1:1190301/55‑1 (MQ=255)
                      cAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGa                                                            >  1:662162/1‑56 (MQ=255)
                      cAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGa                                                            >  1:1854469/1‑56 (MQ=255)
                      cAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGa                                                            >  1:4550408/1‑56 (MQ=255)
                       aaGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAGCGCGAGCAAGACGCAGATTTATCACg                                             >  1:1524662/1‑70 (MQ=255)
                       aaGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACgg                                            <  1:980298/71‑1 (MQ=255)
                         gggTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCt                                          >  1:4642408/1‑71 (MQ=255)
                                      ggCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAtt                                                    >  1:1181066/1‑48 (MQ=255)
                                         aTCAGCCTGCTTGAGCTGTACCTGTACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGAttt                                 <  1:630182/64‑1 (MQ=255)
                                         aTCAGCCTGCTTGAGCTGTACCTGAACGCGAGCaa                                                              >  1:4645605/1‑35 (MQ=255)
                                         aTCAGCCTGCTAGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCg                           >  1:3481961/1‑70 (MQ=255)
                                            aGCCTGCTTGAGCTGTACCTGCACGCGAGCAAGACGCAGATTTATCa                                               >  1:3783472/1‑47 (MQ=255)
                                            aGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAtt                                                    >  1:1147457/1‑42 (MQ=255)
                                               cTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCgtgt                    >  1:420127/1‑71 (MQ=255)
                                                tGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGAtt                                                    <  1:4064443/38‑1 (MQ=255)
                                                tGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCa                                               >  1:1912843/1‑43 (MQ=255)
                                                tGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTc                   >  1:3964827/1‑71 (MQ=255)
                                                 gCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTga                                      <  1:1489544/51‑1 (MQ=255)
                                                 gCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCt                  >  1:3103621/1‑71 (MQ=255)
                                                  cTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTa                                                  >  1:726027/1‑38 (MQ=255)
                                                  cTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTa                                                  >  1:2306140/1‑38 (MQ=255)
                                                  cTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTg                 >  1:2842935/1‑71 (MQ=255)
                                                    tGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCt                  >  1:2080684/1‑68 (MQ=255)
                                                     gAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGc                <  1:2246451/69‑1 (MQ=255)
                                                          gTACCTGAGCGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGccc          >  1:302814/1‑70 (MQ=255)
                                                          gTACCTGAGCGCGAGCAAGACGCAGATTTATCACGGCTTTAAGATTTTCATCGCTCGTGTCTGCCAGccc          >  1:1666310/1‑70 (MQ=255)
                                                          gTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTga                                      <  1:324709/42‑1 (MQ=255)
                                                           tACCTGAGCGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGcccc         <  1:188992/70‑1 (MQ=255)
                                                           tACCTGAGCGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGcccc         <  1:2499157/70‑1 (MQ=255)
                                                           tACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCa                              >  1:87227/1‑49 (MQ=255)
                                                           tACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGcccc         <  1:469388/70‑1 (MQ=255)
                                                           tACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGcccc         <  1:483601/70‑1 (MQ=255)
                                                             ccTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGCCCCtgt      >  1:1107811/1‑71 (MQ=255)
                                                                 aaCGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGCCCCtgttgt   <  1:735821/70‑1 (MQ=255)
                                                                  aCGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGCCCCTGTTGTa  <  1:4530641/70‑1 (MQ=255)
                                                                  |                                                                     
CACCCTTCACCTTTGATGGTATCAAGGGTTTTCAGCGCGGCATCAGCCTGCTTGAGCTGTACCTGAACGCGAGCAAGACGCAGATTTATCACGGCTTTGAGATTTTCATCGCTCGTGTCTGCCAGCCCCTGTTGTA  >  W3110S.gb/2637473‑2637608

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: