| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,733,821 | 0 | T | G | 55.0% | 16.0 / 22.5 | 21 | T79P (ACC→CCC) | punR | DNA‑binding transcriptional activator PunR |
| Reads supporting (aligned to +/- strand): ref base T (3/6); new base G (11/0); total (15/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.17e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTCTGGCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGT > NZ_CP009273/1733741‑1733898 | ttCTGGCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGGCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGcct > 2:719697/1‑88 (MQ=255) tGGCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCAtttt > 2:290631/1‑90 (MQ=255) gACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGac > 1:422696/1‑90 (MQ=255) gACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGac < 1:525628/90‑1 (MQ=255) atCCACTGCGATAGCTAACTGACCGCGCCAGCCGCTCGCTATCTGCTGCCATTGCTGGCGGGCGAACCGGGTTTTTTTGACAACAAAgcg > 1:478859/1‑90 (MQ=255) tCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTGACAACAGAgcgc > 1:244797/1‑90 (MQ=255) aCTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTGACAACAGAGCGCCct > 1:674705/1‑90 (MQ=255) cTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGGTCTGCATTTTTTTTACAACAAAGCGCCctt > 2:135105/1‑90 (MQ=255) cTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCCTTTTTTTTGCAACAGAGCGCCctt > 2:451768/1‑90 (MQ=255) gCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTTACAACAGAGCGCCcttct > 2:298887/1‑90 (MQ=255) gCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTGTCAACAGAGCGCCcctct > 2:685541/1‑90 (MQ=255) gCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTGACAACAGAGCGCCcttct > 1:444401/1‑90 (MQ=255) gACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTTACAACAGAGCGCCCTTCTTTGAGAAACCACg > 2:273022/1‑90 (MQ=255) gcgcCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTcc < 2:17327/90‑1 (MQ=255) aGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCcagcag < 1:201705/90‑1 (MQ=255) ttCGCTATCTGCTGGCATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTc < 2:29775/90‑1 (MQ=255) ttCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTc < 1:779546/90‑1 (MQ=255) tCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCa > 1:539654/1‑90 (MQ=255) gACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCcac > 1:536414/1‑85 (MQ=255) gACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCcac < 2:536414/85‑1 (MQ=255) aTTGCTGGCGGGGGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGt > 2:111871/1‑90 (MQ=255) | TTCTGGCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGT > NZ_CP009273/1733741‑1733898 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |