| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 43.5% | 55.5 / 21.5 | 46 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (18/8); new base A (0/20); total (18/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.65e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
CGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT > NC_000913/436372‑436660 | cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGAGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATg > 1:366766/1‑150 (MQ=255) tGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgcc > 2:332020/1‑151 (MQ=255) gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgccg > 2:484883/1‑151 (MQ=255) ggcgcggcgcGGGGCGCTCTCTGGGGGTGAAGATTACGAGTTGTGTTTCCCTGGGCCGGAACTGACCCGG‑GGCGCGTTGGTTGTGGCTCTCGGCCCCCG‑GGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATC‑GAAGGGCt < 2:933337/150‑1 (MQ=255) gcgcggGGCGCTCTCGGGGGGTGAAGTTTACGATTTGTTTTTCCCTGTGCCGGAACTGACCCGT‑GGCGCGGTGGATGTGGCTCCCGGACACGT‑GGGGGTCCTTTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGt < 1:277354/148‑1 (MQ=255) gggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGACCCGT‑GGCGCGCTGGATGTGGCTCTCGGACACTT‑GGGCGTCCCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTAtt < 2:797845/149‑1 (MQ=255) ggCGCTCTCGGGGGGTGAAGTTTAGGAGTTGTGTTTCACTGGGCGGGAACTGAACCGG‑GGCGCGCTGGATGTGGCTCCCGCCCACGT‑GGGAGTCCTTTTTCCT‑TGTATCGGGAAAATAACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGt < 1:95094/151‑1 (MQ=255) gcgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTg > 2:595571/1‑151 (MQ=255) tctcgGGGGGTGAAGATTACGGGTTGTGTTTCACTGTGCGGGAACTGAACGGG‑GGCGCGCGGGAGGGGGCTCCCGGACCCGG‑GGGGGTCCCGTTTCCT‑TGTAGCGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGa < 2:95389/148‑1 (MQ=255) ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGc > 2:565140/1‑151 (MQ=255) tggggCGGGGATAATTACGGTTTGTGTTTATCTGT‑CCGGGAATGAAGCGTGGGCGGGTGGTTTGTGGCTCCCGATCACCT‑GGGCGATCCGTTATTC‑GCTTTCGGGAAAGGGCCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGc < 2:251502/147‑1 (MQ=255) ggggggTGAAGATTACGAGTTGTGTTCCCCTGTGCCGAAACTGACCCGG‑GGCGCGTTGGATGTGGCTCCCGGACCCGG‑GGGCGTACCTTTTACC‑TGTAGCGGGAAAAGGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACgg < 1:198024/143‑1 (MQ=255) gggggTAAAGATTACGAGTTGTTTTTCACGGTGCCGGAACTGACCGGG‑GGCGCGGTGGATGTGGCTCGCGGACACCT‑GGGGGTACCTTTTACC‑TGTATCGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGc < 2:834310/144‑1 (MQ=255) ggCGGGGAAGATTACGATTTGTGTTTCGCGGGGCCGAACTTGACCCGT‑GGCGCGGTGGATGGGGCTCGCGGCCACCG‑GGGCGTTCTGTTTTCC‑TGTAGGGGGCAAAGGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCg < 1:865403/148‑1 (MQ=255) ggggAAGATTACGGGTTTGGTTTCTCGGGGCGGGACCTGACCCGG‑GGCGGGCTGGTTGTGGCTCCCGGCCCCGG‑GGGCGTCCTTTTTCCT‑TGTACCGGGCAAATAACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGaa < 1:741121/144‑1 (MQ=255) aaGATTACGTGTTGTGTTTCCCTGTGCCGGAACTGAACGGG‑GGCGCGCGGGTTGTGGCTCTCGGACACCG‑GGGCGTCCCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa < 2:656092/150‑1 (MQ=255) ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt > 2:74953/1‑150 (MQ=255) tACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt > 2:578459/1‑149 (MQ=255) gggTTGTGTTCCCGGGTGCGGGAACGGACCGGG‑GGCGCGGTGGTTGTGGTTTTGGGACACC‑‑GGGGGTTCCTTTTTCCTTGTAGGGGGAAAATACCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg < 1:611195/146‑1 (MQ=255) gtTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg > 2:287323/1‑151 (MQ=255) tttCGCTGTGCGGGACGTGAACCGG‑GGCGCGCTGTATGGGGCTCCCGGACCCCG‑GGGGGTCCCGTTTACC‑TGTATGGGGAAAATGCCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa < 1:298101/146‑1 (MQ=255) tttCACTGTGCCGGAACGGAACCGT‑GGCGCGCGGGATGTGGCTCTCGGACCCCG‑GGGCGTCCCGTTTACC‑TGTATCGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg < 2:814798/149‑1 (MQ=255) ttccccTGTGCCGGAAGTGAGCGGT‑GCCGCGCTGGTGTTGGCTCCCGGACCCTT‑GGGCTTCCCTTTTCCT‑TTTATGGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa < 2:360681/142‑1 (MQ=255) gCCGGAACGGACCGGG‑GGCGCGTTGGAGGTGGCTCTGGGACCCGG‑GGGCGTCCCTTTTCCT‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc < 2:311/148‑1 (MQ=255) ccGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt > 2:349776/1‑151 (MQ=255) cGGAACTGAACCGT‑GGCGCGCTGATTGTGGCTCGCGGACACGT‑GGGCGTCCCGTTTCCC‑TTTATCGGGCAAATAACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt < 1:185084/149‑1 (MQ=255) acctaaccgggggCGCGTTGGTGGGGGCTCTCGCACCCGT‑GGGCGTCCCTTTTCCT‑TGTAGGGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg < 1:305721/139‑1 (MQ=255) cTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc > 2:369350/1‑150 (MQ=255) acccGT‑GGCGCGGTGGATGTGGCTCTGGGACCCGG‑GGGCTTCCCTTTTCCT‑TTTATGGGCAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc < 2:781860/149‑1 (MQ=255) acccGG‑GGCGCGCTGGATGTGGCTCCCGGACCCGG‑GGGCGTCCCGTTTCCC‑TGTTTGGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg < 2:465408/143‑1 (MQ=255) aaCCGT‑GGCGCGCTGGATGTGGTTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc > 1:93417/1‑151 (MQ=255) cGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt > 2:280169/1‑151 (MQ=255) gT‑GGCGCGCTGGATGTGCCTCTCGGCCACCT‑GGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt < 2:221322/150‑1 (MQ=255) ggggcgtTGGTTGTGGCTCGC‑GCCGCCTGGGGCGTCCCTTTTTCC‑TGTAGCGGGAAAATGCCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg < 1:787778/141‑1 (MQ=255) tGGATGTGGCTCTCGGACACTT‑GGGCGTACCGTTTACC‑TGTATCGGCCAAACGCCCGCCG‑TACCGGACGGGCTTTGTGTTATTCCTGTCGGCGAACCCGTTACATTTTACTGGACAGGCTATGCCCATTTGGCCACGCAATAAAGTTGt > 1:23820/1‑149 (MQ=255) tgtgGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGa > 2:1007505/1‑150 (MQ=255) tgGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa > 2:336721/1‑150 (MQ=255) ggACACTG‑GGGCGTGCGTTTTCCC‑TGTGGGGGGAAAGTGCCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt < 1:116531/150‑1 (MQ=255) ccT‑GGGTGTCCTTTTTCCC‑TGTGGCGGGCAAATGCCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt < 1:44258/151‑1 (MQ=255) gggCGTCCTTTTTCCT‑TGTATCGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt < 2:315720/148‑1 (MQ=255) ggCGTCCCGTTTCCC‑TGTAGGGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt < 2:467764/151‑1 (MQ=255) cGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGt > 2:677625/1‑151 (MQ=255) tttCCC‑TGTAGCGGCAAAATGCCCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg < 2:163844/150‑1 (MQ=255) tttccttGTATGGGGAAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTaa < 2:93407/138‑1 (MQ=255) tACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc > 1:50529/1‑150 (MQ=255) taccgttgcagaacAAATAACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc < 1:258487/137‑1 (MQ=255) gggCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAAACTGTTACATTAGACTGGAAAGGATATAACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGCGGCATCTACtt > 2:131574/1‑148 (MQ=255) | CGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT > NC_000913/436372‑436660 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |