| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,522,720 | T→A | 100% | G128G (GGA→GGT) | pdxB ← | erythronate‑4‑phosphate dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,522,720 | 0 | T | A | 83.8% | 63.4 / 7.2 | 37 | G128G (GGA→GGT) | pdxB | erythronate‑4‑phosphate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (4/2); new base A (31/0); total (35/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.25e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-02 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCG > minE/1522663‑1522783 | 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:571084/1‑66 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGccaa > 1:481213/1‑66 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGccaa > 1:1506477/1‑66 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:3094252/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2899832/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2925274/1‑69 (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: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:2856503/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2651470/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2422330/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2097926/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1798215/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1774046/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1590387/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1545244/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:14929/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1435703/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1427004/1‑69 (MQ=255) aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1297121/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) gACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCg > 1:3044599/1‑68 (MQ=255) | CGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCG > minE/1522663‑1522783 |
| 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: ‑ |