| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,657,262 | C→T | 50.0% | E273E (GAG→GAA) | wzyE ← | predicted Wzy protein involved in ECA polysaccharide chain elongation |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,657,262 | 0 | C | T | 50.0% | ‑5.0 / 23.4 | 18 | E273E (GAG→GAA) | wzyE | predicted Wzy protein involved in ECA polysaccharide chain elongation |
| Reads supporting (aligned to +/- strand): ref base C (7/2); new base T (4/5); total (11/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.35e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCA > W3110S.gb/3657196‑3657326 | gcgGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATtctc > 1:195871/1‑69 (MQ=255) gACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACg > 1:158866/1‑71 (MQ=255) cAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGc < 1:1315417/71‑1 (MQ=255) gCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgc > 1:831399/1‑71 (MQ=255) cTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAg > 1:644997/1‑71 (MQ=255) gAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga < 1:694797/71‑1 (MQ=255) tCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTgaga < 1:206921/69‑1 (MQ=255) tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa > 1:1127849/1‑59 (MQ=255) tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga < 1:441225/71‑1 (MQ=255) tgttgtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAgaga < 1:1697292/71‑1 (MQ=255) gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACg > 1:361388/1‑71 (MQ=255) tCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGaa < 1:1349111/67‑1 (MQ=255) tGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGgt > 1:1704041/1‑38 (MQ=255) tGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTcgcg < 1:323375/44‑1 (MQ=255) aGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATc > 1:1755480/1‑71 (MQ=255) gCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa > 1:1109412/1‑46 (MQ=255) gCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATa > 1:588924/1‑46 (MQ=255) cAACGCCAGATTTTCCCACG‑‑GAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGaa > 1:448651/1‑58 (MQ=255) cAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTca > 1:515695/1‑71 (MQ=255) | GCGGACAATTGGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCATCGCCGCTCA > W3110S.gb/3657196‑3657326 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |