Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 63.6% | 16.0 / 34.4 | 33 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (0/12); new base G (21/0); total (21/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.82e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC > NZ_CP009273/822870‑823037 | gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAgg < 1:1181710/90‑1 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑Ggggggggg > 2:1276872/1‑90 (MQ=255) tCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑Ggggggggggg > 2:1190765/1‑90 (MQ=255) aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑Ggggggggggggc > 2:990677/1‑88 (MQ=255) cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAAc < 2:1304801/90‑1 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGAAg > 1:395047/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGAAg > 2:950661/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGAAg > 2:1080567/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGGAg > 2:351086/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGGAg > 1:1518724/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGAAg > 1:93073/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGGAAAGCCGGGAg > 2:208314/1‑90 (MQ=255) cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGGAGGCGGgtcgtc > 2:1445690/1‑90 (MQ=255) tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGCAACGCCGGAAGTCGGGtcgtcgg > 2:907373/1‑89 (MQ=255) cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaaa < 1:1414410/90‑1 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGGGaaaa > 2:684942/1‑89 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGCACGCCCGGAGGCGGTTCGTCGAGaaaaaa > 1:742667/1‑87 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAg > 2:80628/1‑90 (MQ=255) cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa < 2:431053/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:351086/90‑1 (MQ=255) aCCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCa < 1:684942/90‑1 (MQ=255) aCCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAACAGAATGTCCGGGTCa > 2:384348/1‑90 (MQ=255) cATGCTGTTGATGTGCAG‑CAAAAATCACGGCGGGGGGGGGGGGGGAAACCCGGAAGCTGGTTTTCCGGAAAAAAAAATGTCCGGTTcggt > 2:1219856/1‑87 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAATGTCCGGGTCATGc > 2:1002062/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGGAACGCCGGGAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGc > 2:1007888/1‑90 (MQ=255) aTGCCGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGCAACGCCGGGAGGCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGc > 2:382235/1‑90 (MQ=255) tGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCa > 2:18716/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGTTCATGCATCAg > 2:467697/1‑90 (MQ=255) tGATGTGCAGCCAAAATTCACGGC‑GGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGCATCAGCg > 2:703847/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:990677/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:703847/90‑1 (MQ=255) tgCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAAAGCCGGAAGACGGGTCGTCGGGAAAAAGAAAGTCCGGGTTTTGCATCAGCGAACAg > 2:1286438/1‑90 (MQ=255) aTTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcca < 1:1061471/90‑1 (MQ=255) cggc‑ggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:1527875/90‑1 (MQ=255) cggc‑ggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:678348/90‑1 (MQ=255) cggc‑ggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:478385/90‑1 (MQ=255) | GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC > NZ_CP009273/822870‑823037 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |