Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,781,385 | 0 | T | G | 53.8% | 10.1 / 13.1 | 13 | G396G (GGT→GGG) | envC | murein hydrolase activator EnvC |
Reads supporting (aligned to +/- strand): ref base T (1/5); new base G (7/0); total (8/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTT > NZ_CP009273/3781297‑3781472 | gCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGC‑GGCCAGCCGATTCCACTGGTGGGCAGCAGTGGCGGTCAGGGTCg < 1:312381/90‑1 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCct > 1:381961/1‑90 (MQ=255) cACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTAtt > 1:34001/1‑90 (MQ=255) gAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAAt > 1:12184/1‑90 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTc > 1:136211/1‑90 (MQ=255) ttGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCAGTGGTGGGCAGCAGTGGCGGGCCGGGGCGGCCTTCCCCCTATTTCGAAAATtgcc > 2:68841/1‑90 (MQ=255) gTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCgccgcc < 1:377729/90‑1 (MQ=255) ggTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGGCTTCCCTCTTTTTTGAAATTTGCCCCCCGGGTCa > 2:127001/1‑90 (MQ=255) gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 1:392032/90‑1 (MQ=255) ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCCTCCCCCTATTTTGAAATTCGCCGCCCGGGGCAGGCGGGCAAt > 2:185157/1‑90 (MQ=255) gCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTgg < 1:323760/90‑1 (MQ=255) gcagcaGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCCGGGGCAGGCGGGCAATCCACAGCCGTGGGTGGGGAGAt > 1:293837/1‑90 (MQ=255) ggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAttt < 2:136211/90‑1 (MQ=255) | GCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTT > NZ_CP009273/3781297‑3781472 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |