Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,655,636 | C→A | 25.0% | 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,636 | 0 | C | A | 25.0% | 43.8 / 11.2 | 24 | intergenic (‑91/+61) | epd/tktA | D‑erythrose 4‑phosphate dehydrogenase/transketolase A |
Reads supporting (aligned to +/- strand): ref base C (10/8); new base A (3/3); total (13/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01 |
CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGA > NC_002947/5655488‑5655783 | cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGc < 7:123012/149‑1 (MQ=255) acaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAACCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGc > 1:2918/1‑148 (MQ=255) cGATACGGCCATAACCGTTGAGTGCATCTTTGTAGGGACGGTGGTGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAT‑GGTTTCATGCAACCTCtgt < 6:63021/149‑1 (MQ=255) tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCT‑TCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGa < 5:270134/149‑1 (MQ=17) ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGAAACC‑CTGTGGGAGCg < 8:83815/149‑1 (MQ=255) ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg < 6:155845/149‑1 (MQ=255) aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc < 7:38000/149‑1 (MQ=255) aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc < 4:286001/149‑1 (MQ=255) aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTg < 1:171563/149‑1 (MQ=255) cTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAACCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAgg < 2:2918/149‑1 (MQ=255) cccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGaa > 1:62800/1‑149 (MQ=255) gCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACaa > 2:218479/1‑149 (MQ=38) gCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACaa > 2:229613/1‑149 (MQ=255) tCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATg > 4:62678/1‑149 (MQ=255) ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGc > 6:8804/1‑149 (MQ=17) cacaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGcc > 4:68014/1‑148 (MQ=255) acagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAAt < 1:301703/139‑1 (MQ=38) acagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAAt > 2:301703/6‑144 (MQ=38) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCATCCAGGGTGAAGCCAAACTCTTCGAACAATGc < 1:41098/147‑1 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGc > 2:41098/1‑146 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 3:209971/1‑148 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGAAGATC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 2:330004/1‑148 (MQ=255) tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAgg < 2:323220/114‑1 (MQ=255) tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAgg > 1:323220/1‑114 (MQ=255) tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCgg > 5:113280/1‑149 (MQ=255) gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCAGa > 5:263417/1‑149 (MQ=255) | CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTC‑TCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGA > NC_002947/5655488‑5655783 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |