Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,628,625 | 0 | A | C | 28.6% | 51.3 / 10.5 | 28 | V105G (GTC→GGC) | pdeF | cyclic di‑GMP phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (20/0); new base C (0/8); total (20/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA > NC_000913/2628486‑2628764 | gCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGa < 2:291564/140‑1 (MQ=255) tAATCCTACAAATGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGCGGGGAAATAAACCGGGTACCCCCATGTAAAAC‑GTTTTTACGCGGGGGAGAAAAGCCCCTGTAAc < 2:51744/140‑1 (MQ=255) aaacccTAAAAATGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAAAATATCCGCGGGAAAATTAACCGGGTACCCCCATGTGAACC‑TTTGTTACGCCGGGGAGTAAAGCCCCTGAAAc < 2:46483/136‑1 (MQ=255) aaaaTGCAGAAAACTGAAAAAGTATCAAAAACAGCGTTGAAGAAAACACAATCTGCCAGAATATCCGCGGGAAAATAAAACGGGTACCCCCATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCaa < 2:1038705/139‑1 (MQ=255) aaaaTGAAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCAAGAATATACGCGGGAAAATTAACCGGGTACCCCCATGTAAAAC‑GTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCaa < 2:471871/139‑1 (MQ=255) aaggcagaaAACGGAAAAAGTATAAGAAACAGCTTGGCAGAAAACAAAATCTGCAAGAATATCCGCGGGAAAATTAAAGGGGTATCCCCATGTAAAAC‑GTTGTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACa < 2:494704/137‑1 (MQ=255) aaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACa > 1:219300/1‑140 (MQ=255) aaTGCAGAAAACTGGAAAAGTATCAGAAACGGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCGGGGAAATAAAGCGGGTCCCCCCATGTAAAAC‑GTTTTTCCGCGGGGGAGAAAAGCCCCTGAACCCGCCCCAACa < 2:30157/140‑1 (MQ=255) aaTGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGCAGGAAACACAATCTGCAAAATTATACGTGGGTAAATAAAAGGGCTTTCACCAGGTGAATT‑GTTGTTCCCCGGGGGAGTAAACCCCCTGTAACCGCCCCAACa < 2:677973/140‑1 (MQ=255) tGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACaga > 2:647419/1‑140 (MQ=255) gcagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACagag > 2:857480/1‑140 (MQ=255) cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga > 1:1088310/1‑140 (MQ=255) aagcaAACGGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGCTGGAAAATTAAACGGGTATCACCATGTGAAAC‑TTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAgaga < 2:867412/139‑1 (MQ=255) agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaa > 2:862346/1‑140 (MQ=255) agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaa > 1:1053193/1‑140 (MQ=255) aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt > 2:474273/1‑140 (MQ=255) aaGCGTGGCAGAAAACACAATCTGCAAGAATATCCGCGGGAAAATTAACCGGGTATCACCATGTAAAAC‑GTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg < 2:1088310/139‑1 (MQ=255) aGCGTTGCAGAAAACACAATCTGCAAAAATATACGCGGGAAAATAAAACGGGTACCCCCATGTAAAAC‑TTTTTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt < 2:239809/140‑1 (MQ=255) gCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGa > 1:763803/1‑140 (MQ=255) tGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa > 2:881152/1‑140 (MQ=255) gCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaaa > 1:770370/1‑140 (MQ=255) aaTATACGCGGGAAAATTAAACGGGTATCACCATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg < 1:857480/140‑1 (MQ=255) tatcCGCTGGAAAATTAAGCGGTTATCCCCATGTGAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct < 2:219300/140‑1 (MQ=255) gggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat > 2:724309/1‑140 (MQ=255) ggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt > 2:639301/1‑140 (MQ=255) ggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt > 2:1138722/1‑140 (MQ=255) aaTTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc > 1:1037059/1‑140 (MQ=255) ggTATCACCATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaaa < 2:763803/140‑1 (MQ=255) ccccATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt < 1:324906/138‑1 (MQ=255) ttggAAACGGTTGTTACGCGGGGGAATAAAGCCCCTGTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg < 2:427271/137‑1 (MQ=255) gtgAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg < 2:160430/140‑1 (MQ=255) tgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa > 2:1148787/1‑140 (MQ=255) gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg > 1:160430/1‑140 (MQ=255) gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg > 2:241610/1‑140 (MQ=255) gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCa > 1:675598/1‑140 (MQ=255) ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa > 1:928702/1‑140 (MQ=255) ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa > 1:1090233/1‑140 (MQ=255) aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa > 1:497890/1‑140 (MQ=255) aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAAGAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa > 2:964534/1‑140 (MQ=255) | GCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA > NC_000913/2628486‑2628764 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |