| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,611,069 | A→G | 100% | C223C (TGT→TGC) | eutB ← | ethanolamine ammonia‑lyase, large subunit, heavy chain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,611,069 | 0 | A | G | 80.2% | 48.2 / 5.4 | 35 | C223C (TGT→TGC) | eutB | ethanolamine ammonia‑lyase, large subunit, heavy chain |
| Reads supporting (aligned to +/- strand): ref base A (7/0); new base G (28/0); total (35/0) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGT > minE/1611023‑1611136 | cgcCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGa > 1:584053/1‑56 (MQ=255) acgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAAttt > 1:600452/1‑69 (MQ=255) acgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAAttt > 1:1517998/1‑69 (MQ=255) acgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAAttt > 1:2051415/1‑69 (MQ=255) acgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAAtt > 1:3010642/1‑68 (MQ=255) acgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAAtt > 1:1594466/1‑68 (MQ=255) cgacgGCTCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:422/1‑69 (MQ=255) cgacgGATCGGTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2105075/1‑69 (MQ=39) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGa > 1:838295/1‑58 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGAt > 1:2589391/1‑59 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAAttt > 1:1985584/1‑68 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:497518/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2938789/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:3147674/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:3253424/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:3311079/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:3313810/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:435250/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2709199/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:505689/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:74365/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:825123/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:84554/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:911242/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:953665/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:1131801/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:1331655/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:1744334/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2803169/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:1908383/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:1970415/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2051299/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2167320/1‑69 (MQ=255) cgacgGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACGCAGCCCTGAGTTGGGATGTTGAATTTg > 1:2296128/1‑69 (MQ=255) gACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCaa > 1:3395265/1‑69 (MQ=255) acacAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAg > 1:3023935/1‑69 (MQ=255) acacAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAg > 1:127477/1‑69 (MQ=255) acacAGCCCTGAGTTGGGATGTTGAAATTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAg > 1:1525118/1‑69 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATa > 1:3396690/1‑43 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATa > 1:2954622/1‑43 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATa > 1:1711011/1‑43 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATa > 1:1625498/1‑43 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATa > 1:1268656/1‑43 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATa > 1:1367999/1‑43 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGt > 1:1903935/1‑69 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGt > 1:1488107/1‑69 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGt > 1:1096701/1‑69 (MQ=255) cacaGCCCTGAGTTGGGATGTTGAATTAGTCGATCACGCCAt > 1:1366061/1‑42 (MQ=255) | CGCCACGACGGATCGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGT > minE/1611023‑1611136 |
| 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: ‑ |