Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 320,985 A→G 31.2% I52V (ATC→GTC)  ykgE → predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb320,9850AG31.2% 37.2 / 25.6 32I52V (ATC→GTC) ykgEpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (5/5);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGTCGCGTAAATTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATT  >  W3110S.gb/320919‑321053
                                                                  |                                                                    
ggCTGTCGCGTAAATTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATc                                                                    <  1:3128124/69‑1 (MQ=255)
ggCTGTCGCGTAAATTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCa                                                                   >  1:3435041/1‑70 (MQ=255)
        cgTAAATTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGa                                                          <  1:4002617/71‑1 (MQ=255)
            aaTTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTc                                                       >  1:3679713/1‑70 (MQ=255)
            aaTTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTc                                                       >  1:336154/1‑70 (MQ=255)
                  ccGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCGGGATg                                               <  1:1551044/71‑1 (MQ=255)
                    ggAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCaa                                                                  <  1:965841/51‑1 (MQ=255)
                     gagaAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGa                                              >  1:723157/1‑70 (MQ=255)
                       gaAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGaaaa                                           <  1:4881305/71‑1 (MQ=255)
                       gaAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGaaaa                                           <  1:4238419/71‑1 (MQ=255)
                       gaAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGaaaa                                           <  1:2750917/71‑1 (MQ=255)
                          aCAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATc                                        <  1:1004657/71‑1 (MQ=255)
                               gATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGa                                     >  1:5022302/1‑69 (MQ=255)
                               gATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGa                                     >  1:3914754/1‑69 (MQ=255)
                                aTGCTGCGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATcg                                  <  1:3655279/71‑1 (MQ=255)
                                  gctgcGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGa                                                          >  1:3763139/1‑45 (MQ=255)
                                  gctgcGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGa                                                          >  1:1864117/1‑45 (MQ=255)
                                      cGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATcgccgc                              <  1:904546/69‑1 (MQ=255)
                                      cGGTCAGCCTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATcgccgc                              <  1:4646654/69‑1 (MQ=255)
                                      cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATcgccgc                              <  1:2026562/69‑1 (MQ=255)
                                      cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCAc                            >  1:4221711/1‑71 (MQ=255)
                                         tCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATg                                               <  1:3551376/49‑1 (MQ=255)
                                         tCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATg                                               <  1:1690287/49‑1 (MQ=255)
                                         tCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATcgc                                 >  1:901895/1‑63 (MQ=255)
                                          cAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTg                          >  1:707840/1‑69 (MQ=255)
                                              cTGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGAAAAATCTg                                      >  1:4618751/1‑53 (MQ=255)
                                               tGCGATCAATAGCGGTTATGTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGAt                    >  1:2191141/1‑70 (MQ=255)
                                               tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAggg                                                  <  1:945615/40‑1 (MQ=255)
                                                 cGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAac                 >  1:4905873/1‑71 (MQ=255)
                                                     cAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTggagg                      >  1:1340248/1‑62 (MQ=255)
                                                          gCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATcc         >  1:1500884/1‑70 (MQ=255)
                                                             gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGa                                     >  1:853602/1‑39 (MQ=255)
                                                               tatgTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCatta    <  1:3004184/70‑1 (MQ=255)
                                                               tataTCAAAGAAGCGATTCCAGGGATGAAAAATCTGAt                                    >  1:473340/1‑38 (MQ=255)
                                                               tataTCAAAGAAGCGATTCCAGGGATGAAAAATCTGAt                                    >  1:3741420/1‑38 (MQ=255)
                                                                atatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCattatt  >  1:3234925/1‑71 (MQ=255)
                                                                atatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCattatt  >  1:1964212/1‑71 (MQ=255)
                                                                  |                                                                    
GGCTGTCGCGTAAATTTCCCGGAGAAACAGGGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATT  >  W3110S.gb/320919‑321053

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: