Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,697,316 | G→C | 8.6% | intergenic (+26/‑51) | hslV → / → hslU | ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,697,316 | 0 | G | C | 8.6% | 176.3 / 6.9 | 58 | intergenic (+26/‑51) | hslV/hslU | ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit |
Reads supporting (aligned to +/- strand): ref base G (31/22); new base C (2/3); total (33/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.43e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01 |
CAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGC > NC_002947/5697176‑5697457 | cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAggg > 5:155229/1‑149 (MQ=255) aaGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAgggg > 1:307217/1‑149 (MQ=255) gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAgg < 1:64613/145‑1 (MQ=255) gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAgg > 2:64613/1‑145 (MQ=255) gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAGGGGcg < 5:341517/149‑2 (MQ=255) gaccgaccTCTCGCCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCt < 4:237966/149‑1 (MQ=255) accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCtt > 1:251877/1‑149 (MQ=255) ccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAGGGGcgcc < 6:241427/149‑4 (MQ=255) cgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAGGGGcgcca < 2:307217/149‑5 (MQ=255) ctctCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCAc < 6:155229/149‑1 (MQ=255) tctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACg > 8:324893/1‑149 (MQ=255) ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa < 3:111240/149‑1 (MQ=255) ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa < 4:64512/149‑1 (MQ=255) gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGccc > 5:362407/1‑147 (MQ=255) ccGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAg > 2:198964/1‑149 (MQ=255) agaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGAc < 4:88945/149‑1 (MQ=255) gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa > 8:108103/1‑148 (MQ=255) gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa < 7:108103/148‑1 (MQ=255) tGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAATGACCGTATTGAt < 2:251877/149‑1 (MQ=255) tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc > 7:265002/1‑149 (MQ=255) cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCt < 3:299064/111‑1 (MQ=255) cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCt > 4:299064/1‑111 (MQ=255) cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATg < 6:362407/149‑1 (MQ=255) gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc > 5:202406/1‑143 (MQ=255) cACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc > 3:335090/1‑149 (MQ=255) cACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc > 5:158146/1‑149 (MQ=255) accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAggg < 5:62550/79‑1 (MQ=255) accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAggg > 6:62550/1‑79 (MQ=255) accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTcc > 6:224944/1‑149 (MQ=255) caaccaCACACTGACCATCGAGGATCAGGACCTGGCAGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCGCATGGACCGTATTTATCATGTCCTTGTCCCCCAGCGAGTTCGTCCAc < 8:66168/149‑1 (MQ=255) aaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc > 8:20613/1‑148 (MQ=255) cacaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGaa > 8:223612/1‑148 (MQ=255) cacaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAAGACCCCCCGCGAGATCGTCCACGAAc < 1:198964/149‑1 (MQ=255) ccTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa > 6:350653/1‑83 (MQ=255) ccTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa < 5:350653/83‑1 (MQ=255) cTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGTTCGTGGACGAtctcgcgg < 7:324893/149‑5 (MQ=255) aCCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTCGTTATGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGCCTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGc > 1:65549/1‑149 (MQ=255) ccATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCca < 4:335090/149‑1 (MQ=255) aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCcaca < 7:223612/149‑1 (MQ=255) aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGACATCGTCCACGAACTCACCCGCcac > 4:7013/1‑148 (MQ=255) cGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCGAcat > 6:126035/1‑148 (MQ=255) aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat < 6:158146/149‑1 (MQ=255) aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat > 2:361662/1‑149 (MQ=255) aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat > 6:89427/1‑149 (MQ=255) aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat > 1:56860/1‑149 (MQ=255) gCAGGACCTTGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTc < 8:208745/131‑1 (MQ=255) gCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTc > 7:208745/1‑131 (MQ=255) cTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGgac > 2:85045/1‑148 (MQ=255) tGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGCTCATGTCCATGACCCCCCGCGAGATCGTCCCCGAACTCAACCGCCACATCATCGGCCAGgacgc > 8:282972/1‑148 (MQ=255) cGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGc > 7:257372/1‑148 (MQ=255) gATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcg > 1:2090/1‑149 (MQ=255) gATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcg > 5:319365/1‑149 (MQ=255) aTTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGACCTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcgc > 4:297923/1‑149 (MQ=255) cT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg < 1:361662/149‑1 (MQ=255) tGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAtt > 7:262449/1‑149 (MQ=255) aaaagcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAtt > 1:351163/5‑149 (MQ=255) cgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc < 7:20613/149‑1 (MQ=255) ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc < 8:262449/149‑1 (MQ=255) | CAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGC > NC_002947/5697176‑5697457 |
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: ‑ |