Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 71,575 | 0 | C | T | 56.0% | 3.1 / 26.2 | 25 | K131K (AAG→AAA) | thiB | thiamine ABC transporter substrate binding subunit |
Reads supporting (aligned to +/- strand): ref base C (4/7); new base T (14/0); total (18/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.86e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAA > NZ_CP009273/71451‑71709 | tttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATTAACg > 1:477704/1‑137 (MQ=255) tCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTGTTTTTTTCCTATACCACGGCg > 2:74496/1‑137 (MQ=255) gCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGc < 1:552211/138‑1 (MQ=255) cAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTGTCATAAAAGAAAGCGAAGTAGCCCTa > 1:114498/1‑138 (MQ=255) gCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGGTTGTTTTTGTCATAAACAAAGGAGAAGTAGACATAAACacaaag > 2:329291/1‑134 (MQ=255) cgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTTATAAACCAAGGGGGAGGAGGCATAATCAAAAGGGACGGAAACTGTCttttat > 2:564958/1‑136 (MQ=255) cgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGGAC‑GAAAGTGTCATTATTc > 1:515077/1‑139 (MQ=255) aTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTTTTTTTTTCATAAAAGAAGGCGGAAGAACCATAAACCAAAGGGAC‑GAAAGTGTCATTTTTTCAg > 1:400171/1‑139 (MQ=255) cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTGTCTTAAACTGAGGGGGAGTAGCCACTATCTAACGGTGC‑GTAAGGGGCTTTTTTGCAGGCGCCCGGGACgcttt > 1:263311/1‑135 (MQ=255) cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAGTTTGTTCTTTTCATAAAAACAGGAGAAGAAGCCATAATAAAAAGGTAC‑GAAAGTGTCATTTTTCCAACCGCCGGGGACgctt > 1:379066/1‑135 (MQ=255) ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATc < 1:198022/138‑1 (MQ=255) tcAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTGTTTTTAATAAAGGGGCAGTAGGGAGTATTAAAAATAAA‑GAAAGGGTCATTTATCCAACCGGCGGGGAGGTTAAACGGATCGGCGGGCac‑ac > 2:237489/1‑139 (MQ=255) aCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGc < 1:321898/139‑1 (MQ=255) cTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATAAAAGAAGGGGAAGTAGGCATAATCAAACGGTAC‑GAAAGAGGCATTATTTCAGCCGCCGGGGACGTTAAAGGCATCCGCTGCCAC‑ACCGC‑tttttg > 2:321898/1‑136 (MQ=255) ttCAGGCTTTGTGGCGGGGTTTTTTGGTTGTTTTTGTCATAAACAAAGGGGAAGGAGCCATAATCAAACGGTAA‑CAAAATGTCATTATTTCAGCCCCCCGGGACGTTAACAGCGTCCGCTGCCCC‑ACCCC‑TTTTTGCa > 1:347850/1‑138 (MQ=255) ttCAGGCTTTGTGGCGGGGTTTTTAGGTTGTTTTTGTCATAAAAGAAGGCGAAGGAGCCATAATTAAAAGGTAC‑GAAAGGGTCATTATTTCAGCCGCCGGGGACGGTAACGGCATCCGCTGCCACAACCCC‑TTTTgca > 1:37095/1‑136 (MQ=255) tCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGCCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCa < 2:263311/137‑1 (MQ=255) tCAGGCTTTGTGGCGGGGTTTTTAGGTTTTTCTTTTTGTAAACCAAGGAGGAAGAGGCCTAATAAAAAGGGAG‑GGAAGAGTCCTTATTTCAAGCGCCGGGGAAGGTAACGGGAACCGCCGCCCC‑ACCCCGTTTTggaa > 1:185568/1‑136 (MQ=255) tgtgGCGGGTTTTTTAGTTTTTTTTTTTTATAAAAGAAAGCGGAGTAGCCATAATTAAAAGGTAC‑GGAAGGGTGATTTTTTTAGCCCCCGGGGACGTTAACGGCATCCGCCGCCAC‑ACCCC‑TTTTGGGAAACAGTGCgg > 2:31325/1‑139 (MQ=255) tgtgGCGGGTTTTTTAGTTTTTTTTTTTCATAAACGAAAGGGGAATAGCCCTAATCAAAAGGGAC‑GAAAGTGTCATTATTTCAGCCGCCGGGAACGTTAACGGCAACCGCCGCCCC‑ACCCC‑CTTTTGCAAAAAGTcccg > 1:488916/1‑137 (MQ=255) ttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGGTTTTACTggcg < 2:515077/138‑1 (MQ=255) cAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGGTTTTACTggcggc < 1:442705/138‑1 (MQ=255) tGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTaa > 1:369225/1‑139 (MQ=255) tCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTc < 2:111965/114‑1 (MQ=255) tCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTc > 1:111965/1‑114 (MQ=255) | TTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAA > NZ_CP009273/71451‑71709 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |