Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 955,076 | Δ1 bp | 6.9% | intergenic (+30/‑283) | 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,076 | 0 | C | . | 6.9% | 182.1 / 7.2 | 58 | intergenic (+30/‑283) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (29/25); new base . (2/2); total (31/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01 |
AGCAAGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAA > NC_002947/954930‑955221 | aGCAAGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCa > 3:137332/1‑149 (MQ=255) aGCAAGTGGGCCCGCAAGGTGTTCGGCTTCTCCATACTAACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCa > 6:316973/1‑149 (MQ=255) gTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACtgt‑g > 2:203671/1‑149 (MQ=255) cccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCggg < 3:258834/149‑1 (MQ=255) ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCgggg < 4:229979/149‑1 (MQ=255) ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCgggg > 6:125861/1‑149 (MQ=255) gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ttttttttt > 3:49050/1‑149 (MQ=255) gtCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGGGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTg < 2:83388/148‑1 (MQ=255) cGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgg < 1:203671/148‑1 (MQ=255) ggCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAACCCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgg < 1:241520/148‑1 (MQ=255) ggCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgggg < 5:252247/149‑1 (MQ=255) cTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCCCGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGt < 4:49050/149‑1 (MQ=255) tCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTcc > 3:217389/1‑149 (MQ=255) cTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCt < 5:125861/148‑1 (MQ=255) tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCtgt > 7:133415/1‑149 (MQ=255) tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCtgt > 4:57649/1‑149 (MQ=255) tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttc > 5:244527/1‑127 (MQ=255) tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttc < 6:244527/128‑2 (MQ=255) aCCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCtgtg < 4:299850/148‑1 (MQ=255) aCCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCtgt < 4:217389/147‑1 (MQ=255) aCTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGcc > 4:232429/1‑148 (MQ=255) aGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg < 3:22064/148‑1 (MQ=255) gCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGt > 6:295750/1‑149 (MQ=255) gTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTa > 3:262873/1‑148 (MQ=255) atgatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAc < 3:232429/148‑1 (MQ=255) gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACt > 7:254927/1‑89 (MQ=255) gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACt < 8:254927/89‑1 (MQ=255) gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTATTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTcc > 1:214514/1‑148 (MQ=255) gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttcg > 2:295228/1‑97 (MQ=38) gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttcg < 1:295228/99‑3 (MQ=38) gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACag > 7:86430/1‑148 (MQ=255) aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg > 7:9436/1‑149 (MQ=255) aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGCATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg > 7:6879/1‑149 (MQ=255) gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 3:57649/149‑1 (MQ=255) cGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACACCAACTACAACGTTccaggag > 4:33797/1‑143 (MQ=255) cGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc < 4:262873/147‑1 (MQ=38) cGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg < 5:251289/149‑1 (MQ=255) ggATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg < 8:86430/149‑1 (MQ=37) cTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc < 8:9436/149‑1 (MQ=255) tgatgaCCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCAACTAGAAGGTTCCGGTAGGCGCGGGTTTAGCcg > 7:168521/1‑149 (MQ=255) ccTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACACCGCCTCCAAGGGTGCGGGAGGAGCGGGGTTAACCGGGaa > 6:224986/1‑147 (MQ=255) gcacgcTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa > 2:111348/1‑122 (MQ=255) gcacgcTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa > 4:260971/1‑141 (MQ=255) gcacgcTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa < 3:260971/141‑1 (MQ=255) acgcTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc < 8:210569/104‑1 (MQ=255) acgcTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc > 7:210569/1‑104 (MQ=255) cTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTAGTTTTCCTGTGCCGGACTCTTCGCGGGTAAACCCGCACCAAACAAGTTTCTGGATGACCCGGTTTACCCCCCGATAGGCCAGcac > 7:203660/1‑148 (MQ=255) tGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc < 1:135871/131‑1 (MQ=25) tGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc > 2:135871/1‑131 (MQ=25) cTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTACTACAAGGTTGCTGTAGGAGCGGGTTTTCCCGCGGAGAGGGCCGCACAGaaaaaaaaaaa > 8:244119/1‑148 (MQ=37) tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc < 5:158257/114‑1 (MQ=18) tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc > 6:158257/1‑114 (MQ=18) tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt < 8:203768/96‑1 (MQ=255) tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt > 7:203768/1‑96 (MQ=255) acACGAAACCCCGGCCCCCGGG‑GCCTGGG‑TTTTTTTTTGGGGGTTTTTTACCTGTGCCGGCCTCTTGGCGGGAAAACCCGCTCCTACACGGTTCCTGTAGTTGGGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCt < 8:168521/147‑1 (MQ=25) gAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACACGCAACTACAAGGTTCCTGTAGGGACGGGTTTAACCGCGAAGAGGCCAGCCCACAAAAAAAAAAATAACTGAatt > 2:217971/1‑149 (MQ=37) gAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGGTTGCTGTTGGAGCGGGGTTTGCCGCGAAGCGGCCCGCACAGAAAAATAAAAATAACTGAttt > 1:17962/1‑149 (MQ=38) cGGCCACTGG‑GCCGGGG‑TTTTTTTTTTGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa < 3:33797/149‑1 (MQ=255) | AGCAAGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAA‑ACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAA > NC_002947/954930‑955221 |
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: ‑ |