Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,088,662 | C→A | 6.7% | intergenic (+158/‑96) | pmbA → / → PP_0943 | regulatory protease/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,088,662 | 0 | C | A | 6.7% | 186.9 / 3.2 | 60 | intergenic (+158/‑96) | pmbA/PP_0943 | regulatory protease/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (27/29); new base A (2/2); total (29/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 = 5.44e-01 |
TGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGC > NC_002947/1088514‑1088802 | tGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGGTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTc > 5:195122/1‑149 (MQ=255) gCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCgg < 7:274323/149‑1 (MQ=255) cgGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAAACCCGGCGCGAGGGCCGGGtcttcgga < 1:145359/149‑5 (MQ=37) aaaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCAGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATACGAAGACCCGGACCTCGCGCCGGGTTTTTATTtcat > 8:16987/1‑148 (MQ=255) aaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTtcatca < 2:215686/149‑1 (MQ=255) aaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTtcatc < 2:115295/148‑1 (MQ=255) ccGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTGAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCAcc > 1:241616/1‑149 (MQ=255) ccTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGtttt > 2:307416/1‑131 (MQ=255) ccTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGtttt < 1:307416/131‑1 (MQ=255) cTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTACAAATCCGAAGACCCGGCGCGAGGGCCGGGtcttc > 3:146020/1‑130 (MQ=25) cTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTACAAATCCGAAGACCCGGCGCGAGGGCCGGGtcttc < 4:146020/131‑2 (MQ=25) acacGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGttt > 1:54478/1‑127 (MQ=255) acacGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGttt < 2:54478/127‑1 (MQ=255) cacGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTg < 2:317250/148‑1 (MQ=255) ggtggtGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGt > 3:174886/1‑149 (MQ=255) gACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTa > 7:237510/1‑149 (MQ=255) gggTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTta < 2:216131/149‑1 (MQ=255) ggTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTAATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaa > 5:206698/1‑149 (MQ=255) acacAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATg > 7:23430/1‑148 (MQ=255) ggCCGTTGT‑AGGTGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATc < 2:256841/149‑1 (MQ=255) gTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTa > 1:189889/1‑149 (MQ=255) tcttcgCGGGTGAACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATc > 7:152396/6‑147 (MQ=255) gAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACt > 7:93771/1‑124 (MQ=255) gAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACt < 8:93771/124‑1 (MQ=255) gggTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAAccc < 2:241616/149‑1 (MQ=255) cccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCTGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatg < 7:16987/149‑1 (MQ=255) ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaa < 5:16625/115‑1 (MQ=255) ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaa > 6:16625/1‑115 (MQ=255) ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatga > 7:64783/1‑149 (MQ=255) ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatga > 5:130107/1‑149 (MQ=255) cgAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGt < 4:174886/149‑1 (MQ=255) gCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGa > 1:243432/1‑88 (MQ=255) gCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGa < 2:243432/88‑1 (MQ=255) ccGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGaaa < 8:152396/149‑1 (MQ=255) cGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAt > 3:310572/1‑149 (MQ=255) cgAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCCCGCATGGTTTTTATTTCATCACCTTTTTCTTGAAGTCATTCTATTTATCAATTAATAATGAATATAATTATCCAGTACCCCGACGATGATGTTCATGAAATc > 4:62826/2‑149 (MQ=255) gCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGtcct < 6:195122/149‑1 (MQ=255) cACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTg < 1:20373/72‑1 (MQ=255) cACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTg > 2:20373/1‑72 (MQ=255) cACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGtcctc < 8:64783/149‑1 (MQ=255) tACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatg > 2:250814/1‑125 (MQ=255) tACACCTTTCCCTGCAAATCCGAAGACCCGGCCATCGCGCCCGGTTTTTATTTCAGCACCTTTTTCTTGAAGTTATTCTATTAAACACTTAATAAAAAATAACATTATCCATTAAGCCGTCGATCATGTTCATTAAAACCGTCCTCcgc > 8:188565/1‑149 (MQ=255) cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGTTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACt < 8:237510/149‑1 (MQ=255) cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCTCTTAATAATGAATATCATTATCCAGTAACCCgg < 2:270852/115‑1 (MQ=255) cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCTCTTAATAATGAATATCATTATCCAGTAACCCgg > 1:270852/1‑115 (MQ=255) cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACt > 6:201814/1‑149 (MQ=255) gCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcgcg > 5:263797/1‑137 (MQ=255) gCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcgcg < 6:263797/137‑1 (MQ=255) cAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTAc < 4:201072/147‑1 (MQ=255) cAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTAc < 5:201814/147‑1 (MQ=255) aaaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTg < 2:189889/149‑1 (MQ=255) aaaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTg < 6:130107/149‑1 (MQ=255) aaaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGATCATGAAATCCGTCCTCGGCGAACTGCCCTCCCTg > 8:123986/1‑149 (MQ=255) aaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTgg < 7:67547/149‑1 (MQ=255) tCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGaa < 4:310572/149‑1 (MQ=255) ataaaaaCCCGGCGCGAGGGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGaa > 8:267758/7‑148 (MQ=255) gACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTgg > 4:257537/1‑149 (MQ=255) gACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTgg < 3:62826/149‑1 (MQ=255) aCCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGgc < 1:320202/149‑1 (MQ=255) cccGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGgcg > 6:263726/1‑149 (MQ=255) ccGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCATCCGCGAACTGCCCTACCTGGAAAACTGgcgc > 2:237981/1‑149 (MQ=255) | TGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGC > NC_002947/1088514‑1088802 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |