| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 108,722 | C→T | 22.2% | D322D (GAC→GAT) | aceE → | pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 108,722 | 0 | C | T | 22.2% | 79.0 / 14.3 | 58 | D322D (GAC→GAT) | aceE | pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding |
| Reads supporting (aligned to +/- strand): ref base C (18/27); new base T (9/4); total (27/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.13e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCC > minE/108656‑108784 | tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACt > 1:934097/1‑68 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACt > 1:2686916/1‑68 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACt > 1:457522/1‑68 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACt > 1:2243531/1‑68 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACt > 1:2999474/1‑68 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACt > 1:2129966/1‑68 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:2550056/69‑1 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:2507608/69‑1 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:2452802/69‑1 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:2956415/69‑1 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:501623/69‑1 (MQ=255)tGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:807400/69‑1 (MQ=255)tGAACTGCTGCGTAAGGATACAAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:3325769/69‑1 (MQ=255) ctgctgCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAg < 1:3433358/69‑1 (MQ=255) ctgctgCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAg < 1:1550578/69‑1 (MQ=255) ctgctgCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAg < 1:2438770/69‑1 (MQ=255) aaGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:1433166/56‑1 (MQ=255) gCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATg < 1:3143785/69‑1 (MQ=255) gCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATg < 1:3122268/69‑1 (MQ=255) gCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATg < 1:2670585/69‑1 (MQ=255) gCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATg < 1:1192213/69‑1 (MQ=255) gCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATg < 1:1132390/69‑1 (MQ=255) aCTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa < 1:362199/39‑1 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGCCGGCGACTACCAg > 1:2698108/1‑42 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGGGACTACCAGACCTTCAAATCGAAAGATGGTGCGTa > 1:361954/1‑68 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTa > 1:1997774/1‑38 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:1301476/1‑69 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:1369641/1‑69 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:3218466/1‑69 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:1843019/1‑69 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:3122267/1‑69 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:2359416/1‑69 (MQ=255) cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc > 1:2417503/1‑69 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:2832595/67‑1 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:1043745/67‑1 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:3128195/67‑1 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:1348829/67‑1 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:2511369/67‑1 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:449882/67‑1 (MQ=255) ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATAGAAAGATGGTGCGTACGtt < 1:2948808/67‑1 (MQ=255) aaCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc < 1:2926386/69‑1 (MQ=255) aaCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc < 1:766574/69‑1 (MQ=255) aCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc < 1:379544/68‑1 (MQ=255) aCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc < 1:490651/68‑1 (MQ=255) aaaCCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATg < 1:2659848/42‑1 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:1122547/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:177847/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:397446/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:1728159/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:3367725/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:2133685/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:2318796/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:2637920/1‑68 (MQ=255) aaCCGTTGACGGCGATTACCAGACATTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg > 1:778783/1‑68 (MQ=255) gTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTaaa < 1:1409904/69‑1 (MQ=255) ggCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt < 1:820222/42‑1 (MQ=255) ggCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATcc > 1:3261032/1‑68 (MQ=255) ggCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATcc > 1:67398/1‑68 (MQ=255) | TGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCC > minE/108656‑108784 |
| 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: ‑ |