Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 71,561 | 0 | T | G | 58.3% | 9.6 / 26.1 | 24 | N136T (AAC→ACC) | thiB | thiamine ABC transporter substrate binding subunit |
Reads supporting (aligned to +/- strand): ref base T (4/6); new base G (14/0); total (18/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.56e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGG > NZ_CP009273/71432‑71690 | gCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGtt < 2:569582/139‑1 (MQ=255) tCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGttttt > 1:191229/1‑137 (MQ=255) 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) aGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGt > 1:295028/1‑96 (MQ=255) aGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGt < 2:295028/96‑1 (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) | GCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAC‑GAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCAC‑ACCGC‑TTTTGGCAAACAGTCCGG > NZ_CP009273/71432‑71690 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |