| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,151,925 | A→C | 58.5% | P207P (CCA→CCC) | fabD → | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,151,925 | 0 | A | C | 58.5% | 28.8 / 18.0 | 24 | P207P (CCA→CCC) | fabD | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
| Reads supporting (aligned to +/- strand): ref base A (7/3); new base C (6/8); total (13/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGGCGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTG > W3110S.gb/1151858‑1151991 | cggcggGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACccgca < 1:1382329/71‑1 (MQ=255) gcAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACACGCAGCCGACaa < 1:579786/71‑1 (MQ=255) gcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAg < 1:209491/71‑1 (MQ=255) cTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGa < 1:1615003/71‑1 (MQ=255) aCCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCg > 1:1210001/1‑70 (MQ=255) aCCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCg > 1:1468881/1‑70 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:16912/1‑71 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:1023928/1‑71 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:1444971/61‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:1472933/61‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:914494/61‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:1638520/61‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaa < 1:1773510/61‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGa < 1:782483/52‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGa < 1:661103/52‑1 (MQ=255) gTCTCACTGTGCGCTGATGGAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:65723/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:960075/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:949345/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:561141/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:1825673/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:1510436/1‑71 (MQ=255) gTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACg > 1:110848/1‑71 (MQ=255) gtgCGCCGATGAAACCCGCAGCCGACAAACTGGCAGTAGaa > 1:544967/1‑41 (MQ=37) tgCGCTGATGAAACCCTCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAg > 1:530639/1‑71 (MQ=255) aaCCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTtgtg < 1:1168636/71‑1 (MQ=255) | CGGCGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTG > W3110S.gb/1151858‑1151991 |
| 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: ‑ |