Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,963 | 0 | T | G | 60.0% | 7.7 / 41.4 | 35 | D492A (GAC→GCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (3/11); new base G (21/0); total (24/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.72e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-07 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAA > NZ_CP009273/822878‑823051 | cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAAc < 2:1304801/90‑1 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGAAg > 2:1080567/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGAAg > 2:950661/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGGAACGCCGGAAg > 1:395047/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGGAg > 1:1518724/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGGGGGGGGGGCAACGCCGGGAg > 2:351086/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:678348/90‑1 (MQ=255) cggc‑ggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:478385/90‑1 (MQ=255) cggc‑ggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:1527875/90‑1 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa > 2:566664/1‑90 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa > 1:1071940/1‑90 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa > 1:1018972/1‑90 (MQ=255) | CATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAA > NZ_CP009273/822878‑823051 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |