| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 669,029 | C→T | 30.8% | E37K (GAA→AAA) | ybeA ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 669,029 | 0 | C | T | 30.8% | 33.1 / 18.2 | 26 | E37K (GAA→AAA) | ybeA | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (9/9); new base T (6/2); total (15/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01 | |||||||||||
TGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACC > W3110S.gb/668962‑669093 | tGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCa > 1:851310/1‑69 (MQ=255) cACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCaa < 1:3110911/66‑1 (MQ=255) cACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCa < 1:2296285/65‑1 (MQ=255) tGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAg > 1:2467364/1‑68 (MQ=255) gAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCa < 1:1111903/68‑1 (MQ=255) tGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGaa > 1:3689231/1‑58 (MQ=255) tGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCt > 1:220190/1‑68 (MQ=255) tGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATc > 1:1758431/1‑50 (MQ=255) tGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCtt > 1:1194942/1‑69 (MQ=255) tGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCtt > 1:3204644/1‑69 (MQ=255) gTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAac < 1:2550690/68‑1 (MQ=255) tCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAacga > 1:1040037/1‑69 (MQ=255) cGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAacga > 1:2416637/1‑67 (MQ=255) cGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAacga > 1:3424685/1‑67 (MQ=255) cATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGc < 1:270327/68‑1 (MQ=255) cTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCg > 1:3411697/1‑37 (MQ=255) cTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTa > 1:1791816/1‑69 (MQ=255) cTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGaa < 1:773956/56‑1 (MQ=255) ttCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAAc > 1:854108/1‑69 (MQ=255) ttCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAAc > 1:2661081/1‑69 (MQ=255) gccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGt < 1:3694620/68‑1 (MQ=255) gccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGt < 1:2125844/68‑1 (MQ=255) cggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAacgacg > 1:683077/1‑45 (MQ=255) cggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTg < 1:30877/69‑1 (MQ=255) ggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTg < 1:59861/68‑1 (MQ=255) aTTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTAcc < 1:793483/69‑1 (MQ=255) | TGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACC > W3110S.gb/668962‑669093 |
| 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: ‑ |