Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,621,984 | 0 | A | C | 49.1% | 67.2 / 18.1 | 57 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (23/6); new base C (0/28); total (23/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTC > NZ_CP009273/2621897‑2622063 | gTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc > 2:1283184/1‑90 (MQ=255) ttGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAccc < 1:1968556/90‑1 (MQ=255) tGCAGGAAACCAAATCTGCCAGAATATCCGCTGGGAAATTAAACGGG‑TTCCCCCATGTGAACCGTTGTTACGCCGGGGAGTAAAGACCCt < 1:1301417/90‑1 (MQ=255) gaaaaCACAATCTGCCAGAATATCCGCGGGAAAATTAACCGGG‑TACCCCCATGTGAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTaa < 1:1200/88‑1 (MQ=255) gaaaaCAAAATCTGCCAGAATATACGCGGGAAAATAAAACGGG‑TACCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTaa < 1:1185852/88‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAAcc < 2:71525/90‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAAcc < 2:559905/90‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAAcc > 1:2623439/1‑90 (MQ=255) aaaCACAACCTGCCAGAATATCCGGTGGGAAATTAACCGGG‑TACCCCCATGTAAAACGTTTTTACGCGGGGGAGTAAAGACCCTGTAAcc < 2:2512270/90‑1 (MQ=255) aaaCAAAATCTGCCAAAATATACGCTGGGAAATTAAAGGGG‑TACCCCCTTGTGAACCGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAAcc < 1:1720893/90‑1 (MQ=255) caATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGcccc > 1:369435/1‑90 (MQ=255) tCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCCCCATGTGAAACGTTGATCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAAc < 1:22211/90‑1 (MQ=255) tCTGCCAGAATATACGCTGGAAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGACCCTGTAACCGCCCCAAc < 2:977006/90‑1 (MQ=255) tCTGCCAAAATATACGCTGGAAAATTAACCGGG‑TACCCCCATGTGAAACTTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAAc < 2:2411774/90‑1 (MQ=255) tCTGCAAAAATATCCGGGGGAAAATAAAACGGG‑TACCCCCATGGAAAAGGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAAc < 1:2182573/90‑1 (MQ=255) cTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACa > 1:1415624/1‑90 (MQ=255) tcccAGAATATCCGCTGGAAAATAAAACGGG‑TACCACCATGTGAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACag < 1:1888835/88‑1 (MQ=255) tcccAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTAAACCGTTTTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACag < 2:2378590/88‑1 (MQ=255) tGCCAGAATATCCGCGGGGAAATTAACCGGG‑TACCCCCATGTGAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACag < 1:490924/90‑1 (MQ=255) tGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACag > 1:408416/1‑90 (MQ=255) tGCCAGAATATACGCGGGGAAATAAAACGGG‑TACCACCATGTAAAACGTTTTTCCGCCGGGGAAAAAAGACCCTGTAACCGCCCCAACag < 1:1622448/90‑1 (MQ=255) tGCAAGAATATACGCGGGAAAATAAAACGGG‑TATCCCCATGTAAAAGGTTGTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACag < 1:2396533/90‑1 (MQ=255) ccAGAATATCCGCGGGAAAATAAAACGGT‑TACCCCCATGTAAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag < 1:2859581/90‑1 (MQ=255) caaGATTATCCGCTGGAAAATAAACCGGG‑TATCACATGGTAAAACTTTGTTCCGCGGGGGATAAAAGCCCCTGTAACCGCCCCAACagag < 2:1415624/88‑1 (MQ=255) aGAATATCCGGGGGAAAATTAACCGGG‑TACCCCCATGAAAAAGGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGaa < 1:1761304/90‑1 (MQ=255) aGAATATACGCTGGGAAATAAAACGGG‑TATCCCCATGTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGaa < 1:929223/90‑1 (MQ=255) tGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata > 1:1551687/1‑90 (MQ=255) tGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata > 2:2864567/1‑90 (MQ=255) tGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata > 2:111630/1‑90 (MQ=255) gggAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa > 1:1787565/1‑90 (MQ=255) ggAAATTAAACGGG‑TATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATgataatgataat < 2:351827/90‑1 (MQ=255) ggAAAATAAACGGGTTATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTACCCGCCCCAACAGAGAAtgataatgataa < 1:551956/90‑1 (MQ=255) aaaCGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTa < 2:2788196/90‑1 (MQ=255) aCGGG‑TATCCCCATGTAAAACTTTTTTCCGCCGGGGAAAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc < 1:111630/90‑1 (MQ=255) aCGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc > 2:2163316/1‑90 (MQ=255) ggg‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg > 1:2260985/1‑90 (MQ=255) aTCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTaa > 1:201550/1‑90 (MQ=255) tCACCATGTAAAACTTTTTTACGCCGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAAt < 2:1810853/90‑1 (MQ=255) cATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATa > 2:63811/1‑90 (MQ=255) aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAg > 1:1918508/1‑90 (MQ=255) aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAg > 1:1288965/1‑90 (MQ=255) tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc > 2:778960/1‑90 (MQ=255) tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc > 2:1548408/1‑90 (MQ=255) tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc > 1:2599882/1‑90 (MQ=255) aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt > 1:1627676/1‑90 (MQ=255) aaggtttttACGCGGGGGAAAAAAGCCCCTTAACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc < 1:2609072/87‑1 (MQ=255) tttgttCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa < 2:2599882/89‑1 (MQ=255) tgttACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa < 1:166812/90‑1 (MQ=255) gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa > 2:660969/1‑90 (MQ=255) gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa > 1:1452456/1‑90 (MQ=255) cgcgGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc < 1:396250/90‑1 (MQ=255) gggggAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc < 2:1627676/89‑1 (MQ=255) gggggAAAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc < 2:2623439/89‑1 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATTATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 1:1399968/1‑90 (MQ=255) ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct > 2:570140/1‑90 (MQ=255) gggAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc < 1:778960/90‑1 (MQ=255) gggAGTAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc < 2:1164026/90‑1 (MQ=255) | GTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTC > NZ_CP009273/2621897‑2622063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |