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.gb | 373,280 | 0 | A | G | 32.4% | 38.9 / 30.8 | 37 | A63A (GCA→GCG) | mhpE | 4‑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: ‑ |