| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,558,416 | T→C | 59.2% | G201G (GGT→GGC) | gpsA → | glycerol‑3‑phosphate dehydrogenase (NAD+) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,558,416 | 0 | T | C | 59.2% | 15.7 / 27.8 | 44 | G201G (GGT→GGC) | gpsA | glycerol‑3‑phosphate dehydrogenase (NAD+) |
| Reads supporting (aligned to +/- strand): ref base T (7/11); new base C (14/12); total (21/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.73e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCG > minE/2558351‑2558455 | tACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCg > 1:1142503/1‑69 (MQ=255)tACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCg > 1:1006206/1‑69 (MQ=255)tACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCg > 1:909745/1‑69 (MQ=255)tACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCg > 1:3107881/1‑69 (MQ=255)tACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCg > 1:2285213/1‑69 (MQ=255) aGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCggg < 1:300722/68‑1 (MQ=255) aGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCggg < 1:2399783/68‑1 (MQ=255) cGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTcc > 1:973697/1‑68 (MQ=255) cGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCg > 1:169307/1‑69 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:146161/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:1168969/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:859027/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:783754/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:51249/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:276753/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:268259/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:2674745/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:2132274/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:2077693/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGCTTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:2076349/69‑1 (MQ=255) cGGATTTCATCGGTGTGCAGATTGGCGGTGCGGTGAAAAACGTCATTGCCATTGGCGCGGGGATGTCCg < 1:1235502/69‑1 (MQ=255) cTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCg < 1:662853/69‑1 (MQ=255) tGGCGGCGCGTTGAAAAACGTTATTGCCATTGGTGCggg < 1:2438418/39‑1 (MQ=255) tGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCg < 1:34887/67‑1 (MQ=255) aaaaCGTCATTGGCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:1652191/1‑48 (MQ=255) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATGGGTTTTg > 1:1795249/1‑48 (MQ=38) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCgg > 1:3326106/1‑43 (MQ=38) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:2216837/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:1205480/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:1697609/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:3139948/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:3111683/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:1218836/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:1752983/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:2615569/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:261481/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:2547126/1‑48 (MQ=39) aaaaCGTCATTGCCATTGGCGCGGGGATGTCCGACGGTATCGGTTTTg > 1:1389854/1‑48 (MQ=39) ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg < 1:2640090/42‑1 (MQ=255) ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg < 1:2542503/42‑1 (MQ=255) ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg < 1:745667/42‑1 (MQ=255) ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg < 1:1674984/42‑1 (MQ=255) ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg < 1:1408664/42‑1 (MQ=255) ggTGCGGGGATGGCCGACGGTATCGGTTTTGGTGCGAATgcg < 1:124379/42‑1 (MQ=255) | TACAGCAATCCGGATTTCATTGGCGTGCAGCTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCG > minE/2558351‑2558455 |
| 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: ‑ |