| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,368,067 | Δ1 bp | 100% | coding (426/2190 nt) | fadB → | fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,368,067 | 0 | G | . | 80.4% | 125.8 / 40.3 | 51 | coding (426/2190 nt) | fadB | fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base G (4/6); new base . (0/41); total (4/47) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.40e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
CGTCTGGCGACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTA > minE/2368020‑2368097 | cGTCTGGCGACTCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGcatcat > 1:1229058/1‑56 (MQ=255)cGTCTGGCGACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGcatcat > 1:192143/1‑56 (MQ=255)cGTCTGGCGACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGcatcat > 1:1743460/1‑56 (MQ=255)cGTCTGGCGACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGcatcat > 1:396103/1‑56 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:617409/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:240388/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:267513/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:305918/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:307833/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:339617/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:34760/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:444824/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:463501/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:531813/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1924601/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:685203/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:735935/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:746557/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:829843/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:861450/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:889426/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:916973/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:943304/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:946245/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:973770/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1651223/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1190083/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1232369/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1348837/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1355162/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1368544/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1387804/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:139480/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1427599/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:163012/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:115314/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:113883/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:176658/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1803171/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1857668/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1874601/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1879814/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1042889/65‑1 (MQ=255) gACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1037824/65‑1 (MQ=255) gCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTa < 1:1395300/67‑1 (MQ=255) gCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTa < 1:916782/67‑1 (MQ=255) gCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTa < 1:1033453/67‑1 (MQ=255) gCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTa < 1:379058/67‑1 (MQ=255) gCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTa < 1:1727906/67‑1 (MQ=255) gCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTa < 1:1972103/67‑1 (MQ=255) cTGCGCATCGGTCTGCCGGAAACCAAACT‑GGCATCATGCCTGGCTTTGGCGGTTc < 1:1955179/55‑1 (MQ=255) | CGTCTGGCGACGCCGGATCTGCGCATCGGTCTGCCGGAAACCAAACTGGGCATCATGCCTGGCTTTGGCGGTTCTGTA > minE/2368020‑2368097 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |