Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,781,385 | 0 | T | G | 51.3% | 43.8 / 56.7 | 39 | G396G (GGT→GGG) | envC | murein hydrolase activator EnvC |
Reads supporting (aligned to +/- strand): ref base T (4/15); new base G (20/0); total (24/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTTC > NZ_CP009273/3781297‑3781473 | gCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTc < 1:1365945/90‑1 (MQ=255) cTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCg < 2:161133/90‑1 (MQ=255) cTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCg < 1:1263973/90‑1 (MQ=255) cTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCg > 1:898856/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCGGGGGGGGCGTTCAct > 1:661681/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCCGGGTTGGCCCTCAct > 2:985900/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTcccc > 2:1029590/1‑87 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCct > 2:162502/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCct > 2:1265709/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCCGGGGCGGGCTTCAct > 2:1263973/1‑90 (MQ=255) gAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTTCGGGGCGGCCCTCCTCCCTTTTCGAAAt > 1:98940/1‑90 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTc < 1:431815/90‑1 (MQ=255) ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTcgccgc > 2:605171/1‑90 (MQ=255) tCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAg < 1:767689/90‑1 (MQ=255) gCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCCGGGGCGGCCCTCCCTCTATTTCGAAATTCGCCCCCAggg > 1:1121678/1‑90 (MQ=255) gCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACCCTATTTTGAAAATCGCCCCCAggg > 1:449747/1‑90 (MQ=255) ttCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAgg < 1:145155/90‑1 (MQ=255) cgcgcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCg < 2:372625/90‑1 (MQ=255) gcgcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCg < 2:511149/49‑1 (MQ=255) gcgcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCg > 1:511149/1‑49 (MQ=255) gcgcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCCCCTTTTTTGAAAATCGCCGCCCGGGGCAGgggg > 1:1281850/1‑90 (MQ=255) gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 1:875978/90‑1 (MQ=255) gggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCaa < 1:1142350/90‑1 (MQ=255) ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTTGCCGCCCGGGGCAGGCGGGCAAt > 2:734543/1‑90 (MQ=255) gCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGGCCTCCCCCTTTTTTGAAATTTGCCCCCCGGGGCAGGCGGGCAATc > 1:1224209/1‑90 (MQ=255) aGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCcac < 1:726823/90‑1 (MQ=255) ttGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGt < 2:1016338/90‑1 (MQ=255) ttGCACTGGTGGGCAGCAGTGGCGGTCAAGGGCGGCCTTCACTCTCTTTCGAAATTCGCCCCCAGGGGCAGGCGGTCAATCCCCAGCCgc > 2:1142350/1‑89 (MQ=255) cTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCCGGGGGAGGCGGTCAATCCACAGCCGTGGGTg > 1:919053/1‑90 (MQ=255) gcagcaGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGAt > 2:602475/1‑90 (MQ=255) gcagcaGTGGCGGTCAGGGGCGGCCTTCACTCTTTTTTGAAATTTGCCGCCAGGGGGAGGCGGGGAAACCACAGCCGTGGGTGGGGAGAt > 1:855028/1‑90 (MQ=255) agTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGGAGATAAGtt > 2:315942/1‑90 (MQ=255) agTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGTGGGTGGGAAGATAAGtt > 2:434840/1‑90 (MQ=255) ggCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAg < 2:1140451/45‑1 (MQ=255) ggCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAg > 1:1140451/1‑45 (MQ=255) ggCGGTCAGGGGGGGGCTTTACTCTATTTCGAAATTTGCCGCCAGGGGCAGGGGGGCAATCCACAACCGCGGGTTGGGAGGGAAGTTTTg > 1:826005/1‑90 (MQ=255) gggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAtt > 2:942176/1‑90 (MQ=255) ggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAttt < 1:869547/90‑1 (MQ=255) gTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTtc < 2:855028/90‑1 (MQ=255) | GCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTTC > NZ_CP009273/3781297‑3781473 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |