Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,141,015 | 0 | A | C | 30.9% | 125.6 / 14.1 | 55 | T30P (ACC→CCC) | yijC | DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base A (32/6); new base C (0/17); total (32/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.48e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGC > REL606/4140965‑4141065 | gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:4486877/51‑1 (MQ=255) gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:2889076/51‑1 (MQ=255) gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:4705334/51‑1 (MQ=255) cAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc > 1:6103546/1‑51 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc > 1:4843457/1‑51 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc < 1:6383200/51‑1 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc > 1:3326666/1‑51 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc < 1:3164627/51‑1 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc < 1:7163924/51‑1 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc > 1:4985059/1‑51 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc < 1:3011923/51‑1 (MQ=255) ggATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg < 1:2053364/51‑1 (MQ=255) gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc < 1:1118040/51‑1 (MQ=255) gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc > 1:5908607/1‑51 (MQ=255) aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc > 1:1749153/1‑51 (MQ=255) tCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGt > 1:7158192/1‑51 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGtt > 1:3618202/1‑51 (MQ=255) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTc > 1:2390751/1‑51 (MQ=255) cATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt > 1:6628913/1‑51 (MQ=255) aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg > 1:4492387/1‑51 (MQ=255) tCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg > 1:4031843/1‑51 (MQ=255) gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg > 1:6178353/1‑51 (MQ=255) tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg > 1:2972884/1‑51 (MQ=255) gATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGa > 1:3939862/1‑51 (MQ=255) aTGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGaa > 1:3485353/1‑51 (MQ=255) ggCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc > 1:1362539/1‑51 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:4218898/1‑51 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:1674538/1‑51 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:5751230/1‑51 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg > 1:3935870/1‑38 (MQ=255) agagCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt > 1:1028190/1‑51 (MQ=255) agagCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt > 1:3749466/1‑51 (MQ=255) gagCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTa > 1:5602310/1‑51 (MQ=255) gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 1:4655324/1‑51 (MQ=255) gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 1:6072478/1‑51 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt > 1:1165121/1‑51 (MQ=255) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGcc < 1:1076260/35‑1 (MQ=25) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCAttt < 1:4037848/41‑1 (MQ=37) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:910126/42‑1 (MQ=37) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:5159917/42‑1 (MQ=37) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg < 1:4506949/43‑1 (MQ=37) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa < 1:5132461/47‑1 (MQ=38) aaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGa < 1:609787/31‑1 (MQ=21) caaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa > 1:5778113/1‑51 (MQ=255) aaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:476836/42‑1 (MQ=37) aaaaaaaaaaaaCCCGCCGTTCGCTGGTGGAAGCCg < 1:6258820/36‑1 (MQ=255) aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGcc < 1:5274969/35‑1 (MQ=25) aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGc < 1:7023889/34‑1 (MQ=25) aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCAttt < 1:2888842/41‑1 (MQ=38) aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCAt < 1:4780038/39‑1 (MQ=37) aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:4682201/42‑1 (MQ=38) aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:5179521/42‑1 (MQ=38) aaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTaa > 1:1498713/1‑51 (MQ=255) aaaaaaaaaaCCCCGCCGTTCGCTGGTGGAAGCCg < 1:4502080/35‑1 (MQ=25) aaaaaaaaaaCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAg < 1:1287066/51‑1 (MQ=39) aaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt < 1:6785157/51‑1 (MQ=39) aaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt > 1:1371667/1‑51 (MQ=255) aaaaaaaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg < 1:3152482/48‑1 (MQ=39) aaaaaaaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg < 1:2063240/48‑1 (MQ=39) aGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGc > 1:204476/1‑51 (MQ=255) gAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt > 1:1156895/1‑51 (MQ=255) aaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACg > 1:4657899/1‑51 (MQ=255) aCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc > 1:1471144/1‑51 (MQ=255) | GCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGC > REL606/4140965‑4141065 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |