| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,522,714 | T→A | 100% | R130R (CGA→CGT) | pdxB ← | erythronate‑4‑phosphate dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,522,714 | 0 | T | A | 86.1% | 70.4 / 4.8 | 36 | R130R (CGA→CGT) | pdxB | erythronate‑4‑phosphate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (3/2); new base A (31/0); total (34/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.59e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAA > minE/1522663‑1522770 | 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) | CGGATCGCAAAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAA > minE/1522663‑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: ‑ |