Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,790,961 | C→T | 100% | D269D (GAC→GAT) | gabD → | succinate‑semialdehyde dehydrogenase I, NADP‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,790,961 | 0 | C | T | 77.5% | 79.4 / 17.4 | 62 | D269D (GAC→GAT) | gabD | succinate‑semialdehyde dehydrogenase I, NADP‑dependent |
Reads supporting (aligned to +/- strand): ref base C (13/1); new base T (35/13); total (48/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.59e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATTCCGC > minE/1790911‑1791003 | tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGaa > 1:286710/1‑69 (MQ=255) tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGaa > 1:1160102/1‑69 (MQ=255) tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGaa > 1:1180497/1‑69 (MQ=255) tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGaa > 1:1460439/1‑69 (MQ=255) tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGaa > 1:1570853/1‑69 (MQ=255) tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGaa > 1:829803/1‑69 (MQ=255) tCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGa > 1:845870/1‑68 (MQ=255) tGGAGCTGGGCGGCTACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:2097898/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:2599106/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:2504501/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:625603/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:806034/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:1896726/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:1610090/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:3042160/1‑69 (MQ=255) tGGAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:975901/1‑69 (MQ=255) ggAGCTGGGCGGTAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:393170/1‑68 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcgc > 1:919578/1‑69 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcgc > 1:770560/1‑69 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcgc > 1:3118033/1‑69 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcgc > 1:2887420/1‑69 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcgc > 1:2379560/1‑69 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcgc > 1:3099176/1‑69 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGgcg > 1:408654/1‑68 (MQ=255) ggAGCTGGGCGGCAACGCACCGTTTATTGTCTTTGACGATGCCGATCTAGATAAAGCCGTGGAAGgcgc > 1:1922696/1‑69 (MQ=255) aaCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaaa > 1:3267708/1‑69 (MQ=255) aaCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaaa > 1:2921916/1‑69 (MQ=255) aaCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaaa > 1:2219268/1‑69 (MQ=255) aaCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaaa > 1:2591422/1‑69 (MQ=255) aaCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaaa > 1:754711/1‑69 (MQ=255) aaCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaa > 1:1302270/1‑68 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:3318204/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:565898/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:3272219/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:719934/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:2893443/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:998311/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:307270/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:1030367/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:1392563/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:1924411/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:2435857/69‑1 (MQ=255) cgcaCCGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCCCTGGCCTCGAAAtt < 1:2152383/69‑1 (MQ=255) gcaccGTTTATTGTCTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt < 1:2294241/65‑1 (MQ=255) ccGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATTCCGc < 1:805299/69‑1 (MQ=255) cTTTGACGATGCCTATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:944581/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCg > 1:1058242/1‑49 (MQ=39) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:1056930/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2012657/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2136553/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2202319/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2240253/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2288381/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:563526/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:454292/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:412916/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2580538/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:3334390/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2720145/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:2963559/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:3181/1‑50 (MQ=255) cTTTGACGATGCCGATCTCGATAAAGCCGTGGAAGGCGCGCTGGCCTCGa > 1:3021882/1‑50 (MQ=255) | TCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATTCCGC > minE/1790911‑1791003 |
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: ‑ |