Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,621,984 | 0 | A | C | 55.0% | 18.3 / 27.6 | 20 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/11); total (9/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTC > NZ_CP009273/2621905‑2622063 | aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAAcc > 2:207480/1‑90 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATAAAACGGGTATAACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGCCCCTGTAAcc < 1:150932/90‑1 (MQ=255) acaATCTGCCAAAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGccc < 2:692720/90‑1 (MQ=255) tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCGTGTAAAACG‑TTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAAc < 2:196457/90‑1 (MQ=255) tGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACag > 2:852573/1‑90 (MQ=255) ccAAAATACCCGCGGGAAAATTAACGGGGTACCCCAGTGTAAAACG‑TTGTTACGCCGGGGAGCAACGCCCCTGTAACCGCCCCAACagag < 1:460034/90‑1 (MQ=255) cAGAATATGCGGTGGGAAATTAAAGGGGTAACCCCAGGAGAAACG‑TTTTTACGCGGGGGGGTAAAGCCCTTGTAACCGCCCCAACAgaga < 2:467070/90‑1 (MQ=255) aGAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGaa < 2:179672/90‑1 (MQ=255) cgcgGGAAAATAAAACGGGTATCCCCATGTAAAACG‑TTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGaatgataatg < 1:207480/90‑1 (MQ=255) tGGGAAATTAAACGGGTATCACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata > 1:257172/1‑90 (MQ=255) tGGGAAATTAAACGGGTATCACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata > 2:750301/1‑90 (MQ=255) ggTATCACCATGTAAAACG‑TTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt < 2:289420/90‑1 (MQ=255) aTGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAg > 1:892670/1‑90 (MQ=255) gtgAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg > 1:704210/1‑90 (MQ=255) gtgAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg > 1:447964/1‑90 (MQ=255) aaaCGTTTGTTCCGCGGGGGAGAAAAGCCCTTATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt < 1:852573/90‑1 (MQ=255) tgttACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa < 2:704210/90‑1 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 2:9662/1‑90 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 1:158742/1‑90 (MQ=255) gggAGTAAAGCCCTTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc < 2:412982/90‑1 (MQ=255) | AAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTC > NZ_CP009273/2621905‑2622063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |