| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 42.9% | 54.2 / 17.4 | 43 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (21/3); new base A (0/18); total (21/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.52e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC > NC_000913/436368‑436662 | ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGc < 1:832434/149‑1 (MQ=255) cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATg > 2:413384/1‑150 (MQ=255) aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCt > 2:705706/1‑151 (MQ=255) gcgctgGGCGTTCTCTGGCGGTGAAGATTACGAGTTTTTTTTCCCTGTGCC‑GGAACTAAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTCCCTTTTCCC‑TGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGttt < 2:740946/150‑1 (MQ=255) gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTAt > 2:325368/1‑151 (MQ=255) gggCGCTCTGTGGGGGTGAAGATTACGAGTTGTTTTTCATTGTGCC‑GAAACAGAC‑CGGTGGCGCGGTGGGTGTGGCTCTGGGACCCCGGGGCGCCCCTTTTCCC‑TGTATCGGCA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTa < 2:796503/147‑1 (MQ=255) gggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg > 1:755320/1‑151 (MQ=255) gggCGCTCTCGGGCGTTGAAGATTACGGGTTGTGTTTCACTGTGCC‑GGACCTGAC‑CCGGGGCGCGGTGGATGTGCCTCTCGGACACTGGGGCGTCCCTTTTACC‑TTTGCCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTg < 1:874225/142‑1 (MQ=255) gcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc > 1:798530/1‑151 (MQ=255) ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc > 2:441357/1‑151 (MQ=255) tcggggggtaaggATTCCGATTTGTGTTTCTTGGTG‑C‑GGGACTGAACCCGGGGCGCGTTGGATGTGGCTCTGGGACACCTGGGCGTCCCGTTTCCC‑TGTTTCGGGA‑AAATGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa < 2:540255/134‑1 (MQ=255) gggggTGAAGATTACGAGTTGTGTTTCACTGTGCG‑GGAACTGAA‑CGGTGGCGCGCTGGATGTGGCTCTGGGACCCGTGGGCGTCCCTTTTCCC‑TGTATGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc < 2:567515/148‑1 (MQ=255) ggCGGGGAATATTACGAGTTGTTTTTCAC‑GTGCCTGGAAGTAAA‑CCGTGGTGAGTGGGTGGTGCCTCTGGAACACCTGGGGTTTCTTTTTTCC‑TTTACGGGGA‑AAATCACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg < 1:362731/146‑1 (MQ=255) gCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc > 1:905710/1‑151 (MQ=255) ggggAAGATTACGTGTTGTTTTCCCCTGTGCC‑GGACATGAA‑CCGTGGCGCCTTGGATGTGCCTCTCGGCCCCCGGGGCGTCCCGTTTACC‑TTTATGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg < 2:357732/148‑1 (MQ=255) aaGATTGCGAGTTGTGTTTCACTGTGCC‑GGAACTGAC‑CCGGGGCGCGCTGGTTGTGGCTCTCGGACACCTGGGCGTCCCGTTTCCC‑TGTATGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa < 2:565178/150‑1 (MQ=255) ttACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt > 1:642226/1‑150 (MQ=255) ttGTTTTTCGCTGTGCG‑GGATCTAAC‑CGGGGGCGCGTTGGTTGTGGCTCGCGGACCCCTGGGCGTCCCTTTTCCC‑TGTATCGGGC‑AAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg < 1:64216/150‑1 (MQ=255) ttGTTTTCCACTGTGCC‑GGAACTGAC‑CGGTGGCGCGTTGGATGTGGCTCGGGGCCCCCGGGGC‑TTCCGTTTTCCTTGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa < 1:354321/146‑1 (MQ=255) tgtgTTTCCTTGTGCC‑GGAACTGAC‑CGGTGGCGCGCTGGTTGTGGCTCTCGGACCCGGGGGCGTACCGTTTACC‑TTTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa < 2:761091/150‑1 (MQ=255) tgtTTCACTGTGCC‑GGACCTGAC‑CCGGGGCGCGTTGGATGTGGCTCTCGGCCCCCTGGGCTTCCCTTTTCCC‑TGTATGGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa < 1:772711/150‑1 (MQ=255) tgtTTCACGGTGCC‑GGAATTGAC‑CCGGGCCGCGCTGGTTGTGGCCCCCGGCCCCTTGGGCGTCCTTTTTCCC‑TGTACGGGGA‑AAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa < 1:926415/148‑1 (MQ=255) cACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat > 1:833634/1‑151 (MQ=255) tgtgCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc > 2:690051/1‑151 (MQ=255) gtgCG‑GAAACGGAC‑CGGTGGGGCGGTGGATGGGCCTCTCGGACCCGGGGGCGTACCTTTTACT‑GGTATGGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa < 2:770861/144‑1 (MQ=255) gtgCC‑GAAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCGTTTACT‑TGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc < 1:433396/150‑1 (MQ=255) ggAACTGAC‑CCGGGGCGCGCTGTGTGGGGCTCTGGAACCCGGGGGTTTACCTTTTGCT‑TGTATCGGGT‑AATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt < 2:858299/151‑1 (MQ=255) ggAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt > 2:95682/1‑151 (MQ=255) gAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg > 1:184262/1‑151 (MQ=255) aa‑CCGGGGCGCGTTGGATGTGGCTCTCGGACACCTGGGCGTCCCTTTTACT‑TGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg < 2:550864/150‑1 (MQ=255) gTGGCGCGTTGGTTGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCC‑TGTAGCGGAC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc < 1:407455/148‑1 (MQ=255) gTGGCGCGGTGGTTGGGGCTCCCGAACCCCTGGGCGTCCTTTTTCCC‑TGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa < 1:320807/151‑1 (MQ=255) tGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa > 1:548903/1‑150 (MQ=255) ggCGCGCTGGTTGGGGCTCTGGGACACTGGGGCGTACCTTTTCCC‑TGTTTGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa < 2:465402/149‑1 (MQ=255) tGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg > 2:318889/1‑151 (MQ=255) ggCTCTGGGACCTCGGGGCGTCCCTTTTCCT‑TGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa < 2:698149/150‑1 (MQ=255) ggCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa > 2:99624/1‑149 (MQ=255) gcctcccgaccccgtGGGCGTCCCTTTTCCC‑TGT‑GCGGGCAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg < 1:355254/132‑1 (MQ=255) tcGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc > 2:165640/1‑151 (MQ=255) gCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt > 2:738658/1‑150 (MQ=255) cc‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc > 2:652083/1‑151 (MQ=255) cc‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc > 2:828410/1‑151 (MQ=255) tATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt > 1:350680/1‑151 (MQ=255) cGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc > 2:487257/1‑151 (MQ=255) | TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCC‑GGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC > 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: ‑ |