| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,497 | 0 | A | C | 40.8% | 74.9 / 11.5 | 53 | V288V (GTA→GTC) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base A (20/9); new base C (0/20); total (20/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.04e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGT > NC_000913/436348‑436642 | ctgcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTAc > 2:44286/1‑151 (MQ=255) gcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTAc < 2:830794/149‑1 (MQ=255) gcCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCg < 1:381834/151‑1 (MQ=255) cGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGtt < 2:411900/151‑1 (MQ=255) cGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGtt < 1:854142/151‑1 (MQ=255) tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTAc > 2:591909/1‑151 (MQ=255) tCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTAcc > 1:956010/1‑150 (MQ=255) tGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGt > 2:898148/1‑151 (MQ=255) gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGt < 2:83600/150‑1 (MQ=255) gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGt < 1:480098/150‑1 (MQ=255) gcgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTCGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGTCACCT‑GGGCGTACCGTTTACC‑TGTATCg < 1:251481/151‑1 (MQ=255) gcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCggg > 1:622253/1‑151 (MQ=255) 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) | CTGCCGTTTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGT‑GGCGCGCTGGATGTGGCTCTCGGACACCT‑GGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATC‑GAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGT > NC_000913/436348‑436642 |
| 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: ‑ |