| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 955,155 | T→A | 14.6% | 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,155 | 0 | T | A | 14.6% | 64.8 / 5.2 | 28 | intergenic (+106/‑207) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
| Reads supporting (aligned to +/- strand): ref base T (14/9); new base A (2/2); total (17/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01 | |||||||||||
CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC > NC_002947/955010‑955302 | ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑t < 1:392197/150‑1 (MQ=255) cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑Tg > 1:147638/1‑150 (MQ=255) aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑Tg > 2:1666860/1‑150 (MQ=255) cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGc < 2:1591653/150‑1 (MQ=255) cGGATGGGCTGATTACCTACGCACGCTCAACCTGCTGCCTGAATCTTTAAAATACCTTTCTCTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCggg < 2:50214/150‑1 (MQ=25) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 1:428205/140‑1 (MQ=37) gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg > 2:428205/1‑140 (MQ=37) gaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGaa < 1:1666860/150‑1 (MQ=255) aCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTCTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTCCGCCCTTGAACCCGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCCGCGAaga < 1:987763/150‑1 (MQ=25) aCCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTCAACCCGCTTCTACAAGTTTCC‑TGTAGGAGCGGGTTTGCCCGCGCAGCGGCCAGCACAGACAAATAAAAATACCTGATTTTCaac > 2:641235/1‑149 (MQ=38) ccACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGACAAATTAAAATACCTGATTTTCACAAATTAGa > 1:397471/1‑150 (MQ=255) acagtggCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt < 2:619063/144‑1 (MQ=255) ttCGCGGGAACACCCGCTCCTACACGGGTCA‑TGTCGGAGCGGGATTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGt > 2:1261822/1‑150 (MQ=255) gggTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTATAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg > 2:392029/1‑150 (MQ=255) aCCCGCTCCTACAAGGTTCC‑TGAAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt > 2:1791325/1‑150 (MQ=255) ccGCTCCTACAGGAACCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATATGCAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCGTCCTGTTTCGGCACCAtct > 1:315690/1‑150 (MQ=25) gCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaaa < 2:810934/64‑1 (MQ=255) gCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaaa > 1:810934/1‑64 (MQ=255) gCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc > 2:200221/1‑150 (MQ=255) cTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttct > 2:1585673/1‑150 (MQ=255) tCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg > 1:1697658/1‑150 (MQ=255) tCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg > 1:1437482/1‑150 (MQ=255) tACAAGGTTCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt > 1:259003/1‑150 (MQ=255) ggTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc < 2:1697658/150‑1 (MQ=255) gTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa > 1:627730/1‑90 (MQ=255) gTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa < 2:627730/90‑1 (MQ=255) gTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt < 1:552643/133‑1 (MQ=255) gTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt > 2:552643/1‑133 (MQ=255) | CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC > NC_002947/955010‑955302 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |