Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,740,492 | 0 | G | A | 39.5% | 41.2 / 13.4 | 39 | A17T (GCC→ACC) | lyxK | L‑xylulokinase |
Reads supporting (aligned to +/- strand): ref base G (4/19); new base A (15/0); total (20/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.01e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAG > NZ_CP009273/3740404‑3740576 | aGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGc < 1:232454/90‑1 (MQ=255) gCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGcc < 2:134209/90‑1 (MQ=255) gCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTCAGATTGTGGCGGTAGCT‑GGCTGAAAGcc < 1:297041/90‑1 (MQ=255) tGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCgg < 2:595449/90‑1 (MQ=255) tGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCgg < 1:279029/90‑1 (MQ=255) gATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGGAGCTGGGTTGAAACCCgg > 1:492491/1‑90 (MQ=255) ctgctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCTg > 1:192823/1‑90 (MQ=255) gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTa < 2:110687/90‑1 (MQ=255) gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCTGTa > 1:101728/1‑90 (MQ=255) gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCTGTa > 1:13923/1‑90 (MQ=255) gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCTGTa > 1:28406/1‑90 (MQ=255) aaGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGc < 2:133693/81‑1 (MQ=255) aaGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGc > 1:133693/1‑81 (MQ=255) aGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGGTGAAAACCGGGCGGTATGACCGCGGAGGGCggg > 2:336577/1‑90 (MQ=255) aGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGGAGCT‑GGCTGAAAACCGGGGTGTATGACCCCGAAGGGCggg > 2:170330/1‑90 (MQ=255) gCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGaggc < 2:13923/90‑1 (MQ=255) aTAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGaggcagg < 1:370999/90‑1 (MQ=255) aTGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGGAGGCGGGCAg > 1:48546/1‑90 (MQ=255) aTGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAg < 2:617095/90‑1 (MQ=255) ggctggGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCTGTATGACCGCGAAGGCCGGGGGGGAGGCGGGCAGCGCCCGCCGCtgt > 2:289435/1‑90 (MQ=255) ggctggGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCCGTATTACCCCCAAGGGCGGGGGGGAGGCGGGCCGCGCCCGCCGCtgt > 2:16734/1‑90 (MQ=255) ctggGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTgc < 1:73636/90‑1 (MQ=255) ctggGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTgc < 2:212705/90‑1 (MQ=255) gggTTAGATTGTGGCGGTAGCT‑GGCTGAAAACCGGGCTGTATGACCGCGAAGGCCGGGGGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCa > 1:656145/1‑90 (MQ=255) gggTTAGATTGTGGCGGGAGCT‑GGCTGAAAACCGGGGTGGATGACCGCGAAGGGCGGGAGGGAGGGGTGGAGCGCCCGCCGCCGGGCGCa > 2:73636/1‑90 (MQ=255) gTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCAtt < 1:385178/90‑1 (MQ=255) ttGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGccc < 1:289435/90‑1 (MQ=255) gtgGCGGTAGCT‑GGCTGAAAACCGGGGTGTATGACCGCGAAGGCCGGGGGGGAGGCGGGCGGCGCCTGCCGCTGTGCGCATTAAACCCGc > 2:385178/1‑90 (MQ=255) gtgGCGGTAGCT‑GGCTGAAAACCGGGCTGTATGACCGCGAAGGCCGGGGGGCAGGGGGGCAGCGCCCGCCGCTGTGCGCATTAAGCCCGc > 2:404751/1‑90 (MQ=255) gtgGCGGTAGCT‑GGCTGAAAACCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGGGGGCAGCGCCTGCCGCTGTGCGCATTAAGccccc > 2:95814/1‑88 (MQ=255) gtgGCGGTAGCT‑AGCTGAAAGCCGGGGTGTATGACCGCGAAGGCCGGGGGGCAGGCGGGGAGAGCCTGCCGCTGGGCGCGTTAAGCCCGc > 2:517063/1‑90 (MQ=255) gCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGc < 2:101728/90‑1 (MQ=255) gTAGCT‑GGCTGAAAACCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAg > 1:289722/1‑90 (MQ=255) ct‑ggctgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTg < 2:317029/90‑1 (MQ=255) ggctgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTggg > 2:613735/1‑90 (MQ=255) gctgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGc < 1:465836/90‑1 (MQ=255) tgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCag < 2:506140/90‑1 (MQ=255) aaGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAgcg > 2:126676/1‑46 (MQ=255) aaGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAgcg < 1:126676/46‑1 (MQ=255) | AGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCT‑GGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGTGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAG > NZ_CP009273/3740404‑3740576 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |