| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 54.5% | 8.7 / 28.4 | 23 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (7/3); new base A (0/12); total (8/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.04e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-04 | |||||||||||
TCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT > NC_000913/436371‑436660 | tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATg > 2:137784/1‑151 (MQ=255) cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTCCCGTTTACC‑TGTATCGGGCaaa < 1:91084/148‑1 (MQ=255) gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgcc < 1:74997/150‑1 (MQ=255) ctgGGCGCTCTCTGGGGGTGAAGATTAGGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGTGAGCTGGATGTGGCTCTTGGAAAACTGGGCGTAACGTTTAGC‑TGTATTGGGGAAATGAACGGCGATATGGAAGGGGTTTGTTTTAtt > 1:438109/1‑151 (MQ=255) gggCGCTCTGTGGCGGTGAAGATTGCGAGTTGTTTTTCCCTGTGCCGGAACTGAGCCGTGGCGCGCTGGATGTGGCTCTCGCACCCCGGGGCGTGCCGTTTCCC‑TGTCTCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg < 2:123430/151‑1 (MQ=255) ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc > 1:538508/1‑151 (MQ=255) gggggTGAAGATTGGGAGTTGTGTTTCGCTTGGCTGGACCTGAACCTTGGCGGGGTGATTTTGTCTCCCGGACCTGTGGGCTTTCCTTTTTCC‑TTTTTCGGAAAAGTGCCCGCCGATATCGAATGGCTTTGTTTTATTCGTGACGGCGa < 1:426167/146‑1 (MQ=255) aaGATTAGGATTTGTGTTTCTCTGTGCCGGAACTGAGCCGTGGCGCGCTGGGTGTGGCTCTGGGACCCTTGGGTGTCCTTTTTTCT‑TATTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa < 2:333864/150‑1 (MQ=255) ttACGATTTTTTTTTCTCTGGGCCGGACCTGACCCGTGGCGGGTTGGGGGGGCCTCTCGGACCCGGGGGCTTCCCGTTTCCC‑TGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt < 2:483078/150‑1 (MQ=255) ttttgtTTCTCTGTGCCGAAACTGAGCCGTGGCGCGCTGGTTGTGCCCCCCGGACCCCTGGGCGTCCTTTTTTCC‑TGTACCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg < 1:48817/146‑1 (MQ=255) ttGTGTTCCCTTGTGCCGAAACTAACCCGTGGCGCGTGGGTTGTGGCTCTCGAACCCCGGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg < 2:87253/145‑1 (MQ=255) cactgtcagggAACTGAGCCGGGGCGGGTTGGAGGGGCCCCTCGGCCTCCGGGGCGTCCCTTTTTCC‑TGTTCAGGGAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat < 1:23871/142‑1 (MQ=255) gtgCCGGAACTGACCGGTGGGGCGGTGGATGTGGCTCTCGGCCCCTTGGGCGTCCTTTTTCCC‑TGTAGCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg < 1:300230/148‑1 (MQ=255) ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt > 2:268502/1‑151 (MQ=255) aCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc > 1:402472/1‑151 (MQ=255) ggATGTGCCTCGCGGACCCCGGGGCGTCCCTTTTCCT‑TGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAg < 2:314025/144‑1 (MQ=255) gtggtgtcttcGGCCCCCTGGGGGGTCCGTTTTCCTTGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg < 1:456455/142‑1 (MQ=255) tctcGGCCACGGGGGCTTCCCTTTTCCT‑TGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa < 2:478590/145‑1 (MQ=255) tctcGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg > 2:374254/1‑151 (MQ=255) gggCGTCCCTTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa < 2:241742/147‑1 (MQ=255) tACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc > 2:39594/1‑151 (MQ=255) ttACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc > 2:62444/1‑151 (MQ=255) gggAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt < 2:299629/148‑1 (MQ=255) | TCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT > NC_000913/436371‑436660 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |