| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,631,436 | G→A | 59.0% | A217A (GCG→GCA) | livH → | leucine/isoleucine/valine transporter subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,631,436 | 0 | G | A | 59.0% | 10.1 / 26.5 | 39 | A217A (GCG→GCA) | livH | leucine/isoleucine/valine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base G (9/7); new base A (9/14); total (18/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 | |||||||||||
AACACCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAACCCCTACATC > minE/2631383‑2631487 | aaCACCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGtgctg > 1:1466967/1‑69 (MQ=255)aaCACCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGtgctg > 1:1778434/1‑69 (MQ=255)aaCACCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGtgct > 1:2354997/1‑68 (MQ=255)aaCACCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGGGGCGGGTGtgctg > 1:392033/1‑69 (MQ=255) cacCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTgctgct > 1:359788/1‑69 (MQ=255) cacCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTgctgct > 1:1441347/1‑69 (MQ=255) cacCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTgctgct > 1:1523458/1‑69 (MQ=255) gTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTc < 1:1570669/69‑1 (MQ=255) gTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTc < 1:1922934/69‑1 (MQ=255) gTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTc < 1:1883463/69‑1 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:2691686/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:3328819/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:387782/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:3017668/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:443530/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:608566/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:3193364/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:78340/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:1886195/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:1467324/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:1436733/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:1011825/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCgggcg < 1:1070223/51‑3 (MQ=255) tGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCAGTGACgggcg < 1:2159344/51‑3 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:525814/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:833459/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:860655/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:3137677/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:2534406/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:2447617/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:2347726/1‑69 (MQ=255) cTTTGTGATTGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:1648825/1‑69 (MQ=255) cTTTGTGAATGGCGCGGCGATGGCGGCAGTGGCGGGCGTGCTGCTCGGTCAGTTCTACGGCGTAATTaa > 1:1134306/1‑69 (MQ=255) tgATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAAcccc < 1:2598829/68‑1 (MQ=255) tgATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAAcccc < 1:40916/68‑1 (MQ=255) tgATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAAcccc < 1:44496/68‑1 (MQ=255) tgATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAAcccc < 1:2434720/68‑1 (MQ=255) ggCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAACCCCTACATc > 1:2781440/1‑69 (MQ=255) ggCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAACCCCTACATc > 1:404267/1‑69 (MQ=255) | AACACCGACCGGGTGATTGCGCTGACCTTTGTGATTGGCGCGGCGATGGCGGCGGTGGCGGGTGTGCTGCTCGGTCAGTTCTACGGCGTCATTAACCCCTACATC > minE/2631383‑2631487 |
| 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: ‑ |