| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,985,528 | 0 | C | T | 66.7% | 24.2 / 19.3 | 24 | D201N (GAC→AAC) | araF | arabinose ABC transporter periplasmic binding protein |
| Reads supporting (aligned to +/- strand): ref base C (3/5); new base T (9/7); total (12/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGT > NC_000913/1985452‑1985603 | ccGACGATCAGCCAATGTTTAACTTCCGGATATTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTa < 2:362665/81‑1 (MQ=255) cgacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAg < 1:362666/81‑1 (MQ=255) acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGAt > 1:362690/1‑81 (MQ=255) aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTggta < 1:362667/81‑1 (MQ=255) cAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTggtaggt < 1:362668/81‑1 (MQ=255) gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCt > 2:362694/1‑81 (MQ=255) aCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATa > 1:362709/1‑80 (MQ=255) ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAAt < 2:362669/79‑1 (MQ=255) ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt > 2:362692/1‑81 (MQ=255) tGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGtttt > 2:362701/1‑81 (MQ=255) tGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtttt > 1:362705/1‑81 (MQ=255) ttGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTc < 2:362670/81‑1 (MQ=255) gAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCg > 2:362706/1‑81 (MQ=255) ccAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGa < 2:362671/81‑1 (MQ=255) aTTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCg > 1:362702/1‑81 (MQ=255) aTTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCg > 1:362698/1‑81 (MQ=255) gTTGGCAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGc < 1:362672/81‑1 (MQ=255) cAGCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa < 1:362673/81‑1 (MQ=255) gCGTCAAATGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc < 1:362674/81‑1 (MQ=255) aTGCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa > 1:362700/1‑81 (MQ=255) gCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa > 1:362707/1‑79 (MQ=255) gCCCCCGGGATGTTGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa > 1:362710/1‑79 (MQ=255) gggATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATccc < 1:362675/81‑1 (MQ=255) gATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgt < 2:362676/81‑1 (MQ=255) | CCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGT > NC_000913/1985452‑1985603 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |