Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 955,156 | T→C | 16.0% | intergenic (+107/‑206) | 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,156 | 0 | T | C | 16.0% | 107.7 / 5.6 | 47 | intergenic (+107/‑206) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base T (20/17); new base C (4/2); total (26/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCT > NC_002947/955016‑955290 | cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGa < 1:645400/150‑1 (MQ=255) cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCggg < 1:229824/150‑1 (MQ=37) ggATGTGCTGATGACTTACGCACGCTGAAACTGCTGCCTGTTTCACGAACACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt < 1:981969/150‑1 (MQ=255) aTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTa < 2:2213358/139‑1 (MQ=255) aTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTa > 1:2213358/1‑139 (MQ=255) tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCcgc > 2:1157984/1‑150 (MQ=37) tgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCGGGTTTACCCGCGaa < 1:1157984/150‑1 (MQ=25) cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGt > 1:1546805/1‑125 (MQ=255) cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGt < 2:1546805/125‑1 (MQ=255) tACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg > 2:394219/1‑132 (MQ=255) tACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 1:394219/132‑1 (MQ=255) gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCgg > 2:2141530/1‑129 (MQ=25) cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt > 2:718947/1‑127 (MQ=255) cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt < 1:718947/127‑1 (MQ=255) cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAg > 1:463698/1‑150 (MQ=255) tgctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCggg > 2:1343629/1‑118 (MQ=25) tgctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 1:1343629/118‑1 (MQ=25) tgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaa < 2:463698/150‑1 (MQ=255) ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 1:1665262/112‑1 (MQ=255) ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg > 2:1665262/1‑112 (MQ=255) cGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTATTACTGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGGTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAttt > 1:1918220/1‑150 (MQ=38) cccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGc > 2:629269/1‑114 (MQ=255) cccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGc < 1:629269/114‑1 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑Tg > 2:2271086/1‑86 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑Tg < 1:2271086/86‑1 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAaga < 1:745957/110‑1 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAaga > 2:745957/1‑110 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa > 1:327087/1‑150 (MQ=255) gTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa < 2:327087/150‑1 (MQ=255) gtgCCGGCTTCCTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc < 2:1918220/150‑1 (MQ=255) gCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACaa > 2:108243/1‑149 (MQ=255) ttCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt < 1:1084614/102‑1 (MQ=255) ttCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAAGAGAt > 2:1084614/1‑102 (MQ=255) gggTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc < 1:1634651/134‑1 (MQ=255) gggTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc > 2:1634651/1‑134 (MQ=255) cccGCTCCTAC‑AGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc > 2:1658615/1‑150 (MQ=37) cccGCTCCTAC‑AGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc > 1:1473133/1‑150 (MQ=37) ccGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt > 2:1625279/1‑117 (MQ=25) ccGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt < 1:1625279/117‑1 (MQ=25) ccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 2:1987084/1‑150 (MQ=255) ccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 1:497872/1‑150 (MQ=255) ccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 2:1179136/1‑150 (MQ=255) ccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 1:1789683/1‑150 (MQ=255) aCAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc < 2:1172793/92‑1 (MQ=255) aCAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc > 1:1172793/1‑92 (MQ=255) aaGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAg > 1:876665/1‑118 (MQ=255) aaGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAg < 2:876665/118‑1 (MQ=255) | CGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCT > NC_002947/955016‑955290 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |