| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 165,787 | G→A | 100% | A563A (GCC→GCT) | glnD ← | uridylyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 165,787 | 0 | G | A | 77.9% | 31.9 / 4.2 | 27 | A563A (GCC→GCT) | glnD | uridylyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (5/1); new base A (0/21); total (5/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.43e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
TTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCG > minE/165729‑165844 | tttGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCg < 1:314855/69‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:2839235/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:943188/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:772134/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:647847/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:530697/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:496066/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:3416903/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:3207330/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:3162991/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:2923847/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1179255/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:2835459/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:2729128/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:2289774/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:2272370/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1878615/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1613828/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1526596/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1380794/67‑1 (MQ=255) tcGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1367481/67‑1 (MQ=255) gaCTTCCGGGTCCTGAATATCGCGGCGCTGAGCGGTCACCGACATCAGCa < 1:1902341/50‑1 (MQ=255) cTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCtg > 1:3177712/1‑69 (MQ=255) cTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCtg > 1:430529/1‑69 (MQ=255) tCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTcg > 1:255050/1‑69 (MQ=255) tCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTcg > 1:2273689/1‑69 (MQ=255) tCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTcg > 1:3346643/1‑69 (MQ=255) | TTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCG > minE/165729‑165844 |
| 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: ‑ |