Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 955,151 | G→A | 17.2% | intergenic (+105/‑208) | cyoE → / → alaC | protoheme IX farnesyltransferase/aminotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 955,151 | 0 | G | A | 17.2% | 106.6 / 16.1 | 57 | intergenic (+105/‑208) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base G (22/20); new base A (4/5); total (26/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-02 |
AGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG > NC_002947/955009‑955297 | aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTg < 1:44815/149‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGt < 4:18114/149‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGt < 4:293625/149‑1 (MQ=255) cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCgg < 3:52502/149‑1 (MQ=38) tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAg > 3:94/1‑137 (MQ=255) tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAg < 4:94/137‑1 (MQ=255) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTAAACCCC‑TACTTAAAGGTTGCTGTATGGGCGGGTTTTCCcgc > 3:110820/1‑149 (MQ=255) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCC‑TCCTACAAGGTTACTGTAGGAGCGGGTTTACCcgc > 4:70236/1‑149 (MQ=255) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCC‑TCCGACAACGTTCCTGTAGGAGCGGGGTTTCCcgc > 1:8926/1‑149 (MQ=39) aCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgcg > 4:244503/1‑146 (MQ=255) ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg < 1:60231/135‑1 (MQ=37) ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg > 2:60231/1‑135 (MQ=37) cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCGGGt > 7:45013/1‑140 (MQ=17) cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCGGGt < 8:45013/140‑1 (MQ=17) gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc > 8:225635/1‑147 (MQ=255) cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCC‑ACCTACAAGGTTACTGTTTGAGCGGGTTTTACCGCGAACAGGCCAg > 1:293710/1‑145 (MQ=255) aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCC‑TCCAAAAAGGTTGCGGTATGAGGGGGGTTAACCGCGAATAGGCCAGCACAGaa > 3:109057/1‑148 (MQ=255) ttCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg > 4:69233/1‑108 (MQ=21) ttCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg < 3:69233/108‑1 (MQ=21) actaaaaCCCCGGCAACTGGGCGGGGT‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGTGTAAACCCGC‑TCCTAAAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAt < 4:110820/146‑1 (MQ=37) cGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGC‑TACTACAAAGTTACTGTACGATCCGGTTTAACCCCGAAGACGCCAGCAAAGAAAAATAAAAAAAACTGaat > 8:66980/1‑147 (MQ=255) cGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa < 5:103376/149‑1 (MQ=255) cGGCCACTCTGCAGGGG‑TTTTTTTTTGGGCGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCCCGGTAAAACCCTC‑TCATACAATTTACCTTTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa < 3:65110/149‑1 (MQ=37) gCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTTCCTACAAGGTCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAg < 1:241658/119‑1 (MQ=25) cACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg > 3:31696/1‑91 (MQ=17) cACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg < 4:31696/91‑1 (MQ=17) aCTGTGCCGGGG‑TTTTTTTTTGGGGTTAT‑‑‑‑‑TTTTCTTGTGCGGGCCTCTTCGCGGTTAAACCCGC‑TCTTACATGTTTCCTTTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt < 2:8926/149‑1 (MQ=38) aCTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt < 3:244503/149‑1 (MQ=255) cTGGGCGGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATc < 2:196010/149‑1 (MQ=255) gttttttttttGGGGTTTT‑‑‑‑‑TTTTCCTGGGCCGGCGGCTGCGCGGTTAACCCCGC‑TCAACCACGGTTCCTTTTGGCGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGtt < 2:280024/147‑1 (MQ=255) tttttttttGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCTGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt < 7:225635/149‑1 (MQ=255) tGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGTAAAACCCGC‑TCCTACAAGTTTCCTTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa < 3:127192/148‑1 (MQ=38) gggtTAT‑‑‑‑‑TTTTCTTGTTCCGGCCTCTTCTCGTTTAAACCCTA‑TATTACAAGTTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAg < 5:116719/149‑1 (MQ=37) gtgCCGGCCTCTTCGCGGTTAAACCCGC‑TCCAACAATGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAtt < 5:64885/99‑1 (MQ=25) tgCCGGCCTCTTCGCGGGTAAACCCGA‑TCCTACAAGGTTCCTGTAGTAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt < 1:252110/112‑1 (MQ=255) ctTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTTCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt > 3:242568/1‑149 (MQ=255) gcgGTTAAACCCGC‑TCCTAAAAGGTTCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGtt < 1:120212/148‑1 (MQ=255) gcgGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt > 1:74669/1‑145 (MQ=255) gcgGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt < 2:74669/145‑1 (MQ=255) ggTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg > 6:131063/1‑149 (MQ=255) aaCCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc > 8:179592/1‑148 (MQ=255) aCCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa > 5:41882/1‑149 (MQ=255) ccGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta > 2:237095/1‑102 (MQ=255) ccGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta < 1:237095/102‑1 (MQ=255) ccGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa < 2:112793/115‑1 (MQ=255) ccGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa > 1:112793/1‑115 (MQ=255) ccGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc > 1:117283/1‑149 (MQ=255) c‑TCCTACAAGTTTCCTTTAGTAGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc < 7:276016/149‑1 (MQ=255) cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa > 2:37588/1‑149 (MQ=255) tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt > 2:230955/1‑149 (MQ=255) ctgaaGTTTCCTGTAGGAGCTGGTTTACACGCGAAGAGGCCAGCACAGATAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa < 8:207243/145‑1 (MQ=255) aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt > 1:54031/1‑148 (MQ=255) cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt < 7:179592/147‑1 (MQ=255) aaGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa < 3:70236/149‑1 (MQ=255) aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAATATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt > 5:230987/1‑147 (MQ=255) aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATATGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTACTGAtag > 7:207243/1‑146 (MQ=255) aGGTTCCTGTAGTTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg < 7:293662/149‑1 (MQ=255) ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc > 4:50159/1‑131 (MQ=255) ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc < 3:50159/131‑1 (MQ=255) | AGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG > NC_002947/955009‑955297 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |