| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,102,506 | 0 | T | G | 44.8% | 32.7 / 10.9 | 29 | G88G (GGT→GGG) | yiiR | DUF805 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base T (1/15); new base G (13/0); total (14/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.06e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGG > NZ_CP009273/4102419‑4102595 | ggCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGaa > 1:683248/1‑90 (MQ=255)ggCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGaa > 2:400840/1‑90 (MQ=255) gCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAAc < 2:362865/90‑1 (MQ=255) aGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAATGGGGGGtttttc > 2:68241/1‑88 (MQ=255) cTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCggg < 2:546403/90‑1 (MQ=255) cTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAAGGGGGGGTTTTTCCcgg > 1:440877/1‑90 (MQ=255) tGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGCTGGC‑GGGGGACTGGGGGGTTTTTTCggggg > 2:665577/1‑88 (MQ=255) tGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGCTGGC‑GGGGAACTGGGGGATTTTTCCcgggg > 2:911323/1‑88 (MQ=255) tGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAATGGGGGGTTTTTACggggg > 2:424774/1‑88 (MQ=255) cATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTc < 1:854648/90‑1 (MQ=255) aTGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCt < 2:683248/90‑1 (MQ=255) aTCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGc < 1:400840/90‑1 (MQ=255) ttCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGGAATGGGGGGATTTTTCCGGTGGGTT‑GGAAATgggggg > 2:600497/1‑90 (MQ=255) gcATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAATGGGGGGTTTTACCCGGGGTCTGGGCAAgggggggggggggg > 1:311363/1‑88 (MQ=255) tttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCg < 2:995716/90‑1 (MQ=255) tttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCg < 1:68241/90‑1 (MQ=255) ttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGa < 2:46027/45‑1 (MQ=255) ttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGa > 1:46027/1‑45 (MQ=255) ttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAATGGGGGGTTTTTTCCGGGGTGT‑GGGGAAGGGGGGGGGGGGGGTTTTTccccg > 2:751184/1‑88 (MQ=255) ttCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGGAAATGGGGGGTTTTTCCCGGGGTCTGGGAAAGGGGGGGGGGGGGGTTTTTTccca > 2:852940/1‑89 (MQ=255) cTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACg < 1:600497/90‑1 (MQ=255) gatgatTGTGGCGTGGATGCTGCTGGC‑GGGGAACTGGGGGGTTTTTCCGGGGGGCT‑GGCAAAGGGGGGGGGGGGGGTTTTTCCCGGCGtt > 2:162094/1‑90 (MQ=255) atTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGAtt < 1:305789/90‑1 (MQ=255) tGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGtgatg < 2:517583/90‑1 (MQ=255) ggATGCTGCTGGC‑GGGGAACTGGGGGGTTTTTCCGGGGGTCT‑GGCAAAGGGGGGGGGGGGGGTTTTTCCCCACGGTGATTTTGGgggtga > 1:966497/1‑90 (MQ=255) aTGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGatgatg < 1:1009883/90‑1 (MQ=255) ctgGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATc < 1:424774/90‑1 (MQ=255) ctgGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATc < 1:162094/90‑1 (MQ=255) ctgGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATc < 2:785816/90‑1 (MQ=255) tAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTgg < 1:13649/90‑1 (MQ=255) | GGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGC‑GGGTAACTGGGCGATTTTACCGGGTGTCT‑GGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGG > NZ_CP009273/4102419‑4102595 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |