| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_017635 | 76,754 | 0 | C | T | 33.3% | 44.0 / 13.2 | 27 | L134L (CTG→CTA) | thiB | thiamine ABC transporter substrate binding subunit |
| Reads supporting (aligned to +/- strand): ref base C (7/11); new base T (9/0); total (16/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.02e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTTCTGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAACGGG > NC_017635/76705‑76802 | aGTTCTGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTc < 2:1293658/50‑1 (MQ=255) gTTCTGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTcc > 1:317090/1‑49 (MQ=255) ttCTGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGGTTTTTAg > 1:1708026/1‑50 (MQ=38) ttCTGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGCTTGTTctg > 2:1471778/1‑48 (MQ=38) tGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGGTTTTTTGttt > 1:281461/1‑50 (MQ=37) tGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGGTTTTTAGttt > 1:1193123/1‑50 (MQ=38) ctctCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTTAGTtttttttt > 1:1034972/1‑50 (MQ=38) ctctCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCtt < 1:772744/50‑1 (MQ=255) ctctCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCtt < 2:1655255/50‑1 (MQ=255) tctcAACCAATTCTTTCAGGCTTTGCGGCGGGGTTTTTTGTttttttttt > 2:137608/1‑49 (MQ=21) ctcAACCAATTCTTTCAGGCTTTGCGGCGGGGTTTTTTTTTTGTtttttt > 2:1627462/1‑48 (MQ=21) tcAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTg < 1:447902/48‑1 (MQ=255) aCCAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTCATa < 2:57505/50‑1 (MQ=255) cAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTCATAGa < 1:1600297/50‑1 (MQ=255) cAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTCATAGa < 1:1611224/50‑1 (MQ=255) cAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTAATAGa < 2:147658/50‑1 (MQ=255) aaTTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTCATAGAc < 2:1823405/50‑1 (MQ=255) aaTTCTTTCAGGCTTTGCGGCGGGGTTTTTAGTTTGTTTTTGTCATAGAc > 1:1783739/1‑50 (MQ=38) aTTCTTTCAGGCTTTGCGGCGGGGTTTTTCGTTTTTTTTTGTCATAGAc > 1:1893698/1‑49 (MQ=25) aTTCTTTCAGGCTTTGCGGCGGGGTTTTTAGGTTTTTCTTGTCATAGAc > 1:1787898/1‑49 (MQ=25) gcggcgGGTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGcc < 1:1088496/49‑1 (MQ=255) gcggcgGGTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGcc < 2:1903473/49‑1 (MQ=255) cggGTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAAt > 2:1969128/1‑50 (MQ=255) gggTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAAt > 1:1245746/1‑49 (MQ=255) gggTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAAt > 2:702823/1‑49 (MQ=255) tttttCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCaaa > 1:643538/1‑50 (MQ=255) tttCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAAc < 2:1952503/49‑1 (MQ=255) ttCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAActg > 1:459960/1‑48 (MQ=255) ttCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAACgg > 1:925793/1‑50 (MQ=255) tCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAACgg > 1:819576/1‑49 (MQ=255) cAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAACggg < 1:59168/49‑1 (MQ=255) cAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAACggg < 2:486186/49‑1 (MQ=255) | AGTTCTGATCGCTCTCAACCAATTCTTTCAGGCTTTGCGGCGGGTTTTTCAGTTTGTTCTTGTCATAGACGAAGGCGAAATAGCCGTAATCAAACGGG > NC_017635/76705‑76802 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |