| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,522,738 | C→A | 100% | V122V (GTG→GTT) | pdxB ← | erythronate‑4‑phosphate dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,522,738 | 0 | C | A | 81.3% | 49.5 / 8.1 | 32 | V122V (GTG→GTT) | pdxB | erythronate‑4‑phosphate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base C (4/2); new base A (26/0); total (30/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.02e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
AAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGC > minE/1522672‑1522795 | aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACTAACGCCAAc > 1:1331221/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2899832/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:882151/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:580962/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:345596/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:319134/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:3159194/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:3094252/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2970343/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2945289/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2925274/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:1297121/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2856503/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2715187/1‑69 (MQ=255)aaGTAAGGTTTTGATCCCTAAAGCTTCCAGTCGCGCCTGTAAACGACGACCAACGTTACCAACGCCAAc > 1:2422330/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:1590387/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:2651470/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) gTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATc > 1:983594/1‑67 (MQ=255) gTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCAATCGCGTTCGGCAAGCATc > 1:944134/1‑67 (MQ=255) ttACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATcagc < 1:1643048/69‑1 (MQ=255) ttACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATcagc < 1:946454/69‑1 (MQ=255) | AAGTAAGGTTTTAATCCCTAACGCTTCCAGTCGCGCCTGCAATCGACGTCCAACGTTACCAACGCCCACGATCCCCACCGTACGGTCGTACAGTGAAAATCCATCGCGTTCGGCAAGCATCAGC > minE/1522672‑1522795 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |