Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,621,984 | 0 | A | C | 57.1% | 7.4 / 20.6 | 21 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (7/2); new base C (0/12); total (7/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.10e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT > NZ_CP009273/2621905‑2622062 | aaccacaTTCTGCAAAAATATCCGCGGGAAAATTAAACGGG‑TACCCCCATGTGAACCGTTGTTACGCCGGGGAGTAAAGACCCTGTAAcc < 1:407870/87‑1 (MQ=255) aaaCACACTCTCCAAGAATATCCGCTGGAAAATAAACCGGT‑TCCCACCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAAcc < 2:334497/90‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGACCCTGTAAcc > 1:485407/1‑90 (MQ=255) aaaCAAAATCTGCAAAAATATACGCTGGGAAATTAACCGGG‑TACCCCCATGTAAACCTTTGTTCCGCCGGGGAGTAAAGCCCCTATAAcc < 2:451594/90‑1 (MQ=255) tctcCCAGAATATACGCGGGAAAATTAAACGGG‑TACCCCCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAAc < 2:738353/90‑1 (MQ=255) tccgCCAAAATTTACGCGGGAAAATAAAACGGT‑TCCCCCCATGAAAAAGGTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAAc < 2:621406/87‑1 (MQ=255) tCTGCCAGAATATACGCTGGAAAATTAAACGGT‑TACCACCATGTAAAACGTTTTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAAc < 1:834015/90‑1 (MQ=255) tCTGCAAAAATATCCGC‑GGGAAAATAAACGGGTTACCCCCATGTAAAAGTTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAAc < 2:902734/90‑1 (MQ=255) ccAAAATATACCCGGGAAAATTAAACGGG‑TATCACCATGTGAAACTTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACagag < 1:181247/90‑1 (MQ=255) aGAATATCCGGTGGAAAATAAAACGGG‑TATCACCTTGTAAAACTTTTTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGaa < 1:662706/90‑1 (MQ=255) ataCGGGGGAAAATAAAACGGG‑TACCCCCATGTAAAACGTTTTTACCCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGaatgata < 2:632754/90‑1 (MQ=255) tGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata > 2:326743/1‑90 (MQ=255) agggg‑TACCCCCATGTAAAACGTTGTTACGCGGGGGAAAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc < 2:485407/88‑1 (MQ=255) ggg‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg > 1:474402/1‑90 (MQ=255) tATCCCCATGTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa < 1:453774/90‑1 (MQ=255) aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAg > 2:532041/1‑90 (MQ=255) aaggttgttACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc < 2:843803/87‑1 (MQ=255) tgttACCCCGGGGAAAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa < 2:725954/90‑1 (MQ=255) gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa > 2:719411/1‑90 (MQ=255) cGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc > 1:486173/1‑90 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 1:680971/1‑90 (MQ=255) | AAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT > NZ_CP009273/2621905‑2622062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |