Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,621,984 | 0 | A | C | 61.5% | 6.2 / 10.3 | 13 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/8); total (5/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.77e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT > NZ_CP009273/2621897‑2622062 | gTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc > 2:488302/1‑90 (MQ=255) caATCTGCAAGAATTTGCGCGGGAAAATTAAGCGGGTACCCCAATGTAAAATGTTGTTCCGCCGGGGAAAAAAGCCCCTGTAACCGcccc < 2:566160/90‑1 (MQ=255) tCTGCCAGAATATACGCTGGAAAATTAACCGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAAc < 1:115143/90‑1 (MQ=255) ccAAAATATCCGGGGGAAAATTAAACGGGTATCACCATGTAAAACGTTTTTCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACagag < 2:99115/90‑1 (MQ=255) ccAAAATATACGCGGGAAAATAAACCGGTTACCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag < 2:573163/90‑1 (MQ=255) ataCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgata > 1:35819/1‑90 (MQ=255) gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa > 2:252995/1‑90 (MQ=255) aCGGGTACCACCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc < 1:449178/90‑1 (MQ=255) gggTACCCCCTTGTAAAACTTTTTTACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg < 2:2378/90‑1 (MQ=255) tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc > 1:206249/1‑90 (MQ=255) tgttCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa < 1:186227/90‑1 (MQ=255) cGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc < 1:90574/90‑1 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 1:142203/1‑90 (MQ=255) | GTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT > NZ_CP009273/2621897‑2622062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |