| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,265,887 | 0 | T | G | 53.8% | 11.3 / 16.7 | 13 | S163R (AGT→AGG) | pdeC | c‑di‑GMP phosphodiesterase PdeC |
| Reads supporting (aligned to +/- strand): ref base T (0/6); new base G (7/0); total (7/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGG > NZ_CP009273/4265801‑4265975 | cgcgATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTTTTAACCCCCTCTTTTTGggtgag > 2:594264/1‑89 (MQ=255) gcgATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGagtgaag < 1:525832/90‑1 (MQ=255) gcgATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCCCTCTTCTGGagggaag > 1:266857/1‑90 (MQ=255) gATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGagtgaagtg < 1:478041/90‑1 (MQ=255) ggCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATTAACCCCCTCTTTTGGGGGGAAGGGATGTCTGATGACCCGACa > 1:68524/1‑90 (MQ=255) ggCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGGGGGAAGTGATGTCTGATGACCCGACa > 1:495656/1‑90 (MQ=255) aaaaTGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCaa < 2:545312/90‑1 (MQ=255) aaaaTGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCaa < 2:13188/90‑1 (MQ=255) gACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATgggg < 1:271922/90‑1 (MQ=255) agccTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATgggg < 2:528720/88‑1 (MQ=255) ccTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGgtg < 1:256551/90‑1 (MQ=255) ccTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGgtg < 2:540974/90‑1 (MQ=255) gCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGGGGGGAGTGATGTCTGATGACCCGACATTGCAATGGGGGGTGTATGATAc > 1:631297/1‑90 (MQ=255) gCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGGGAAGTGATGTCTGATGACCCGACATTGCAATGGGGGGGGTATGATAc > 1:205702/1‑90 (MQ=255) gCGGGGAAATTATGTGGCGGTTATCAAACCTCTCTTCTGGGGGGAAGGGATGTCTGATGACCCCACATTTCAATGGGGGGTGTATGAtct > 1:528720/1‑88 (MQ=255) gggAAATTATGTGGCGGGTATTAAACCCCTCTTCTTGAGGGGAGGGATGTTTTGTGGCCCGACATTGCAATTGGGGGTGGATGATAAggg > 2:319432/1‑89 (MQ=255) tgtgGCGGTTATCAACCCCCTCTTCTGGGGGGAAGGGGTGTCTGATGACCCGACATTGCAATGGGGGGTGTATGATACGGGGACGAAAAc > 2:181969/1‑90 (MQ=255) gtgaagtgaTGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGAGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCgg > 2:661807/1‑90 (MQ=255) | CGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGG > NZ_CP009273/4265801‑4265975 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |