| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 280 | 0 | G | C | 32.0% | 98.7 / 11.7 | 50 | intergenic (+25/‑56) | thrL/thrA | thr operon leader peptide/bifunctional aspartokinase I/homeserine dehydrogenase I |
| Reads supporting (aligned to +/- strand): ref base G (18/16); new base C (0/16); total (18/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.31e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTA > REL606/232‑329 | aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGcc < 1:4949492/51‑1 (MQ=255)aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGcc < 1:2951989/51‑1 (MQ=255)aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGcc > 1:3941296/1‑51 (MQ=255) ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGccc < 1:2658110/51‑1 (MQ=255) caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc < 1:1434250/50‑1 (MQ=255) caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCa < 1:5786272/51‑1 (MQ=255) caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCa < 1:5242951/51‑1 (MQ=255) aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCa < 1:4136021/50‑1 (MQ=255) tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTg < 1:2257813/51‑1 (MQ=255) aaCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGa > 1:6772395/1‑51 (MQ=255) aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGccc > 1:3053694/1‑42 (MQ=255) cGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa < 1:4398648/51‑1 (MQ=255) ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGccc > 1:5164282/1‑40 (MQ=255) ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGccc > 1:3445289/1‑40 (MQ=255) gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGt > 1:6773630/1‑51 (MQ=255) tGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGccc > 1:5173321/1‑38 (MQ=255) gggCTGACGCGTACAGGAAACACAGAAAAAAGccc > 1:853308/1‑35 (MQ=255) gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg > 1:6101648/1‑51 (MQ=255) tGACGCGTACAGGAAACACAGAAAAAAGccc > 1:2126395/1‑31 (MQ=255) tGACGCGTACAGGAAACACAGAAAAAAGccc > 1:3768944/1‑31 (MQ=255) gACGCGTACAGGAAACACAGAAAAAAGccc > 1:2448121/1‑30 (MQ=255) aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttt > 1:7034062/1‑51 (MQ=255) cgcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttt > 1:5249613/1‑51 (MQ=255) gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttttttt > 1:6730903/1‑51 (MQ=255) ggAAACACAGAAAAAAGCCCCCACCCGACAGTGCGGGCtttttttt > 1:5892012/1‑46 (MQ=39) gAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt > 1:6620475/1‑45 (MQ=255) aaaCACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCa < 1:5194686/50‑1 (MQ=255) aaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCaaa < 1:6339212/51‑1 (MQ=255) cacaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAgg < 1:7289033/51‑1 (MQ=255) aaaaaaacCCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGt < 1:1549026/47‑1 (MQ=38) aaaaaaacCCCCCCCCCTGACAGTGCGGGc < 1:1454340/28‑1 (MQ=12) aGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTaa < 1:6516062/51‑1 (MQ=255) aaaaaacCCCCCCCCCTGACAGTGCGGGCtttttt < 1:222322/34‑1 (MQ=16) aaaaaacCCCCCCCCCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAAc < 1:2239609/50‑1 (MQ=255) aaaaagCCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCaaa < 1:1878914/44‑1 (MQ=37) aaaaagCCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCAAAgg < 1:4131685/46‑1 (MQ=38) aaaaagCCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCAAAg < 1:2283026/45‑1 (MQ=38) aaaaagCCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCAAAg < 1:2382594/45‑1 (MQ=38) aaaaacCCCCCCCCCTGACAGTGCGGGCTTTTTTTTCGACCaa < 1:2641055/43‑1 (MQ=25) aaaaacCCCCCCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGt < 1:1708198/47‑1 (MQ=38) aaaaaaGCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCaaa < 1:2011004/44‑1 (MQ=255) aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGt < 1:4416817/47‑1 (MQ=255) aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACg < 1:7198498/51‑1 (MQ=255) aaaaaaCCCCGCCCCTGACAGTGCGGGCTTTTTTTTc < 1:1245357/37‑1 (MQ=25) aaaaaaCCCCGCCCCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTaa < 1:4262305/49‑1 (MQ=38) aaaaaaCCCCGCACCTGACAGTGCGGGCtttttt < 1:4627368/34‑1 (MQ=255) aaaaaaCCCCGCACCTGACAGTGCGGGCTTTTTTTTCg < 1:3214620/38‑1 (MQ=255) aaaaaaCCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTaa < 1:3289407/49‑1 (MQ=255) aaaaaaCCCCCCACCTGACAGTGCGGGCtttttttt < 1:3497373/36‑1 (MQ=37) aGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTa > 1:7145574/1‑51 (MQ=255) | ACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTA > REL606/232‑329 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |