Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,773,700 | 0 | A | T | 31.0% | 52.4 / 17.6 | 29 | N474Y (AAC→TAC) | mtlA | mannitol‑specific PTS enzyme IICBA component |
Reads supporting (aligned to +/- strand): ref base A (10/10); new base T (0/9); total (10/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.14e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTT > NC_000913/3773644‑3773761 | gCGGCCTGTAAACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACCCCACACGGCAAACGAagag > 1:207749/1‑65 (MQ=255) cGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga < 1:2854687/65‑1 (MQ=255) cTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGt < 1:529309/65‑1 (MQ=255) acacCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaaa > 1:3694021/1‑65 (MQ=255) acCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa > 1:2439287/1‑65 (MQ=255) ccAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc > 1:3861647/1‑65 (MQ=255) ccAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc > 1:2558480/1‑65 (MQ=255) aGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg > 1:1872259/1‑65 (MQ=255) gCCTGACCGAACGTCTGGTTGCTTCCCAGCGCAACACGGCATACGAAGAGAAAGTAAAAGACAGc < 1:694660/65‑1 (MQ=255) gCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc < 1:1594657/65‑1 (MQ=255) gACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGa < 1:142215/65‑1 (MQ=255) gACCGAACGTCTGGTTGCTGCCCAACGCAACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGa < 1:3517685/65‑1 (MQ=255) aCCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaa < 1:3208286/65‑1 (MQ=255) gTCTGGTTGCTGCCCAACAGCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:1492161/65‑1 (MQ=255) gTCTCGTTGCGGCCCAACGCCACACGGCAGTCGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:1231957/65‑1 (MQ=255) gctgcCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGa > 1:3142228/1‑65 (MQ=255) cccAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTcc > 1:2324942/1‑65 (MQ=255) ccAACGCCACACGGCAAACGAAGAGAAAGCAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCa > 1:3011843/1‑65 (MQ=255) taCGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:3493507/64‑1 (MQ=255) aaCGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:2404525/65‑1 (MQ=255) aaCGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:502167/65‑1 (MQ=255) aaCGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:2278167/65‑1 (MQ=255) aaCGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:1273391/65‑1 (MQ=255) aaCACCACACGGCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:627416/65‑1 (MQ=255) acGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCt < 1:968718/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:3146914/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:2994277/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:2974087/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:2788001/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:244197/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:225292/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:655052/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:1461865/65‑1 (MQ=255) ggCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:2860440/65‑1 (MQ=255) gCAAACGAAGAGAAAGTAAAAGATAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt > 1:755952/1‑65 (MQ=255) gCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt > 1:177434/1‑65 (MQ=255) | GCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTT > NC_000913/3773644‑3773761 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |