| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 955,154 | G→A | 13.3% | 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,154 | 0 | G | A | 13.3% | 149.1 / 11.2 | 60 | intergenic (+105/‑208) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
| Reads supporting (aligned to +/- strand): ref base G (23/29); new base A (6/2); total (29/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 | |||||||||||
CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG > NC_002947/955010‑955300 | caaGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGCGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑t < 1:1524224/148‑1 (MQ=255) aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTTCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGt < 2:2652412/150‑1 (MQ=255) aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGt < 2:50639/150‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTa < 2:1974204/150‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTa < 1:578150/150‑1 (MQ=255) cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTGCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGa < 1:2323748/150‑1 (MQ=255) gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAATCCCGCTCCTTCAAGGTTCC‑TGTAGGAGctg < 2:968545/150‑3 (MQ=255) gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCgg < 1:2452768/150‑1 (MQ=37) gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 1:2415126/150‑1 (MQ=255) gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 2:3017049/150‑1 (MQ=255) cgatctgtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACACGAGCCTAC‑AGGAACCTTGTAGGAGCgg < 1:1950887/145‑1 (MQ=37) tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTACCGTTTTGGTCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCgggagatcg > 2:1950887/1‑144 (MQ=25) gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCggg > 1:1715139/1‑144 (MQ=25) gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 2:1715139/144‑1 (MQ=25) tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTAct > 1:2575051/1‑146 (MQ=255) tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTAcc < 2:2575051/147‑1 (MQ=255) tgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCggg > 2:3061534/1‑138 (MQ=37) tgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 1:3061534/138‑1 (MQ=37) ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTCGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑Tg < 1:1999592/103‑1 (MQ=255) ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTCGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑Tg > 2:1999592/1‑103 (MQ=255) cccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCggg > 2:2514244/1‑97 (MQ=18) cccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 1:2514244/97‑1 (MQ=18) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCggg > 1:2176070/1‑96 (MQ=18) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc > 1:304414/1‑83 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 2:304414/83‑1 (MQ=255) ccGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 2:2176070/96‑1 (MQ=18) cGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑t > 2:2029408/1‑84 (MQ=255) cGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑t < 1:2029408/84‑1 (MQ=255) gggg‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 2:1857754/81‑1 (MQ=255) gggg‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg > 1:1857754/1‑81 (MQ=255) g‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGt < 1:1542697/70‑1 (MQ=255) g‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGt > 2:1542697/1‑70 (MQ=255) attttTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt < 1:1683301/62‑1 (MQ=17) attttTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt > 2:1683301/2‑63 (MQ=17) tggCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaa > 1:1793628/2‑103 (MQ=14) tggCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaa < 2:1793628/102‑1 (MQ=14) ccTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc > 2:1564768/1‑150 (MQ=255) gggTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt > 2:713892/1‑99 (MQ=255) gggTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt < 1:713892/99‑1 (MQ=255) ggTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc > 1:1485596/1‑150 (MQ=255) gTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa < 1:1564768/150‑1 (MQ=255) aaaCCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt > 1:310075/1‑150 (MQ=255) aaCCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGccg < 2:362577/52‑2 (MQ=255) aaCCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGccg > 1:362577/1‑51 (MQ=255) aCCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAAt > 2:266995/1‑87 (MQ=255) aCCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAAt < 1:266995/87‑1 (MQ=255) cccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa < 2:1507180/126‑1 (MQ=255) cccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa > 1:1507180/1‑126 (MQ=255) cccGCTCCTAC‑AGGAACCTTGTAGGAGCGGGTTTACCGGCGAAGAGCTTAACACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc > 1:2364419/1‑150 (MQ=37) cccGCTCCTAC‑AGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc > 1:919271/1‑150 (MQ=37) ccGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 2:842104/1‑150 (MQ=255) gCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc > 1:1966666/1‑98 (MQ=255) gCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc < 2:1966666/98‑1 (MQ=255) cTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttct > 1:1250817/1‑150 (MQ=255) cTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa < 1:2408737/109‑1 (MQ=255) cTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa > 2:2408737/1‑109 (MQ=255) aCAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca > 1:1733172/1‑46 (MQ=255) aCAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca < 2:1733172/46‑1 (MQ=255) aCAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc > 1:2658912/1‑150 (MQ=255) aaGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg > 2:465301/1‑150 (MQ=255) | CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG > NC_002947/955010‑955300 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |