Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 645,914 | 0 | A | G | 51.7% | 28.8 / 37.4 | 29 | T11T (ACT→ACC) | citE | citrate (pro‑3S)‑lyase subunit beta |
Reads supporting (aligned to +/- strand): ref base A (0/14); new base G (15/0); total (15/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGG > NZ_CP009273/645825‑645998 | aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑a < 2:325744/90‑1 (MQ=255) cATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGt < 2:1648309/90‑1 (MQ=255) 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:1773237/90‑1 (MQ=255) gCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGAt < 1:1065853/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) | AACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCG‑AGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGG > NZ_CP009273/645825‑645998 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |