Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,655,687 | C→A | 37.9% | intergenic (‑91/+61) | epd ← / ← tktA | D‑erythrose 4‑phosphate dehydrogenase/transketolase A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,655,687 | 0 | C | A | 37.9% | 72.3 / 54.5 | 58 | intergenic (‑91/+61) | epd/tktA | D‑erythrose 4‑phosphate dehydrogenase/transketolase A |
Reads supporting (aligned to +/- strand): ref base C (21/15); new base A (11/11); total (32/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.94e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01 |
CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCG > NC_002947/5655539‑5655833 | cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCa < 2:1585717/150‑1 (MQ=255) ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccct < 1:2219278/150‑3 (MQ=38) ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt < 1:2592678/150‑5 (MQ=38) gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTgg < 1:27949/150‑1 (MQ=38) gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg < 1:178182/150‑1 (MQ=255) gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg < 2:2571466/150‑1 (MQ=255) gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg > 1:2571466/1‑150 (MQ=255) tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt < 1:2165135/150‑5 (MQ=21) aggCCATACGCGTTTAGTGGAACTTGGTAGGGATGGTGGGGCATTCGGGGCTCACGGGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCTCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGc < 2:1057133/149‑1 (MQ=255) ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGc > 2:2820923/1‑149 (MQ=255) ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGc < 1:2820923/149‑1 (MQ=255) ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCg < 1:1345679/150‑1 (MQ=255) gtcATAACCGTTGATTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGTGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg < 1:397713/148‑1 (MQ=255) gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGTGa < 1:3066216/150‑1 (MQ=21) ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCggg < 2:640690/150‑1 (MQ=37) cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc < 2:634377/150‑1 (MQ=255) aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCtt > 2:612240/1‑136 (MQ=255) aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCtt < 1:612240/136‑1 (MQ=255) tttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCtt > 1:2734897/1‑133 (MQ=25) tttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCtt < 2:2734897/133‑1 (MQ=25) ttCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGa > 2:2859095/1‑105 (MQ=21) ttCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGa < 1:2859095/105‑1 (MQ=21) gCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAg < 1:1362990/136‑1 (MQ=37) gCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAg > 2:1362990/1‑136 (MQ=25) gCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGaa < 2:1189103/54‑1 (MQ=255) gCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGaa > 1:1189103/1‑54 (MQ=255) gCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCAGCGGGTG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa > 1:1615130/1‑150 (MQ=255) ccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCa > 1:1841121/1‑85 (MQ=255) ccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCa < 2:1841121/85‑1 (MQ=255) gCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc > 2:1467219/1‑143 (MQ=37) gCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGagatcg > 1:661209/1‑145 (MQ=37) gCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGACAAACTCTTc < 1:1467219/143‑1 (MQ=37) ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 1:359459/1‑150 (MQ=38) ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 1:929236/1‑150 (MQ=38) ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 1:231799/1‑150 (MQ=38) ccACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 2:1162380/1‑150 (MQ=255) cacagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTGACACAATGCCGa > 2:79054/7‑150 (MQ=37) cacagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACAACGCCCGCAGATATTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACAGAGGATGTTCTCCAGGGTGAAGCCAAACTCGTCGAACAATGCCgg > 2:805133/7‑149 (MQ=37) cacaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa > 1:2040751/1‑150 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt > 2:2703305/1‑45 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt < 1:2703305/45‑1 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa > 1:1229179/1‑150 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa > 2:3052886/1‑150 (MQ=255) aG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc > 2:2278386/1‑150 (MQ=255) g‑ggTTTCATGCAACCTCTGTGGGAGCGGGTTAGCCCGTGAATG‑AAGGTGATACCGCTCGCAGATTTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGACCAATGCCGAAGcc > 2:1274420/1‑150 (MQ=255) ggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCg > 2:1553944/1‑150 (MQ=255) tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc > 2:2482785/1‑150 (MQ=255) gcATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg < 2:1959491/111‑1 (MQ=2) gcATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg > 1:1959491/2‑112 (MQ=9) cATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGcc < 1:2482785/150‑1 (MQ=255) aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcag < 1:2826791/76‑1 (MQ=255) aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcag > 2:2826791/1‑76 (MQ=255) aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg > 2:441571/1‑150 (MQ=255) aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg > 1:3047934/1‑150 (MQ=255) aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg > 2:263371/1‑150 (MQ=255) aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg > 2:1395694/1‑150 (MQ=255) tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAAt < 1:2643077/64‑1 (MQ=255) tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAAt > 2:2643077/1‑64 (MQ=255) | CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCG > NC_002947/5655539‑5655833 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |