Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 373,280 A→G 32.4% A63A (GCA→GCG mhpE → 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb373,2800AG32.4% 38.9 / 30.8 37A63A (GCA→GCGmhpE4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (9/3);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GATTCGATTGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTT  >  W3110S.gb/373215‑373345
                                                                 |                                                                 
gATTCGATTGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACAt                                                                >  1:1385345/1‑69 (MQ=255)
gATTCGATTGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAg                                                              >  1:4320714/1‑71 (MQ=255)
 aTTCGATTGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGc                                                             <  1:583307/71‑1 (MQ=255)
 aTTCGATTGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGc                                                             <  1:1452642/71‑1 (MQ=255)
        tGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCtt                                                       >  1:2568260/1‑70 (MQ=255)
        tGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCtt                                                       >  1:3037649/1‑70 (MQ=255)
         gAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTg                                                      <  1:1324915/70‑1 (MQ=255)
          aaGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGaa                                                    >  1:2457141/1‑71 (MQ=255)
            gTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGaa                                                    >  1:194497/1‑69 (MQ=255)
            gTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATg                                                  <  1:3214507/71‑1 (MQ=255)
            gTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATg                                                  <  1:3088691/71‑1 (MQ=255)
                  cACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTg                                                      <  1:1719429/61‑1 (MQ=255)
                        gACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAgcggc                                      <  1:782366/71‑1 (MQ=255)
                         aCGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGaa                                                    >  1:994319/1‑56 (MQ=255)
                          cGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAgc                                         <  1:1338623/66‑1 (MQ=255)
                           ggTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAgcggcggc                                   >  1:784930/1‑71 (MQ=255)
                              ttGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGcaca                                                                 >  1:1005193/1‑38 (MQ=255)
                              ttGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGcaca                                                                 >  1:999514/1‑38 (MQ=255)
                              ttGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGCGGCGGCGGa                                >  1:1029015/1‑71 (MQ=255)
                               tGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGa                                <  1:1236673/70‑1 (MQ=255)
                               tGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGAt                               >  1:4325738/1‑71 (MQ=255)
                                 cAAGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGAGGCGGCGGAtg                              <  1:534079/70‑1 (MQ=255)
                                   aGGTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATgtgg                           >  1:3557729/1‑71 (MQ=255)
                                     gTTCCAGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATgtggt                          >  1:2830534/1‑70 (MQ=255)
                                     gTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATgtggtg                         <  1:3013034/71‑1 (MQ=255)
                                        ccAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGa                                                >  1:2018158/1‑45 (MQ=255)
                                        ccAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGAt                                               >  1:100401/1‑46 (MQ=255)
                                        ccAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGAt                                               >  1:3068988/1‑46 (MQ=255)
                                          aGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAgcgg                                       >  1:2250079/1‑52 (MQ=255)
                                          aGCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCa                    <  1:3665641/71‑1 (MQ=255)
                                          aGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCa                    >  1:425326/1‑71 (MQ=255)
                                           gCTTTAACTATGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATgtgg                           >  1:273678/1‑63 (MQ=255)
                                                aaCTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAgcggc                                      <  1:1154511/47‑1 (MQ=255)
                                                     tGGTTTCGGCGCGCATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATcg         >  1:1292967/1‑71 (MQ=255)
                                                     tGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATcg         >  1:2355421/1‑71 (MQ=255)
                                                        tttCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGa                                <  1:1255306/45‑1 (MQ=255)
                                                            ggCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACgtt  >  1:1175823/1‑71 (MQ=255)
                                                             gcgcACATAGCGACCTTGAATGGATTGAAgcggcggc                                   <  1:2737648/37‑1 (MQ=255)
                                                               gcACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCaaaa            <  1:2863105/58‑1 (MQ=255)
                                                                 |                                                                 
GATTCGATTGAAGTGGCCCACGGCGACGGTTTGCAAGGTTCCAGCTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTT  >  W3110S.gb/373215‑373345

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: