| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 75,079 | 0 | C | T | 28.6% | 24.5 / 12.5 | 21 | L134L (CTG→CTA) | thiB | thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein |
| Reads supporting (aligned to +/- strand): ref base C (1/14); new base T (6/0); total (7/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.29e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGT > NC_000913/74855‑75280 | cTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGtt > 1:253610/1‑229 (MQ=255)cTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTGTTAGt > 1:896870/1‑228 (MQ=255) cTTTTAAAAACAGCCCGTAGGCTCCGCTCCAGCCTTGGGTGACCGTGCCCGTTTTCTTCGCCAGTTTTGCCCGGGCTGGGGGGGCGTCATCGCCATA‑GCTTTTTTGCATCCATAGCACCAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCt < 2:1458066/228‑1 (MQ=255) aaacctggCCGGAGGTTGCGCGCAAGCTTTTGGTGACCGTGCCCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GCCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATa < 2:253168/223‑1 (MQ=255) gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTT‑TTTAAAAGaa > 1:750770/1‑228 (MQ=255) ggCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGCGTTTTTGGTTTGTTTTTGT‑CCCAAACGCAAGCGaa > 1:97171/1‑228 (MQ=255) cgCCTTTGGGGACCGGGCCCTTTTTTATCCCCATTTCCTACCAGGCTGGTGGGGCGTAATCGCCATA‑GTCTTTTTGCATCCATAGCAACAGACCCAGCCCGGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATc < 2:999837/228‑1 (MQ=255) ttttGGTGCCCGTTCCCTTTTTCTTCCCCAGTTCCTACCAGGTTTGGGGGGCTTCATCGCCATA‑GATTTTTTGCATCCATGGCACCAGACCCAGCCCCGGTGTCCTGGTGCGCGGATCCGAATAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCaaa < 2:750254/228‑1 (MQ=255) ggaacggtgaccttTTCCTCCGCCAGTTCCTGCCAGGCTTGTGGGG‑GTCGCTGCCCTAGGACTTTTTGCACCCAAAGCAACACACCCTGCCCGGGGTTGCTGGGGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGt < 2:610903/222‑1 (MQ=255) ccgtTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTT‑TCTTAACCAAAGCGAAGTAGCCATAATCAAAAGGGACGacaatgg > 1:535871/1‑222 (MQ=255) tCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGCTTTTCATTTTTTTTTTTT‑CATAAACAAAAGAGAAGTAGCCATAATCa > 1:1379872/1‑174 (MQ=255) gcATCGCAATA‑GACTTTGTGCACCCATAGAAACAGCCCCAGCCCGGGTGTCGTGGTGGGGGGTCCTTGAAAATGCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCa < 2:1377500/173‑1 (MQ=255) tttGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTa < 1:630694/165‑1 (MQ=255) aGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTTTTTTTTT‑CATAAACCAAAGCGAAAAAGGCATAATCAAACGGTACGaaaaagtc > 1:1518960/1‑147 (MQ=255) ggtgtctgggTGCGCGACCCATGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGTTTTTTCAGTTTGTGCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcggc < 2:1067882/216‑1 (MQ=255) tGGTGCGCGGTTCTTTATAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACa < 1:525182/225‑1 (MQ=255) cGCCAGTTTTGATCGCTCTCAACCGGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTaa < 1:1279341/197‑1 (MQ=255) cAACCAGTTCTTTGAGGCTTGGTGGGGTGTTTATCAGTTTTTCCTTGT‑CAACACCAAAGGCTAAGTAGCCATAATCAACCGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTATCGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTc < 2:764104/229‑1 (MQ=255) aaCCAGTTCTGTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGgttg < 1:1089183/184‑1 (MQ=255) tgtttcgaGGCTGTGGGGCGGGTTGTTCAGTGTGTT‑TTGTGAATAAACGAAGGCGAAGTAGCCATAATAAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGACCGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGt < 2:1913456/222‑1 (MQ=255) tgtgGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATATTCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcg < 1:1093321/150‑1 (MQ=255) | CTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATA‑GACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGT‑CATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGT > NC_000913/74855‑75280 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |