| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,279,224 | G→A | 47.8% | V335I (GTC→ATC) ‡ | yeeJ → | adhesin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,279,224 | 0 | G | A | 47.8% | ‑3.6 / 16.0 | 23 | V335I (GTC→ATC) ‡ | yeeJ | adhesin |
| Reads supporting (aligned to +/- strand): ref base G (3/9); new base A (9/2); total (12/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.23e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01 | |||||||||||
CCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAA > minE/1279159‑1279285 | ccAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCt > 1:78877/1‑69 (MQ=255) gTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtattat > 1:2401024/1‑69 (MQ=255) gTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta > 1:2509561/1‑68 (MQ=255) tttACCTGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg < 1:2665779/54‑1 (MQ=255) tACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTg < 1:1130252/68‑1 (MQ=255) tgCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg < 1:1989603/49‑1 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:820397/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:2092435/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:1980576/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:1811717/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:50163/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:3003956/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:3040986/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:3063262/1‑68 (MQ=255) ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa > 1:3154212/1‑68 (MQ=255) tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:1583538/69‑1 (MQ=255) tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:403228/69‑1 (MQ=255) tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:3095696/69‑1 (MQ=255) tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:3025185/69‑1 (MQ=255) tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:2666087/69‑1 (MQ=255) tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:2542769/69‑1 (MQ=255) aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:2035081/68‑1 (MQ=255) aaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa < 1:847283/67‑1 (MQ=255) | CCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAA > minE/1279159‑1279285 |
| 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: ‑ |