| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 2,620,788 | A→G | 18.2% | E312E (GAA→GAG) | purM → | phosphoribosylaminoimidazole synthetase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,620,788 | 0 | A | G | 18.2% | 70.9 / 11.0 | 33 | E312E (GAA→GAG) | purM | phosphoribosylaminoimidazole synthetase |
| Reads supporting (aligned to +/- strand): ref base A (14/13); new base G (4/2); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01 | |||||||||||
GAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTAT > W3110S.gb/2620720‑2620847 | gAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAg > 1:1933175/1‑70 (MQ=255) tGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAg < 1:2484020/71‑1 (MQ=255) tGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAg < 1:2421608/71‑1 (MQ=255) aaaTGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGa > 1:2664467/1‑64 (MQ=255) tGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGcc < 1:1362820/68‑1 (MQ=255) tGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTc > 1:1244374/1‑71 (MQ=255) gTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc > 1:973325/1‑68 (MQ=255) gTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc > 1:200827/1‑68 (MQ=255) gTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc > 1:934326/1‑68 (MQ=255) gTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc > 1:3260211/1‑68 (MQ=255) tATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCg < 1:3697516/70‑1 (MQ=255) tCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGc < 1:3401604/63‑1 (MQ=255) gCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGccc > 1:1641663/1‑63 (MQ=255) ccTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTc < 1:1415407/71‑1 (MQ=255) ccTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTc < 1:3330667/71‑1 (MQ=255) ccTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTc < 1:471116/71‑1 (MQ=255) tCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAAt > 1:2171956/1‑71 (MQ=255) tCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCaa > 1:3554817/1‑70 (MQ=255) cAACTGTGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATg > 1:3287574/1‑71 (MQ=255) cAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCa < 1:2418405/68‑1 (MQ=255) tCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGaa < 1:1951430/71‑1 (MQ=255) gatTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGt > 1:3618711/1‑71 (MQ=255) attattGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTg > 1:3696671/1‑71 (MQ=255) attGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCTAACGGTGAAAACGCGTGGa > 1:891612/1‑70 (MQ=255) ttGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGaaaa > 1:973777/1‑61 (MQ=255) ttGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGaaa < 1:3309063/71‑1 (MQ=255) ttGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGaa < 1:2263410/70‑1 (MQ=255) gCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGaaaa < 1:154404/70‑1 (MQ=255) cctgcctgCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCg < 1:2755941/71‑1 (MQ=255) tgcctgcTCCGGAAGTGGACAAATCCCTCGCCCTGCTCAATGCCAACGGTGaa < 1:1408303/53‑1 (MQ=255) gcctgcTCCGGAAGTGGACAAAACCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTa > 1:2795074/1‑71 (MQ=255) cctgcTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTat > 1:3373132/1‑71 (MQ=255) cctgcTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTat > 1:2934942/1‑71 (MQ=255) | GAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTAT > W3110S.gb/2620720‑2620847 |
| 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: ‑ |