Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,163,417 | 0 | A | T | 74.2% | 58.2 / 18.0 | 31 | D53V (GAT→GTT) | thiK | thiamine kinase |
Reads supporting (aligned to +/- strand): ref base A (5/3); new base T (9/14); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.12e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGC > NC_000913/1163339‑1163497 | cgTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTc < 1:214524/81‑1 (MQ=255) gTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTTAGCCGCACGATcc > 1:214559/1‑81 (MQ=255) gTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATgcg < 2:214525/80‑1 (MQ=255) ggCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGcgccgc < 1:214526/81‑1 (MQ=255) gCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCa < 2:214527/81‑1 (MQ=255) cTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTTCGTCAGCCGCACGTTCCTGATGCGCCGCAg > 2:214557/1‑81 (MQ=255) gTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCcg > 2:214564/1‑81 (MQ=255) ggCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCcgcg > 1:214566/1‑81 (MQ=255) ggAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCt < 1:214528/81‑1 (MQ=255) aGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTg > 1:214558/1‑81 (MQ=255) ttttCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTgcg > 1:214556/1‑81 (MQ=255) tctcATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTgcgc > 1:214577/1‑79 (MQ=255) tctcATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGcc < 1:214529/80‑1 (MQ=255) aTTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAg < 1:214532/78‑1 (MQ=255) aTTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGt < 1:214531/79‑1 (MQ=255) aTTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTa < 2:214530/80‑1 (MQ=255) aTTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTAt > 2:214570/1‑81 (MQ=255) aTTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTAt > 1:214569/1‑81 (MQ=255) aTTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTa < 2:214533/80‑1 (MQ=255) tGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCg > 1:214568/1‑81 (MQ=255) cATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCt > 1:214573/1‑81 (MQ=255) gCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATcac < 2:214534/81‑1 (MQ=255) ttgtGCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCgcat < 1:214535/81‑1 (MQ=255) tgCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCgcatgca > 2:214567/1‑81 (MQ=255) tgCGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCgcatgc < 2:214536/80‑1 (MQ=255) cGTCAGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCAtt < 1:214537/81‑1 (MQ=255) aGCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCAc < 2:214538/81‑1 (MQ=255) aGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCAc > 2:214581/1‑81 (MQ=255) gCCGCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCAcc < 2:214539/81‑1 (MQ=255) gCACGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGaa < 1:214540/81‑1 (MQ=255) aCGTTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGc > 2:214578/1‑81 (MQ=255) aTCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGc < 1:214542/81‑1 (MQ=255) | CGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGC > NC_000913/1163339‑1163497 |
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: ‑ |