Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 955,152 | T→A | 31.9% | intergenic (+106/‑207) | 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,152 | 0 | T | A | 31.9% | 35.2 / 16.7 | 28 | intergenic (+106/‑207) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base T (8/11); new base A (5/4); total (13/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01 |
GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT > NC_002947/955006‑955294 | gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 8:166089/149‑1 (MQ=255) gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 3:304137/149‑1 (MQ=255) gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 4:274008/149‑1 (MQ=255) ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc > 8:225850/1‑149 (MQ=255) ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc < 3:84199/149‑1 (MQ=255) cAGACGGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGTTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt < 4:61909/149‑1 (MQ=255) cGCTGCGGATGTGCTGATGACCTACGCACGCTGAATACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCGTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc > 1:65341/1‑143 (MQ=255) cGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 2:65341/142‑1 (MQ=255) gCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGTCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa < 4:64027/149‑1 (MQ=255) cTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCACGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAg > 6:263782/1‑149 (MQ=38) cGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGTCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg < 5:263782/149‑1 (MQ=37) cGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg < 4:372088/149‑1 (MQ=255) ggATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg < 3:23176/149‑1 (MQ=255) tgCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc < 7:225850/149‑1 (MQ=255) tgCTGATGACCTACGCACGCTGAA‑ACTGCGGCCTGTTTCACGAAAACCCCGGCCACTGTACCGGGGTTTTTTTTTGGGGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACGAACCTTGTAGGAGCgg < 1:291053/148‑1 (MQ=17) gTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg > 4:407768/1‑109 (MQ=21) gTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg < 3:407768/109‑1 (MQ=21) tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc > 5:310549/1‑112 (MQ=255) tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc < 6:310549/112‑1 (MQ=255) gTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTGTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCGGATGTGCAAAAAATAGATCTGAAAAGTTCTACTAAAGCAGGAAATTATCATGTACAGGCATCCTggtgggg > 4:319969/1‑142 (MQ=21) aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa > 4:26981/1‑149 (MQ=255) cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt > 3:228560/1‑149 (MQ=255) ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc > 2:85810/1‑149 (MQ=255) gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa > 5:164009/1‑125 (MQ=25) gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa < 6:164009/125‑1 (MQ=25) gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 4:175827/1‑148 (MQ=38) tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt > 2:354567/1‑149 (MQ=255) cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt > 1:17746/1‑148 (MQ=255) | GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT > NC_002947/955006‑955294 |
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: ‑ |