Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 645,920 | 0 | A | T | 46.4% | 20.8 / 13.0 | 28 | R9R (CGT→CGA) | citE | citrate (pro‑3S)‑lyase subunit beta |
Reads supporting (aligned to +/- strand): ref base A (2/13); new base T (13/0); total (15/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.83e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGG > NZ_CP009273/645837‑645998 | tCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGCGGTTTTACgtgt > 1:1337314/1‑88 (MQ=255) tCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGCGTG‑GGGTTTTCgtttt > 1:1261165/1‑88 (MQ=255) cGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTCGTTTTTGCAg > 1:1858865/1‑90 (MQ=255) aGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTTGTTGTTTGAGCGGAGAGGgaa > 1:373774/1‑90 (MQ=255) ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTTGTTGTTGGAGCGAAGGGGaaaa > 1:613185/1‑89 (MQ=255) ttGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTTGTTTTTGGAGCGGAGCGGGAAACAtt > 2:1065854/1‑90 (MQ=255) ttGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCG‑GGGTTTTTGTTTTTGGAGCGAAGAGGGAAACAtt > 2:1231967/1‑90 (MQ=255) gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTCGTTGTTTCAGCGAAGCGGAAAACATTTGcga > 1:226673/1‑88 (MQ=255) ccATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAtc < 2:373774/90‑1 (MQ=255) ccATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAtc < 1:854918/90‑1 (MQ=255) ccATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAtc < 1:962831/90‑1 (MQ=255) atctcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCG‑GGGTTTTCGTTGTTGGAGGGAAGCGGGAATAATTTGGAAtcct > 2:308505/5‑90 (MQ=255) gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTccc < 1:1663182/90‑1 (MQ=255) ggCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCAt < 1:1258472/90‑1 (MQ=255) acagggCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGcc < 2:487769/86‑1 (MQ=255) aCCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGcc < 1:1231966/90‑1 (MQ=255) cAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTTGTTGTTGCAGCGAAGAGGGAAACATTGGCAATCCTCCCCTGGCAGAGGCgg > 1:1328106/1‑90 (MQ=255) cAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTCGTTTTTTGAGCGGAACGGAAAAAAATGGGAATACTCCCCTCGCAGAGcccg > 1:922801/1‑88 (MQ=255) cAGGCACAAACAACATGCTGCGGCGGGGGCG‑GGGTTTTTGTTTTTGGAGCGAAGCGGGAAACATTTGCAATCCTCCCCCGGCAGGGGCgg > 1:1797179/1‑90 (MQ=255) gCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGAt < 1:1065853/90‑1 (MQ=255) gCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGAt < 1:1773237/90‑1 (MQ=255) aacaacATGCTGCGGCGGGGGCG‑GGGTTTTTGTTGTTTCAGCGAAGCGGGAAACATTGGCAAACCTCCCCCGGCCGAGGCCGGGGACCGc > 2:1726258/1‑90 (MQ=255) acaTGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgc < 2:961890/90‑1 (MQ=255) acaTGCTGCGGCGGGGGCG‑GGGTTTTTGTTGTTGCAGCGAAGCGGGAAACATTGGCAATCCTCCCCTGGCAGAGCCGGGGGACCGCTggg > 2:1902072/1‑89 (MQ=255) aTGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcgc < 2:1919709/90‑1 (MQ=255) cggcggGGGGG‑GGGTTTACGTTGTTTCAGCGAAGCGGAAAACAATGGCAAACCTCCCCTGGCAGGGCCGGGGTACCCCCGGGGCGGGccc > 2:1622076/1‑89 (MQ=255) gggTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagcag < 2:399372/90‑1 (MQ=255) gTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAggg < 2:1858866/90‑1 (MQ=255) | TCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGG > NZ_CP009273/645837‑645998 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |