Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 955,118 | C→T | 8.6% | intergenic (+72/‑241) | 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,118 | 0 | C | T | 8.6% | 123.5 / 2.6 | 47 | intergenic (+72/‑241) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (16/26); new base T (2/2); total (18/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.39e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01 |
CGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAG > NC_002947/954971‑955248 | cGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCg < 1:86193/149‑1 (MQ=255) gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCgg > 8:179389/1‑149 (MQ=255) atgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTTACCCGc < 2:209114/149‑1 (MQ=255) atgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCg < 7:179389/148‑1 (MQ=255) gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGC‑TCCTACAAg > 6:103376/1‑149 (MQ=255) gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCCGCCTCTTCGCGGGTAAACCCGC‑TCCTACaa < 1:26638/148‑1 (MQ=255) ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGTGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTcga > 5:324177/1‑122 (MQ=255) ccAGACCGATGCTGATGTGCTGATGACCTACGCACGCTTAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTcgc < 6:324177/123‑1 (MQ=255) aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTg < 1:44815/149‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGt < 4:293625/149‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGt < 4:18114/149‑1 (MQ=255) cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCgg < 3:52502/149‑1 (MQ=38) tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAg < 4:94/137‑1 (MQ=255) tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAg > 3:94/1‑137 (MQ=255) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTAAACCCC‑TACTTAAAGGTTGCTGTATGGGCGGGTTTTCCcgc > 3:110820/1‑149 (MQ=255) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑Tc < 4:265724/116‑1 (MQ=255) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑Tc > 3:265724/1‑116 (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 > 2:60231/1‑135 (MQ=37) ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg < 1:60231/135‑1 (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 < 3:69233/108‑1 (MQ=21) ttCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg > 4:69233/1‑108 (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 < 4:31696/91‑1 (MQ=17) cACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAGGAACCTTGTAGGAGCggg > 3:31696/1‑91 (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) ggg‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACa > 7:208675/1‑62 (MQ=255) ggg‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACa < 8:208675/62‑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) | CGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGC‑TCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAG > NC_002947/954971‑955248 |
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: ‑ |