| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 75,079 | 0 | C | T | 60.0% | 5.3 / 27.4 | 25 | L134L (CTG→CTA) | thiB | thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein |
| Reads supporting (aligned to +/- strand): ref base C (1/9); new base T (15/0); total (16/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-05 | |||||||||||
GTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAC > NC_000913/74940‑75223 | gtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTtttttttt > 2:792022/1‑151 (MQ=255) cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTCTTAAAAGaa > 1:605715/1‑151 (MQ=255) ccATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTTTTATTAAAGaaag > 2:682215/1‑149 (MQ=255) cTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTa < 1:653586/151‑1 (MQ=255) cccATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGTTCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑Tc < 1:128509/150‑1 (MQ=255) aCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTAAAAGAAAGCGGAGAAGCCATAA‑TC‑AAAAggga > 2:150794/1‑148 (MQ=255) aaCCCAGCCCCGGTGTGCTGGAGCGCGGATCTTAATAAATCACCCGCCAATTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGTGTTTGTCAGTTTGTTCTTGTCATAAACGAATGCGAATTAGCCATAA‑TC‑AGACGGTACGaa < 1:662414/150‑1 (MQ=255) aCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAATTTTTTTTTTTCTTAAAAAAAGAGCAAGTAGCCCTAATTC‑CAACGG‑ACGaaa > 1:742637/1‑151 (MQ=255) cccAGCCCGGGTGTACTGGTGCGCGGTTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAg < 1:48846/151‑1 (MQ=255) cccAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAg < 2:323803/151‑1 (MQ=255) ccGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGGTTGTTTTTTTCATAAAAGAAGGGGAAGGAGGCATAA‑TC‑AAAAGGGAAGaaaagggcatt > 2:61988/1‑143 (MQ=255) gtgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTCATAAACGAAGGGGAAGTAGCCCTAA‑TC‑AAAAGGTACGAAAGAGTCattatt > 2:375833/1‑151 (MQ=255) gtgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTTTTTTTTTCATTAACGAAAGGGGAGTAGCCCTAA‑ACAAAACGGTAGGAAAG‑GTCattttt > 2:352785/1‑151 (MQ=255) gtgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAGTTTGTTTTTTTCATAAAAGAAGGCGAAGTAGCCATAA‑TC‑AAAAGGTACGAAAGTGTCattatt > 2:792140/1‑151 (MQ=255) gtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTTTCTTAAACAAAAGGGAAGGAGGCCTAA‑TC‑AAAAGGGACGAAAGAGTCattattt > 2:742254/1‑149 (MQ=255) gtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTCTTTTCATAAAAGAAGGCGAAGTAGGCATAA‑TC‑AAAAGGGACGAAAGAGTCatttttt > 2:778635/1‑149 (MQ=255) gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGGTTTTTCCTTTCATAAAACAAGGCGGAGGAGCCATAA‑TC‑AAACAGGACGAAAGTGTGATTATTTCAGcccccg > 2:725188/1‑151 (MQ=255) gATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTTTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGGACGAAAGTGTCATTATTTCAGCCGCCGGGaaa > 2:712521/1‑150 (MQ=255) cTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTCATTAAAAAAGGCGGAGTAGGCATAA‑TC‑AAACAGGACGAAAGGGTCATTATTTCAGCCGCCgggggaggta > 1:159251/1‑144 (MQ=255) tAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACgg < 2:477790/151‑1 (MQ=255) tctcAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGc < 1:594094/151‑1 (MQ=255) gTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTGTTTTTTTTATAAAAGAAGGCGAAGGAGCTATAA‑TC‑AAACGGGACGAAAGGGTTATTATTTCAACCGCCGGGGAAGGTAACGGCATCCGCCGCCCCACCGCCTTTTGCAAAAAGTc > 2:295851/1‑150 (MQ=255) cTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTTTTTTTTTTATAAAACAAGGGGAAAAAGGCATAA‑TT‑AAAAGGTAAGAAAGTGTCATTTTTTCAGGCGCCCGGGAAGGTAACAGGATCCCCCGCCCCACCCCTTTTTGCCAAAAATCCggg > 1:371218/1‑150 (MQ=255) ttCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGttt < 2:483353/151‑1 (MQ=255) tttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc < 1:150794/150‑1 (MQ=255) | GTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAA‑TC‑AAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAC > NC_000913/74940‑75223 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |