Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,773,700 | 0 | A | T | 32.3% | 53.8 / 18.2 | 31 | N474Y (AAC→TAC) | mtlA | mannitol‑specific PTS enzyme IICBA component |
Reads supporting (aligned to +/- strand): ref base A (11/10); new base T (0/10); total (11/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.90e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGT > NC_000913/3773648‑3773760 | ccTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAg > 1:2312209/1‑65 (MQ=255) cTGTACACCAGCCTGACCGAGCGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGt < 1:1881240/65‑1 (MQ=255) acacCAGCCTGACCGAACGTCTGGCTGCTGCCGAACGCCACACGGCAAACGAAGAGAAAGTaaaa > 1:574291/1‑65 (MQ=255) acCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa > 1:3186227/1‑65 (MQ=255) ccAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc > 1:2301226/1‑65 (MQ=255) ccAGCCTGACCGAACGTCTGGTTGCTGCCCAAAGCCACACGGCAAACGAAGAGAAAGTAAAAGaa > 1:1864670/1‑64 (MQ=255) ccTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGcc > 1:389940/1‑65 (MQ=255) cTGACCGAACGTCTGGTTGTTGACCAACGAAACACGGCAAACGAAGAGAAAGTAAAAGACAGCCt < 1:1740510/65‑1 (MQ=255) cTGAACGAACGTCTGGTTGTTTCCAAACGCAACACGGCATACGAAGAGAAAGTAAAAGACAGCCt < 1:2446195/65‑1 (MQ=255) aCCGAACGTCTGGTTGCTGCCAAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaa < 1:3021771/65‑1 (MQ=255) ccGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaaa > 1:2617347/1‑65 (MQ=255) ccGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaaa > 1:1980890/1‑65 (MQ=255) ccGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGACAAAGTAAAAGACAGCCTGaaa > 1:938295/1‑65 (MQ=255) cGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAg > 1:2393705/1‑65 (MQ=255) gTCTGGTTGCTTCCCAACGCAACACGGCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:238738/65‑1 (MQ=255) gTCTGGTTGCTTCCAAGCAGAACACGGCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:238363/65‑1 (MQ=255) gTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:1376113/65‑1 (MQ=255) gTCTGGTTGCTGCCCAACACCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:1713314/65‑1 (MQ=255) gTCTCGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc < 1:1897796/65‑1 (MQ=255) tgctgcCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACg < 1:181554/65‑1 (MQ=255) tgcCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGAtt < 1:1511676/65‑1 (MQ=255) cccAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTcc > 1:2547449/1‑65 (MQ=255) aaCGCCACGCGGCATACGAAGAGAAAGTAAAAGCCAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:1363257/65‑1 (MQ=255) aaCGCCACACGGCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:1228707/65‑1 (MQ=255) aaCGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:2026444/65‑1 (MQ=255) aaCGACACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:2659565/65‑1 (MQ=255) aaCACCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt < 1:2259799/65‑1 (MQ=255) cGCCACAGGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGc < 1:237582/65‑1 (MQ=255) cacGGCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc < 1:446457/65‑1 (MQ=255) acGGCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCt < 1:21124/65‑1 (MQ=255) acGGCAGACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCt < 1:401245/65‑1 (MQ=255) acGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCt > 1:1386724/1‑65 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:1294126/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:1104915/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:3198065/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:1389925/65‑1 (MQ=255) ggCATACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:1769401/65‑1 (MQ=255) ggCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTgc < 1:872525/65‑2 (MQ=255) ggCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:2951098/65‑1 (MQ=255) ggCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:2806958/65‑1 (MQ=255) ggCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:1560593/65‑1 (MQ=255) ggCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt < 1:706602/65‑1 (MQ=255) | CCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGT > NC_000913/3773648‑3773760 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |