| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 438,503 | C→T | 28.9% | intergenic (+4/‑210) | pgm → / → ybfP | phosphoglucomutase/hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 438,503 | 0 | C | T | 28.9% | 28.8 / 9.3 | 31 | intergenic (+4/‑210) | pgm/ybfP | phosphoglucomutase/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (13/9); new base T (3/6); total (16/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.52e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCT > minE/438437‑438568 | cATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAAcaca > 1:1749279/1‑68 (MQ=255)cATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAAcaca > 1:743236/1‑68 (MQ=255)cATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAAcaca > 1:604671/1‑68 (MQ=255) cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAtt > 1:3286793/1‑67 (MQ=255) cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAtt > 1:3291950/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:1492657/1‑67 (MQ=255) tGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATaa‑gaaaaag > 1:3335822/1‑67 (MQ=255) tGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATaa‑gaaaaag > 1:1109989/1‑67 (MQ=255) gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTtaa > 1:1879168/1‑58 (MQ=255) gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAggg > 1:1375096/1‑69 (MQ=255) gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAggg > 1:1357749/1‑69 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:795937/69‑3 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:2778334/69‑3 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:1257676/69‑3 (MQ=255) gCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:1874269/69‑3 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:1460281/68‑3 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAg < 1:174884/69‑1 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAg < 1:3026723/69‑1 (MQ=255) cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAg < 1:3269572/69‑1 (MQ=255) ggTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACATATCAAATAAGAAAAAGGGCGGTCgga < 1:3118976/67‑3 (MQ=255) aTTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGccc < 1:757561/68‑1 (MQ=255) ttGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCtt < 1:2045583/69‑1 (MQ=255) ttGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCtt < 1:2409325/69‑1 (MQ=255) gTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCttt > 1:1372543/1‑68 (MQ=255) gTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCttt > 1:1568046/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:2568863/69‑1 (MQ=255) aaCACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGcc < 1:2806864/69‑1 (MQ=255) cacaTTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCt < 1:259803/68‑1 (MQ=255) | CATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAA‑AAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCT > minE/438437‑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: ‑ |