| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,621,984 | 0 | A | C | 52.9% | 14.2 / 24.5 | 17 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
| Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (0/9); total (7/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.11e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT > NZ_CP009273/2621910‑2622062 | caATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGcccc > 1:215670/1‑90 (MQ=255) aaTCTGCCAGAATATACGCTGGAAAATAAAACGGGTATCCCCATGTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCa < 1:276417/90‑1 (MQ=255) tcggCCAAAATATACGCTGGGAAATTAAACGGGTACCCCCATGTAAAAGTTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAAc < 1:857184/87‑1 (MQ=255) aaaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTa > 1:204373/1‑90 (MQ=255) gggggTACCCCCATGTAAAACTTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc < 1:374416/88‑1 (MQ=255) ccGGTTATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc < 2:912929/89‑1 (MQ=255) aCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc > 1:648985/1‑90 (MQ=255) gTATCCCCATGTGAAATGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt < 2:373000/90‑1 (MQ=255) ccATGTAAAACTTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt < 2:648985/90‑1 (MQ=255) cATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACag < 2:158465/53‑1 (MQ=255) cATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACag > 1:158465/1‑53 (MQ=255) aaCGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc < 1:67288/90‑1 (MQ=255) cGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt > 2:867602/1‑90 (MQ=255) cgcgGGGGAATAAACCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc < 1:868600/90‑1 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 2:685859/1‑90 (MQ=255) gACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc < 1:917386/66‑1 (MQ=255) gACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc > 2:917386/1‑66 (MQ=255) | CAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT > NZ_CP009273/2621910‑2622062 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |