| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 38.2% | 39.7 / 10.6 | 34 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (17/4); new base A (0/13); total (17/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.13e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 | |||||||||||
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT > NC_000913/436367‑436663 | tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGc > 1:611450/1‑150 (MQ=255) ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGa > 2:435379/1‑150 (MQ=255) gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGt > 1:608367/1‑151 (MQ=255) tgcggggggCGCGCTCGGGCTGTAAAGATTACGGGTT‑TGTTGTCACTTTGCCGAAAGTGAACCGTGCCGCGGGGGCTGTGTC‑CCCGGCCCCCTGGGCCGTACC‑TTTTCCTGTTTGCGGCAAAATGACCGCCGATATCGAAGGGCTTTGt < 1:246984/143‑1 (MQ=255) ctgGGCGCTCTCTGGCGGTAAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCTTTTACCTGCAT‑CGGGCAAATGACCGCCGAGTTCGAAAGGCTTTGtttt > 1:405838/1‑148 (MQ=255) gggCGCTCTCGGGGGGTGAAGATTACGGGTTGTGTT‑TCCCTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGG‑CGTCCCCTTTCCCTGTAT‑CGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc < 2:678636/150‑1 (MQ=255) gggCGCTCTCGGGCGGTGAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGG‑CGTCCGTTTTACCTGTGT‑CGGAAAAATGACCGCCGATATCGAAGGGCTTTg < 2:453825/142‑1 (MQ=255) ggggggaaaGATTAGGTTTTTTTTT‑TCACTGGCCGGGAAAGGACCCGTGGCGGGCGGGTGGGGGCCCTGGGACCGCGGGG‑GGCGCGTTTTTCTTGTTG‑CGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa < 1:157750/143‑1 (MQ=255) aaGATTACGAGTTGTGTT‑TAACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa > 1:74563/1‑150 (MQ=255) ttACGAGTTGTGTT‑TCACTGTGCGGGAACTGACCGGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGG‑CGTACCGTTTCCCTTTTT‑CGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt < 2:95485/149‑1 (MQ=255) tACGGGTTGTTTT‑TCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGG‑CGTCCCTTTTCCCTGTAT‑CGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa < 2:149436/150‑1 (MQ=255) cGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc > 1:416673/1‑151 (MQ=255) cGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc > 2:230287/1‑151 (MQ=255) gAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa > 1:445875/1‑149 (MQ=255) ttttgtT‑TCGTTGTGCCGGAACTGACCCGTGGCGCGGTGGATGTGCCTCTGGGACCCTTGGG‑CTTTCCTTTTCCCGGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg < 1:521402/147‑1 (MQ=255) tgttgcggtgtgCGGGAACGGAACGGTGGGTCGGTGTATGTGGCTCTGGGCCCTTTGTG‑TGTTCTTTTTCCTTATGT‑CGGGCAGTTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa < 1:605962/141‑1 (MQ=255) t‑tCGCTGGGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCGCGGCCACCTGGG‑CGTCCTGTTTCCTTGTAG‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa < 1:531955/150‑1 (MQ=255) ccGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGCTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt > 1:53702/1‑151 (MQ=255) aaCTGACCCGTGGCGCGGGGGGTGGGGCTGTCGGCCCCGGGGG‑CGTCCCGTTTCCTTGTTG‑CGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt < 1:14518/149‑1 (MQ=255) aaCGGAACCGGGGCGCTTTGGTTGTGGTTCTCGGACCCCTGGG‑CTTCCCTTTTGCTTATAG‑CGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt < 1:671725/148‑1 (MQ=255) cTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGGTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc > 1:151038/1‑150 (MQ=255) gcgcgcTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa > 2:665601/1‑149 (MQ=255) gcgcTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa > 2:81087/1‑147 (MQ=255) cTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt > 2:398271/1‑150 (MQ=255) ggAGGTGGCTCTCGGCCCCCGGGG‑CGTCCTTTTTCCCTTTCT‑GGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg < 2:638481/150‑1 (MQ=255) tgGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa > 2:90008/1‑150 (MQ=255) ggctctgcgccccggggg‑TGTACCTTTTCCCTGTGA‑GGGAAAAGTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACGTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg < 1:211891/136‑1 (MQ=255) ggg‑CGTCCCTTTTACCTTTAG‑GGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAg < 2:608856/151‑1 (MQ=255) gTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg > 1:296336/1‑151 (MQ=255) tttccttttagggggCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc < 1:31111/134‑1 (MQ=255) tGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGACGAACCTGTTACATTAGACTGGAAAGGATATGACCATATTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc > 1:87994/1‑149 (MQ=255) tAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt > 2:185067/1‑151 (MQ=255) aT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGACGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt > 1:396218/1‑151 (MQ=255) gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt < 2:385622/148‑1 (MQ=255) gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt < 1:39204/151‑1 (MQ=255) | TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTT‑TCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACCTGTAT‑CGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT > NC_000913/436367‑436663 |
| 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: ‑ |