| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 324,514 | 0 | T | G | 53.6% | 20.1 / 42.5 | 28 | T178P (ACC→CCC) | betI | transcriptional regulator BetI |
| Reads supporting (aligned to +/- strand): ref base T (2/11); new base G (15/0); total (17/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.63e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAG > NZ_CP009273/324380‑324650 | cGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAATAACCGGGGGttta > 2:257588/1‑138 (MQ=255)cGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCA‑GGGGATTTACGCGGGgtttt > 1:368716/1‑138 (MQ=255) gACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGGATTAGAGGGGGGTTTAt > 1:35147/1‑139 (MQ=255) gacagggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTg < 1:54195/134‑1 (MQ=255) cGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCC‑GGGGAATAGAGCGGGGGTTTTCCAGCGG‑GTTgt > 2:206820/1‑138 (MQ=255) gAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTGC‑GGGGATTAGAGCGGGGTTTTTTCAGCGG‑GTTGCCg > 1:359654/1‑139 (MQ=255) aGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGACA‑GGGAA‑TAGGCCGGGGTTTATACAACGGATTTGcgc > 1:373605/1‑137 (MQ=255) caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCC‑GGGAGTTAGGGCGGGGTTTTTCCAGCGG‑GTTTGCGCTCCGGGcctggc > 2:54195/1‑135 (MQ=255) caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCA‑GGGGATTAACGCGGGGTTTTTTCCGCGG‑TTTTGCGCTCCGAGCcgcgc > 1:224885/1‑139 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCcgcgcg < 1:292416/139‑1 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGGAATAGCGCGGGGGTTTTTCAGCGG‑TTTTGCGCCCCGAGCCgcgcgc > 2:283675/1‑139 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCgcgcgc < 1:283675/139‑1 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCgcgcgc < 2:316906/139‑1 (MQ=255) taAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCa < 1:180324/138‑1 (MQ=255) tgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTTA‑CGGGATTAGCGCCGGGTTTTTTCAGCGG‑CTTCTCCACCACAGCCCCGCGGAGCCATAAACCa > 1:275185/1‑139 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTGC‑GGGGAGTTGCGCGGGGTTTTTTCCGCGG‑TTTGTCCCTCAGACCCGCGCGCCGACCTAATCCctc > 1:28824/1‑139 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCa < 1:257588/137‑1 (MQ=255) cATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCA‑GGGAATTAGAGCGGGGTTTTTTCAGCGG‑GTTTCCGCGCAGAGCCGCGCGGGGCCAAAATACATCAATCAGCGCGGCCAGGcc > 1:103154/1‑139 (MQ=255) tctcCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGccc < 2:268953/139‑1 (MQ=255) ctcctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGGA‑GGGGATAAACGCGGGGTTTTTCCAGCGG‑GTTTGCGCCCAAAAGCGCGGCCAGACCCAATTCATCAACCCACGCGGGCCGGCCGCAGCCCg > 2:358040/1‑139 (MQ=255) ctccgcGTTAAATCGGGGGGTAGATGCTGATTGAAAAAGGGGC‑GGGGCA‑GGGAATTACCGCGTGTTTTATCCAGCGA‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGTCCGTAGCCCg < 1:433824/139‑1 (MQ=255) ctcGGTTAATCGGTGGGTAGTGGCTGAGTGATAAAGGGGC‑GGGTCG‑GGGAAGTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGCGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCg < 1:181880/136‑1 (MQ=255) ctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCg < 1:311688/136‑1 (MQ=255) tGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGAGCGGGGTTTTTTCAGCGG‑GTTTGCGCCCAGAGCCGCGCGCAGGCCTAATTCATCAATCAGCGCGGCCCGGGCGTAGGCCGCCTCCCGCGccccttc > 1:94455/1‑134 (MQ=255) tGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGAC‑GGGAGTTAGGGGGGGGTTTTTACAACGG‑GATTCCCCTCCCAGCCGCCCGCGCCCCGAATACCTCAAATAGCCCGGCCCGGCCGGGGCCCCCCTCCTGTGCCTgcct > 2:57372/1‑136 (MQ=255) aGTGATAAAGTGGC‑GGGGCA‑GGGGATTATCGCGGGGTTTTTTCTCCGG‑GTTGGCCCGCACAGCCCCGCGCCGCCCAAATTCCTCAAACAGCGCGGCCCCGCCGTAACCCCCCTCCCCTGCCTCTTTGcgccgccaatct > 1:60576/1‑133 (MQ=255) aGTGATAAAGTGGC‑GGGGCA‑GGGGAATAGCGCGGGGGTTTTTCACCGG‑TTTGTCGCTCCGAGACGCGCGCGCCCATAAAACCTCAAACAGCGCGGCCCGGCCGGCGGCCGCCTCCTCTGCCTGTTCGCGCGGCCACTcg > 2:76250/1‑139 (MQ=255) gTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCa < 2:359654/139‑1 (MQ=255) ggTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 1:402387/1‑139 (MQ=255) ggTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 1:282947/1‑139 (MQ=255) | CGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCA‑GGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAG > NZ_CP009273/324380‑324650 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |