| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 438,506 | T→C | 28.9% | intergenic (+7/‑207) | pgm → / → ybfP | phosphoglucomutase/hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 438,506 | 0 | T | C | 28.9% | 26.7 / 10.2 | 31 | intergenic (+7/‑207) | pgm/ybfP | phosphoglucomutase/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (9/13); new base C (3/6); total (12/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
CGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCT > minE/438440‑438568 | cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAtt > 1:3291950/1‑67 (MQ=255)cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAtt > 1:3286793/1‑67 (MQ=255)cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAtt > 1:738533/1‑67 (MQ=255) gCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAttt > 1:297891/1‑67 (MQ=255) gCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTta > 1:269295/1‑68 (MQ=255) tGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATaa‑gaaaaag > 1:3335822/1‑67 (MQ=255) tGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATaa‑gaaaaag > 1:1492657/1‑67 (MQ=255) tGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATaa‑gaaaaag > 1:1109989/1‑67 (MQ=255) gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTtaa > 1:1879168/1‑58 (MQ=255) gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAggg > 1:1357749/1‑69 (MQ=255) gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAggg > 1:1375096/1‑69 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:1257676/69‑3 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:2778334/69‑3 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:795937/69‑3 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:1874269/69‑3 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:1460281/68‑3 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAg < 1:3269572/69‑1 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAg < 1:174884/69‑1 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAg < 1:3026723/69‑1 (MQ=255) ggTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:3118976/67‑3 (MQ=255) aTTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGccc < 1:757561/68‑1 (MQ=255) ttGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCtt < 1:2409325/69‑1 (MQ=255) ttGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCtt < 1:2045583/69‑1 (MQ=255) gTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCttt > 1:1568046/1‑68 (MQ=255) gTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCttt > 1:1372543/1‑68 (MQ=255) gTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCttt > 1:3224179/1‑68 (MQ=255) gTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCttt > 1:2594602/1‑68 (MQ=255) cTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACa > 1:25963/1‑69 (MQ=255) cgTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAAt < 1:3268941/69‑1 (MQ=255) aaCACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGcc < 1:2806864/69‑1 (MQ=255) aaCACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGcc < 1:2568863/69‑1 (MQ=255) cacaTTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCt < 1:259803/68‑1 (MQ=255) aTTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACa < 1:2357238/50‑1 (MQ=255) aTTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACa < 1:1383829/50‑1 (MQ=255) aTTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACa < 1:849871/50‑1 (MQ=255) tttAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACa < 1:2833887/49‑1 (MQ=255) | CGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCT > minE/438440‑438568 |
| 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: ‑ |