| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,906,011 | 0 | T | G | 54.5% | 11.1 / 13.5 | 11 | T275P (ACC→CCC) | htpX | protease HtpX |
| Reads supporting (aligned to +/- strand): ref base T (1/4); new base G (6/0); total (7/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905922‑1906092 | cAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtgggt < 1:597419/90‑1 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa < 1:181363/90‑1 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa < 1:593072/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAAATCACTGAGCGATTTTGACTTACCGTTAATGCAGAGGGcc > 2:462644/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTTATGCAGAGAGcc > 1:367905/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAACGCAATGTGCGGTTTTGTCTTTCCGGTAATTGAGAGAGCcctcattacg > 1:527313/1‑86 (MQ=255) tCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTgggt > 2:873446/1‑88 (MQ=255) tCGTTTATCCAGCGGCGGGGGGGGGATGAACAAATCCCTGAGCGATTTCGACTTTCCGTTAATGCAGAGAGACCTCATGCTGGGTGgttt > 1:350874/1‑89 (MQ=255) gTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt > 2:252823/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGGGATGAAAAAACCACTCAGCGATTTTGCCTTTCCCTGTATGGAGAGGGCCATCACGCTGGGTGcttctt > 1:455772/1‑90 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:252823/90‑1 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 2:527313/90‑1 (MQ=255) gggtgggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa > 1:591042/1‑46 (MQ=255) gggtgggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa < 2:591042/46‑1 (MQ=255) | CAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905922‑1906092 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |