Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 3,819,131 | 0 | C | G | 56.9% | 19.5 / 65.3 | 51 | intergenic (‑142/‑116) | gmk/ligB | guanylate kinase/DNA ligase, NAD(+)‑dependent |
Reads supporting (aligned to +/- strand): ref base C (4/18); new base G (29/0); total (33/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.62e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCG > W3110S.gb/3819084‑3819192 | tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGGAAt < 1:2278180/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:1144271/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:956637/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:3492984/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:288313/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:2098365/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:1918751/60‑1 (MQ=255) tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt < 1:1792576/60‑1 (MQ=255) gcgcAATGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTAACTGCAATaa < 1:988399/61‑1 (MQ=255) gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa < 1:3097497/61‑1 (MQ=255) gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa < 1:4197838/61‑1 (MQ=255) gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa < 1:2687438/61‑1 (MQ=255) gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa < 1:474770/61‑1 (MQ=255) cgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAg > 1:1456174/1‑62 (MQ=255) cAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGAc > 1:3013843/1‑62 (MQ=255) cAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGAc > 1:2926115/1‑62 (MQ=255) cAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGAc > 1:4330535/1‑62 (MQ=255) ttgttgCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa < 1:335654/53‑1 (MQ=255) gCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATACAGACGATTTTCAGGAAGCGt < 1:1089014/61‑1 (MQ=255) gCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGt < 1:1235007/61‑1 (MQ=255) gCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGt < 1:4366890/61‑1 (MQ=255) cGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGt < 1:2646656/60‑1 (MQ=255) gTAGCTGCTTCCTCATTTGCCTGTGATTCACTGCAATAAAGACg > 1:3508874/1‑44 (MQ=38) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACg > 1:2622653/1‑44 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:553055/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:3105082/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:3537547/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:1189299/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:947686/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:3941448/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:4044024/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:84386/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:656492/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:4647466/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:4379188/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:487358/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:4388242/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:4482382/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:4523667/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:1685809/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:3209070/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:1665868/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:1752751/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:1951153/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:2999852/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:2962097/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:1194838/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:2821064/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:2071781/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:2241661/1‑45 (MQ=255) gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa > 1:3237381/1‑45 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTc > 1:2305222/1‑39 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:2718133/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:2225101/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:213247/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:2080887/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:2718149/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:1988692/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:3711458/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:152128/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:1594118/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:3639992/1‑62 (MQ=255) cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg > 1:119440/1‑62 (MQ=255) | TGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCG > W3110S.gb/3819084‑3819192 |
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: ‑ |