Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 436,516 | 0 | C | A | 41.5% | 66.0 / 27.7 | 54 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
Reads supporting (aligned to +/- strand): ref base C (27/4); new base A (0/22); total (27/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.23e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01 |
CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC > NC_000913/436366‑436662 | cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGc > 1:418175/1‑151 (MQ=255) tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGAAACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑aa > 2:220474/1‑150 (MQ=255) cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑Ga < 2:28845/151‑1 (MQ=255) ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GAcc < 2:53046/151‑1 (MQ=255) gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GAccgcc < 1:813288/150‑1 (MQ=255) cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCg > 2:740228/1‑150 (MQ=255) ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGa > 1:654476/1‑150 (MQ=255) ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑GACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGa > 2:692103/1‑150 (MQ=255) ggggcggcggggggCGCTCTCGGGGGGAGAAGATTGCGGTTTGTTTTTCCCTGTCCCGA‑AACTGACCCGGGGCGCGTT‑GGGGGTGCCTCCCGCCCCCTG‑GGGGG‑TCCT‑T‑TTTACCTTTTGGGGAA‑AAAT‑GACCGCCGATATCGAAGGGCtt < 2:1180320/140‑1 (MQ=255) gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCttt > 1:402729/1‑151 (MQ=255) cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGtttt > 1:979128/1‑150 (MQ=255) cgcggGGCGCCCTCGGGGGGTGAAGATTCGGGTTTGTTTTTCCCTGTGCCGG‑AACTGACCCGGGGCGCGCT‑GGATGTGGCTCCCGGACCCTG‑GGGCG‑TCCT‑T‑TTTTCTTGTACCGGAA‑AATT‑GACCGCCGATATCGAAGGGCTTTGtttt < 2:1106924/150‑1 (MQ=255) ggtgGGCCTTCTCGGGCGGTGAAGTTTCCGATTTGTGTTTCCCTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCCCTCGGCCACCT‑GGGGG‑TCCC‑T‑TTTCCTTTTATGGGGA‑AAAA‑GACCGCCGATATCGAAGGGCTTTGTTTTAt < 1:817194/149‑1 (MQ=255) gggcctctgcggggggTGAAGATTAGGGTTTTTGTTTCCCTGGGCCGG‑AACTGAACGGTGGCGCGTT‑GGTTGTGCCCCGCGCAC‑CCG‑GGGGG‑TTCCTT‑TTTCCCTTTAGCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAtt < 2:1078737/135‑1 (MQ=255) gggCGCTCTGTGGCGGGAAAGTTTACGTGTTTTTTTTCTGTGTGCGGG‑AACTGACCGGGGGCGCGCT‑GGGTGTGGCTCTCGGACCCTG‑GGGCG‑TCCT‑G‑TTTCCCTGTGGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCttt < 2:254210/141‑1 (MQ=255) cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg > 1:1168517/1‑150 (MQ=255) ctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc > 1:804978/1‑149 (MQ=255) ggggttGAAGATTCGGATTTGTTTTGTATGGTCCCGG‑AACTGCTCCGGGGCGCGCG‑GGTTGTCTCCTTCGG‑CGACG‑GGGCT‑TTCC‑GTTTTCCCTTTGTCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa < 2:784792/145‑1 (MQ=255) ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGACCCGGGGCGCGTT‑GGTTGTGGCTCTCGGACACTG‑GGGCG‑TACC‑G‑TTTCCTTTTATGGGGA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc < 1:1032022/147‑1 (MQ=255) ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa > 2:203604/1‑150 (MQ=255) ggtgaagtttcgggtttGTGTTTCACTGTGCGGG‑AACTGACCGGGGGCGCGTG‑GGATGGGTCTCTCGGACCCCG‑GGGCG‑TCCT‑T‑TTTCCTTGTAGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg < 2:815475/136‑1 (MQ=255) ggTAAAAATTTCGATTTGTTTTTCACTGTGCGAG‑AACTGACCTGCGGCCCGTG‑GGGTGGGGCTCGCGGCCACGT‑GGGCG‑TACT‑T‑TTTACCTGTACCGGGA‑AAGT‑CACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt < 2:349134/150‑1 (MQ=255) aGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc > 2:791576/1‑150 (MQ=255) taggAGTTTTTTTTCCCTGGGCCGG‑AACTGACCGGGGGCGCGTT‑GGTTGTGGCTCTGGGACCCCG‑GGGCG‑TCCC‑T‑TTTCCTTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa < 1:1050362/147‑1 (MQ=255) tacgttttGTGTTTCCCTGTGCCGG‑AACGGACCGGGGGCGCGTT‑GGTGGTGCCTCCGGGACACCT‑GGGCT‑GCCC‑T‑TTTCCCTTTGCGGGAA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa < 2:435075/144‑1 (MQ=255) tACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGCCACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt > 2:2904/1‑149 (MQ=255) cGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc > 2:140274/1‑151 (MQ=255) ggttgttTTTCCTTGTGCCGA‑AACGGACCCGTGGCGCGGT‑GGGTGGGGCTC‑CGGGAACCTGGGGCG‑TCCC‑T‑TTTTCCTGTACCGGGA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg < 1:634064/146‑1 (MQ=255) tgtgTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa > 1:415923/1‑150 (MQ=255) tttttttcccgggCCGGAAACTGACCCGGGGCGCGTT‑GGTTGTGCCTCTCGGCCCCCG‑GGGCG‑TCCC‑T‑TTTCCTTGTTTCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa < 1:1086262/138‑1 (MQ=255) ctcTGTGCGGG‑AAC‑GAACCGGGGGGCGGTGGGGTGTGCCTC‑CGGGCCCCT‑GGGGGTTCCC‑T‑TTTACCTGTATGGGAA‑AAAG‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata < 2:189651/148‑1 (MQ=255) cACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat > 2:788386/1‑151 (MQ=255) tgCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc > 1:380946/1‑150 (MQ=255) gaaaCTGACCCTTGGCGGGTT‑GGTTGTGTCTCCGGAACCCCT‑GGGCG‑T‑CC‑C‑TTTTCCTTTTCCGGGG‑AAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt < 1:396090/147‑1 (MQ=255) accggcccggTGGCGCGTT‑GGATGTGGCTCCCGGCCACCT‑GGGCG‑TACC‑T‑TTTCCCTGTATGGGGA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt < 2:262218/140‑1 (MQ=255) aaCTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc > 1:468079/1‑151 (MQ=255) tGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc > 2:200081/1‑151 (MQ=255) gAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg > 2:261005/1‑151 (MQ=255) cgtgcggcgcT‑GGAGGTGCTCCTCGTACCCTG‑GGGCG‑TACT‑T‑TTTCCTTTTGTCGGAA‑AAAT‑AACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg < 2:1222531/141‑1 (MQ=255) ggCGGGTT‑GGGTGTGCCTCTCGTCCACCT‑GGGTT‑TCCC‑G‑TTTCCTTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc < 1:703850/143‑1 (MQ=255) cgcT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGTCGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg > 2:322051/1‑151 (MQ=255) ggt‑GGATGTGTCTCTCGGCCCCTG‑GGGCG‑TCCT‑T‑TTTCCTTGTAGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg < 2:243903/148‑1 (MQ=255) gcT‑GGATGTGGCTATCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg > 2:651959/1‑150 (MQ=255) t‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg > 1:1188573/1‑151 (MQ=255) gttggtgcctccGGACCCCTGGGGCG‑TCCC‑T‑TTTACCTTTA‑CGGGCAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc < 2:95360/138‑1 (MQ=255) tgGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa > 2:937685/1‑150 (MQ=255) ggCTCTGGGCCCCCG‑GGGCG‑TCCC‑T‑TTTCCTTTTAGGGGGA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa < 1:284147/150‑1 (MQ=255) ctcGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc > 1:851242/1‑151 (MQ=255) cgggggCG‑TCCC‑T‑TTTCCTTGTAGGGGAA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa < 1:1132430/147‑1 (MQ=255) tttCCCTGTATCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg < 1:920961/148‑1 (MQ=255) ttACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc > 1:433743/1‑151 (MQ=255) cTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt > 1:1273077/1‑151 (MQ=255) tATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt > 1:637848/1‑151 (MQ=255) cGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc > 2:1178118/1‑151 (MQ=255) gggA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt < 1:140276/148‑1 (MQ=255) | CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC > NC_000913/436366‑436662 |
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: ‑ |