Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 1,157,021 | 0 | G | C | 73.6% | 121.2 / 49.9 | 72 | A107A (GCG→GCC) | tmk | thymidylate kinase |
Reads supporting (aligned to +/- strand): ref base G (5/14); new base C (0/53); total (5/67) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.31e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAT > W3110S.gb/1156975‑1157075 | aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2708415/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1853681/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1872965/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1961219/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:203177/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2066051/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2078770/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1035904/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:214545/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2187698/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2212505/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2502867/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2559430/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2614016/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:268098/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:211117/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2715471/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2743403/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2873608/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:360822/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:364498/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:436927/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:58640/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:587018/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:59086/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:813543/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:859246/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:954377/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:984111/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1390686/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1098367/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1098988/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1162980/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1206543/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1226990/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1272960/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1308487/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:131409/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1343814/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1387430/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:179723/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1436029/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1749001/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1748495/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1725947/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1721534/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1716729/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1612303/50‑1 (MQ=255) aCGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1619913/50‑1 (MQ=255) cGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:748990/49‑1 (MQ=39) aCCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:2385936/45‑1 (MQ=38) gTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:134155/39‑1 (MQ=37) gTGATTGGCGATCGCCACGATCTCTCCACTCAAGCCTAt < 1:1146643/39‑1 (MQ=37) tGATTGGCGATCGCCACGATCTGTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:928528/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:42628/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:1137135/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:1561469/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:63736/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:2540694/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:1486596/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:374192/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:2171759/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:1671561/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:327539/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:1662836/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:274617/60‑1 (MQ=255) tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTg < 1:2160819/60‑1 (MQ=255) cTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAt > 1:2280766/1‑62 (MQ=255) cTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAt > 1:2478606/1‑62 (MQ=255) cTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAt > 1:2050768/1‑62 (MQ=255) cTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAt > 1:126322/1‑62 (MQ=255) cTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAt > 1:2560636/1‑62 (MQ=255) | ACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGAT > W3110S.gb/1156975‑1157075 |
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: ‑ |