| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 71,561 | 0 | T | G | 52.0% | 31.6 / 39.2 | 25 | N136T (AAC→ACC) | thiB | thiamine ABC transporter substrate binding subunit |
| Reads supporting (aligned to +/- strand): ref base T (0/12); new base G (13/0); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.92e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAA > NZ_CP009273/71485‑71645 | aCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGgttttttgtttttt > 1:18300/1‑90 (MQ=255) ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTtttttttt > 2:593017/1‑90 (MQ=255) ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTtttttttt > 1:427441/1‑90 (MQ=255) ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTtttttttt > 1:966084/1‑90 (MQ=255) ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGGTCtt > 2:449058/1‑90 (MQ=255) gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTttttttttttt > 2:745375/1‑87 (MQ=255) gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTttttttttttt > 2:247943/1‑87 (MQ=255) gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTGtt > 1:87382/1‑89 (MQ=255) gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTttttttt > 2:552404/1‑87 (MQ=255) ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGcc < 1:659187/90‑1 (MQ=255) tttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTAAAAAAAAGGGGAAG‑AGCCATAATCaa > 1:24965/1‑90 (MQ=255) tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAAc < 2:966084/90‑1 (MQ=255) tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAAc < 1:97721/90‑1 (MQ=255) tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAAc < 1:552404/90‑1 (MQ=255) tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAAc < 1:247943/90‑1 (MQ=255) gATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACg < 2:356563/90‑1 (MQ=255) gATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTTTTTTTGTCAT‑AAAAGAAGGCGAAGTAGCCATAATaaaaag > 2:36341/1‑88 (MQ=255) aTCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACgg < 2:389288/90‑1 (MQ=255) ctctCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTTTTTGTTTTTTTTAT‑TAAAAAAGGGGAAGGAGGGATAATTAAAAAGTACg > 1:464000/1‑90 (MQ=255) aaCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgt < 2:778031/90‑1 (MQ=255) gTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCat < 2:13410/90‑1 (MQ=255) gTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTTTCAT‑AAAAGAAGGGGAAGGAGCCATAATCAAAAGGGACGAAAGTGTCat > 2:361069/1‑90 (MQ=255) ttCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTTTCAT‑AAAAAAAGGGGAAGGAGCCATAATCAAAAGGGACGAAAAAGTCattattt > 2:456427/1‑89 (MQ=255) aGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAg < 1:158353/90‑1 (MQ=255) ggCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGa < 1:449058/90‑1 (MQ=255) gCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGaa < 1:257065/90‑1 (MQ=255) | ACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAA > NZ_CP009273/71485‑71645 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |