| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 360,129 | 0 | T | G | 57.1% | 7.4 / 10.8 | 14 | intergenic (‑60/+226) | lacZ/BW25113_RS23320 | beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base T (0/6); new base G (8/0); total (8/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.33e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC > NZ_CP009273/360052‑360193 | ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGcccccccc > 1:18739/1‑88 (MQ=255) gCAGGGGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCTcccccccccccccc > 1:301570/1‑90 (MQ=255) ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCTGCccccccccccccccg > 2:182839/1‑87 (MQ=38) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCccgccgct > 2:158086/1‑88 (MQ=255) tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCCCCCCCCCCCCCCCCGGCCggt > 1:109417/1‑88 (MQ=255) gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCCGCCCGTGCTTTGCCAAGTTc > 1:241384/1‑90 (MQ=19) cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGGTGCTTTTCAAAGGTCaa > 1:378262/1‑90 (MQ=21) agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCCCCCGGCCGGTGCTTTGCAACCTTTAAAAcccccc > 2:143801/1‑86 (MQ=12) tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc < 2:100139/90‑1 (MQ=255) tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc < 2:241384/90‑1 (MQ=255) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:477151/90‑1 (MQ=25) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 1:474252/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:150606/90‑1 (MQ=35) aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:186508/90‑1 (MQ=35) | CCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC > NZ_CP009273/360052‑360193 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |