| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 75,079 | 0 | C | T | 46.5% | 29.4 / 40.9 | 43 | L134L (CTG→CTA) | thiB | thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein |
| Reads supporting (aligned to +/- strand): ref base C (4/19); new base T (20/0); total (24/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.65e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACA‑CCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAC > NC_000913/74935‑75223 | ggCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGtt < 2:851093/149‑1 (MQ=255) gCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTTGttttt > 1:838780/1‑149 (MQ=255) gtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑ATCTTccatttttt > 2:728348/1‑140 (MQ=255) ggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTATTTTGTttttttt > 1:270798/1‑147 (MQ=255) gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATa < 1:908520/150‑1 (MQ=255) tCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATaa < 1:463122/150‑1 (MQ=255) cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTATTTTGTTTTTGTTAt > 2:50716/1‑147 (MQ=255) cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTAGTttttttttttttttaa > 2:154699/1‑142 (MQ=255) cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑a < 2:496059/150‑1 (MQ=255) ccATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAgg < 1:231498/151‑1 (MQ=255) cATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTCGGTTTTTTTTTTTCTTAAACa‑aa > 2:492213/1‑149 (MQ=255) ttgACTTTTTGCATCCATAGCACCAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCg < 2:526740/148‑1 (MQ=255) tAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTTAGTTTGTTCTTGTTATATAAG‑AAGGCgg > 2:558613/1‑150 (MQ=255) tAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTAGGTTGTTTTTTTTATTAACG‑AAGGCg > 2:320825/1‑150 (MQ=255) gACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAg < 1:245100/151‑1 (MQ=255) tttttGCACCCATAGCAACAGACCCGGCCCCGGTGTACTGGTGCGGGGATCCGTATAAATCTCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTa < 2:826396/150‑1 (MQ=255) ttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTa < 1:788628/149‑1 (MQ=255) tttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTAGTTTGTTCTTGTTATAAAAG‑AAGGCGAAGGA‑g > 2:848096/1‑149 (MQ=255) aTCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTTTGTTTTTTTTTTTCATAAAAA‑AAGGCGGAGTA‑GCCATaa > 1:328550/1‑150 (MQ=255) ccATAGCAACAGCCCCAGCCCCGGTGTCCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAAt < 2:6968/149‑1 (MQ=255) ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAAt < 1:417810/149‑1 (MQ=255) aTAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCaa < 1:400930/150‑1 (MQ=255) aGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAAc < 1:163980/150‑1 (MQ=255) gCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTAGTTTTTTTTTTTTCTAAAAG‑AAAGGGGAGGA‑GCCATAATCaaa > 1:810197/1‑148 (MQ=255) gCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTAGTTTTTTTTTTTCCTAAACA‑AAAGGGGAGGA‑GCCACAATCaaa > 1:616262/1‑148 (MQ=255) cAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACgg < 2:166759/150‑1 (MQ=255) aaCAGACCCAGCCCCGGTTTACTGGTGCGCGGTTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAAc < 1:56572/147‑1 (MQ=255) aaCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTTGTTTTTTTTTTTCTTAAAAG‑AAAGCGGAGTA‑Gccccaaactaaccg > 2:510470/1‑137 (MQ=255) aaCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTTGTTTGTTTTTTTTCTAAACG‑AAAGCGGAATA‑GCCATAATCaaaaag > 1:909854/1‑146 (MQ=255) aGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTTGTTTTTTTTTTTTCTAAAAA‑AAGGGGAAGGA‑AACAAAAAAAAACGGGACg > 2:396110/1‑150 (MQ=255) gACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTAGTTTGTTTTTTTCTTAAAAG‑AAAGGGAAGTA‑GGCACAATCaaaacg > 2:634012/1‑142 (MQ=255) cccGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTTGTTTGTTTTTTTTATAAAAG‑AAGGCGAAGTA‑GGCATAATCAAAAGGGACGAAAGAGTcct > 1:404044/1‑149 (MQ=255) cgcgAATCCGGATAAACCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTCTCAGTTTGTTCTTGTCATAACCG‑AAGGCGAATTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCgg < 2:647163/151‑1 (MQ=255) ggATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTTTGTTTTTTTTTGTCCTAAAAG‑AAGGCGGAGGA‑GCCATAATCAAAAGGGAAGAAAAAGGCATTATTTCAGCCGCCcgg > 1:483107/1‑149 (MQ=255) gATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTATTTTTTTTTTTTTTAAAAAG‑AAGGGGGAGGA‑GCCCTAAACAAAAGGTACGAAAGGGGCATTATTTCATCCGCCggg > 1:635820/1‑148 (MQ=255) cccccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATcc < 2:517880/149‑1 (MQ=255) cAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑GTTTTTTGTTTGTTTTTTTTTTAAAAGAAAGGGGAAG‑A‑GCCCTAATCAAAAGGGTCGGAAATGTCATTTTTTCAGCCCCCGGGGACGGTTAAGGCATccccctcc > 1:804570/1‑144 (MQ=255) ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTTTGTTTGTTTTTTTCATAAAAA‑AAGGCGAAGAA‑GCCATAATCAAAAGGGACGAAAGGGTTATTATTTCAGCCGCCCGGGACGGTTAAGGCATCCGCCGCCACA‑ccc > 1:696532/1‑150 (MQ=255) cAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACA‑CCGCTTTTGGCaa < 1:478636/150‑1 (MQ=255) aGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTTTTGTTATAAACC‑AAGGCGGAGGAGGCCCCAATAAAACGG‑ACGAAAAGGGCCTTATTTCTGCCGCCCGGGACGTTAAAGGGACACGCCACCCCC‑CCCCCTTTTGCAAACaag > 1:564409/1‑148 (MQ=255) gTTCTTTCAGGCTTTGTGGCGGG‑TTTTTTAGTTTGTTTTTGTCATAAAAG‑AAGGGGGAGGA‑GCCATAATAAAAAGGGACGAAAATGGCATTTTTTCAACCCCCCGGGAAGGTAAAGGCATCCGCCGCCCCACCCCCTTTTGCCaaaaatcc > 1:930179/1‑146 (MQ=255) gCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACA‑CCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTc > 1:39603/1‑151 (MQ=255) cGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACA‑CCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt < 2:401705/151‑1 (MQ=255) g‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACA‑CCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc > 1:935683/1‑151 (MQ=255) | GGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTA‑GCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACA‑CCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAC > NC_000913/74935‑75223 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |