Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 955,083 | T→A | 9.8% | intergenic (+34/‑279) | 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,083 | 0 | T | A | 9.8% | 137.4 / 5.7 | 51 | intergenic (+34/‑279) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base T (17/28); new base A (2/3); total (20/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 |
GGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTT > NC_002947/954941‑955217 | ggCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCgg < 2:167768/150‑1 (MQ=255) ggCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCgg < 1:955948/150‑1 (MQ=255) gCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGCTGACCTACGCCAGCTGAAGATGCTGCCGCGGTCCCGTAAACGCCGGCCACGGT‑GCCggg > 2:589985/1‑150 (MQ=255) cAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGtttt < 1:1368470/150‑1 (MQ=255) cAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGCGCATATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGgttt > 1:101697/1‑150 (MQ=255) tCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTtgg < 2:739887/150‑1 (MQ=255) ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTAt > 1:1591482/1‑150 (MQ=255) aTCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTcc < 1:1341063/150‑1 (MQ=255) tCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCCGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCt < 1:589985/150‑1 (MQ=255) tCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCt < 2:1591482/150‑1 (MQ=255) tCACCTTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCtgtg > 2:541600/1‑150 (MQ=255) cGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttttc < 1:1345648/124‑2 (MQ=255) cGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttttc > 2:1345648/1‑123 (MQ=255) cACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCATGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCct < 1:344468/150‑1 (MQ=255) ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc < 1:1121256/150‑1 (MQ=255) ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc < 1:1494138/150‑1 (MQ=255) ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc > 1:501372/1‑150 (MQ=255) ccTGAGCGTGATGACGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGATGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTGTTGGGGGTATTTGTACTGTTCCGGCCTCTAcgc > 1:625610/1‑150 (MQ=255) cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgcg < 1:1194916/150‑1 (MQ=255) tGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg < 1:966503/150‑1 (MQ=255) gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAACCCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc < 2:1210681/150‑1 (MQ=255) ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg < 1:541600/150‑1 (MQ=255) gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGTTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTGTTGTTTTGGGGGTATTTTTCCTGTGCCGGACTCTTCGCGGGTAAACCAGCTCCTGCagga > 1:786089/1‑147 (MQ=255) aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc < 2:764693/150‑1 (MQ=255) aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc < 1:2200/137‑1 (MQ=255) aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc > 2:2200/1‑137 (MQ=255) cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGc < 2:501372/150‑1 (MQ=255) gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGTTTGTTTTACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGTCTCTTCGCGGGTAAACCCGCTTCTACAAGGTTCC‑TGTAGGAGCgg < 2:625610/150‑1 (MQ=255) gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACAAAACCCCCGGCTATTGT‑GCCGGGGTTTTTTTTTTGTGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 2:183542/150‑1 (MQ=39) cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCggg < 2:101697/150‑1 (MQ=255) cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCggg < 1:304201/150‑1 (MQ=37) cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCggg > 1:859708/1‑63 (MQ=255) cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCggg < 2:859708/63‑1 (MQ=255) atgaACTACACACGCTGAAACTGCTGCCTGTTTCACGAAAACACCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGa > 2:61985/1‑150 (MQ=255) gaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCgg < 1:1186107/136‑1 (MQ=37) gaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AGGAACCTTGTAGGAGCgg > 2:1186107/1‑136 (MQ=37) gAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTgggg < 2:1421922/58‑1 (MQ=255) gAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTgggg > 1:1421922/1‑58 (MQ=255) aaCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTc > 1:84085/1‑65 (MQ=255) aaCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTc < 2:84085/65‑1 (MQ=255) gcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 1:1526246/113‑1 (MQ=255) gcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg > 2:1526246/1‑113 (MQ=255) tCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg > 2:782703/1‑106 (MQ=255) tCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCgg < 1:782703/106‑1 (MQ=255) gAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGCtttt > 2:736484/1‑150 (MQ=255) aaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTNTGCCGGCCTCTTCGCGGGTAAACCCGCTc > 1:66055/1‑77 (MQ=25) aaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc < 2:66055/77‑1 (MQ=25) ccccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt > 2:287849/1‑99 (MQ=255) ccccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGt < 1:287849/99‑1 (MQ=255) ccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc < 1:879070/70‑1 (MQ=255) ccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc > 2:879070/1‑70 (MQ=255) | GGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCC‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTT > NC_002947/954941‑955217 |
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: ‑ |