Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,621,984 | 0 | A | C | 58.3% | 3.8 / 10.6 | 12 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/7); total (5/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.26e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTT > NZ_CP009273/2621905‑2622072 | aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAAcc > 1:647207/1‑90 (MQ=255) ccAAAATATCCGCGGGAAAATAAAACGGGTATCCCCATGTAAAACTTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag < 1:139502/90‑1 (MQ=255) aCGCTGGAAAATAAAACGGGTACCACCATGAAAAACTTTTTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGaatgataat < 2:35853/90‑1 (MQ=255) aCGGTTACCCCCATGTAAAACTTTTTTACGCGGGGGAAAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc < 2:309055/90‑1 (MQ=255) aCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc > 2:35219/1‑90 (MQ=255) tATCCCCATGTAAAACTTTTTTACCCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa < 1:410047/90‑1 (MQ=255) aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAg > 2:618681/1‑90 (MQ=255) gttgttACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAATAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa < 1:183753/90‑1 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 2:380014/1‑90 (MQ=255) gggAGAAAAGCCCCTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc < 2:59699/90‑1 (MQ=255) gAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg < 1:380014/90‑1 (MQ=255) gACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt > 1:59699/1‑90 (MQ=255) | AAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTT > NZ_CP009273/2621905‑2622072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |