| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,522,711 | C→T | 100% | L131L (TTG→TTA) | pdxB ← | erythronate‑4‑phosphate dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,522,711 | 0 | C | T | 75.7% | 46.3 / 11.2 | 41 | L131L (TTG→TTA) | pdxB | erythronate‑4‑phosphate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base C (8/2); new base T (31/0); total (39/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.49e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CCCACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAA > minE/1522645‑1522770 | cccACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCaa > 1:235482/1‑69 (MQ=255)cccACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCaa > 1:2222831/1‑69 (MQ=255)cccACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCaa > 1:551743/1‑69 (MQ=255)cccACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCaa > 1:3316364/1‑69 (MQ=255)cccACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCaa > 1:480876/1‑69 (MQ=255) cGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTAcc < 1:253607/69‑1 (MQ=255) cGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTAcc < 1:2326346/69‑1 (MQ=255) cGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACg > 1:2822016/1‑68 (MQ=255) cGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACg > 1:3064402/1‑68 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGAcc > 1:1952338/1‑51 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTAc > 1:1227501/1‑59 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACTAACGCCAAc > 1:1331221/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGccaa > 1:1506477/1‑66 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGccaa > 1:481213/1‑66 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGccaa > 1:571084/1‑66 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2945289/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2970343/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2715187/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:3094252/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:3159194/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:319134/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:345596/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:580962/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:882151/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1590387/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1297121/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1427004/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1435703/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:14929/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1545244/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2925274/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1774046/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1798215/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2097926/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2422330/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2651470/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2856503/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2899832/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACAAACGCCAAc > 1:479238/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTACAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1375094/1‑69 (MQ=255) cgcgCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGaaaa > 1:390993/1‑68 (MQ=255) | CCCACGGTCGGCGCGAGGCGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAA > minE/1522645‑1522770 |
| 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: ‑ |