Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,588,010 T→G 61.5% D64A (GAC→GCC)  ydeP ← predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,588,0100TG61.5% 65.5 / 32.8 44D64A (GAC→GCC) ydePpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base G (20/7);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTA  >  W3110S.gb/1587944‑1588075
                                                                  |                                                                 
aaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCg                                                                 <  1:972464/69‑1 (MQ=255)
aaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCg                                                                 <  1:1792554/69‑1 (MQ=255)
   gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAAt                                                              <  1:1694467/69‑1 (MQ=255)
   gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAAt                                                              <  1:1549631/69‑1 (MQ=255)
   gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAAt                                                              <  1:155610/69‑1 (MQ=255)
           tACCTGCTTATCCGTGACTTCCCAGGCGATTTCTTTTGCGCCGTTTTCACAAATGGCg                                                                 <  1:3554129/58‑1 (MQ=255)
           tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACt                                                      >  1:2006426/1‑69 (MQ=255)
           tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                 >  1:1109298/1‑58 (MQ=255)
           tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                 <  1:456885/58‑1 (MQ=255)
              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACtgtg                                                   >  1:107941/1‑69 (MQ=255)
              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACtgtg                                                   >  1:672714/1‑69 (MQ=255)
              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACtgtg                                                   >  1:1268463/1‑69 (MQ=255)
              cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACtgtg                                                   >  1:1566308/1‑69 (MQ=255)
                gCTTATCCGTGACTTCCCAGGTGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACt                                                      >  1:2753244/1‑64 (MQ=255)
                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACtgtg                                                   >  1:2893020/1‑67 (MQ=255)
                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACt                                                      >  1:443591/1‑64 (MQ=255)
                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACt                                                      >  1:871293/1‑64 (MQ=255)
                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACt                                                      >  1:3365272/1‑64 (MQ=255)
                gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACt                                                      >  1:2324615/1‑64 (MQ=255)
                     tCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGa                                            >  1:657219/1‑69 (MQ=255)
                      ccGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGa                                            <  1:3670046/68‑1 (MQ=255)
                      ccGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGa                                            <  1:1241320/68‑1 (MQ=255)
                       cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGa                                            <  1:503348/67‑1 (MQ=255)
                       cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATc                                          >  1:2164812/1‑69 (MQ=255)
                           aCTTCCGAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCg                                                                 <  1:2179089/42‑1 (MQ=255)
                           aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                      <  1:1342995/69‑1 (MQ=255)
                                   ggCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACg                                                          >  1:273607/1‑41 (MQ=255)
                                   ggCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACg                                                          >  1:3439138/1‑41 (MQ=255)
                                    gCGATTGCTTTTGCGCCGTTTTGACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACaa                             >  1:2828708/1‑69 (MQ=255)
                                    gCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACaa                             >  1:1319902/1‑69 (MQ=255)
                                    gCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACaa                             >  1:703538/1‑69 (MQ=255)
                                      gATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAgg                                             >  1:2646755/1‑51 (MQ=255)
                                      gATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAgg                                             >  1:2186905/1‑51 (MQ=255)
                                      gATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAgg                                             >  1:2544457/1‑51 (MQ=255)
                                         tGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCt                                                 >  1:1974198/1‑44 (MQ=255)
                                          gCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCAc                                                       >  1:2665131/1‑37 (MQ=255)
                                                gcgcCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                  <  1:2975752/68‑1 (MQ=255)
                                                 cgcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCGATGCACAACCCGGACAGTCa                 <  1:2693043/68‑1 (MQ=255)
                                                      ttttCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGccc           <  1:1122153/69‑1 (MQ=255)
                                                        ttCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATg                                  >  1:2768837/1‑44 (MQ=255)
                                                         tCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCctct        >  1:3124197/1‑69 (MQ=255)
                                                             aaaTGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCt    <  1:401369/69‑1 (MQ=255)
                                                               aTGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTa  >  1:1669392/1‑69 (MQ=255)
                                                               aTGTCGAAAGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTa  >  1:1888016/1‑69 (MQ=255)
                                                                  |                                                                 
AAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTA  >  W3110S.gb/1587944‑1588075

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: