| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 44.0% | 32.1 / 10.6 | 25 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (13/1); new base A (0/11); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.69e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 | |||||||||||
TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC > NC_000913/436368‑436662 | ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCa < 2:13247/150‑1 (MQ=255) ccGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCg > 2:384135/1‑151 (MQ=255) gcgctgcgctGGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCC‑GTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCttt < 2:252397/151‑1 (MQ=255) tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt > 2:203385/1‑151 (MQ=255) tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc > 1:637919/1‑151 (MQ=255) cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg > 2:329390/1‑150 (MQ=255) tctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg > 2:460973/1‑151 (MQ=255) gggggTGAAGATATGGAGTTGTTTTTCTTTGTCCGGGATCTAACCCGGGGCGCGTTGGTGGTGCCCCTCGGAC‑CCTGGGGTTTCCTTTTTCCCTTTATGGGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa < 2:252535/147‑1 (MQ=255) gCGTTGAAGTTTACGGGTTGTGTTTCCCTGTGCCGGACCTAACCCGTGGCGGGTTGGTTGTGGCTCTCGGACACCGGGGCGTACC‑TTTTCCCTTTAGCGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc < 1:167015/151‑1 (MQ=255) gattcggAGTTGTGTTTCCCTGGGCCGGAACTGAACCGGGGGGCGGTGGTTGTGCCTCTCGCACACGTGGGCGTCCC‑TTTTCCCGGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa < 1:10254/142‑1 (MQ=255) gTTGTGTTTCACTGTGCCGGAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATCTCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt > 2:370799/1‑149 (MQ=255) gTTGTGTTCCGGTGTGCCGAAATTGAACCGTGCCGCGCTGGGTGTGGTTCCCGGACCCCTGGGCGTACC‑GTTTCCCTTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc < 2:156978/148‑1 (MQ=255) tgttTTTCACTGTGCGGAAACTGACCGGTGGCGCGGTGGTTGTGGCTCTCGAACCCCGGGGCGTCCC‑TTTTACCTGTATGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg < 1:516217/149‑1 (MQ=255) cgggACCTGACCCGGGGCGCGCTGGATGTGCCTCTCGGACACGTGGGGGTCCC‑TTTTCCGTTTTTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat < 1:592611/142‑1 (MQ=255) gTGGCGCGCTGGATGGGGCTCTCGGCCCCCGGGGCGTACC‑TTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa < 2:368458/149‑1 (MQ=255) gcgcTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa > 2:482453/1‑147 (MQ=255) gtggtgtgggcctctCGCCCCCCGGGGCGTCCC‑TTTTCCTTGTAGGGGGAAAGTGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt < 1:459968/138‑1 (MQ=255) gcTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg > 2:12678/1‑150 (MQ=255) gtgGCTCTCGGCCGCGGGGGCGTCCT‑TTTTCCCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGa < 1:199976/149‑1 (MQ=255) gCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt > 2:283068/1‑151 (MQ=255) tcGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc > 1:541755/1‑151 (MQ=255) acacCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa > 1:397774/1‑150 (MQ=255) ccctggggCGTCCC‑TTTTCCCTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAg < 2:395770/146‑1 (MQ=255) tATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt > 2:362048/1‑151 (MQ=255) cGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc > 2:232126/1‑151 (MQ=255) | TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC > NC_000913/436368‑436662 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |