Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,266,830 | C→A | 8.0% | D22E (GAC→GAA) | PP_4642 → | cyclic diguanosine monophosphate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,266,830 | 0 | C | A | 8.0% | 182.0 / 4.7 | 61 | D22E (GAC→GAA) | PP_4642 | cyclic diguanosine monophosphate‑binding protein |
Reads supporting (aligned to +/- strand): ref base C (28/28); new base A (2/3); total (30/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.27e-01 |
TTGCTTCGTTCTAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGA > NC_002947/5266684‑5266976 | ttGCTTCGTTCTAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACtt > 7:17150/1‑149 (MQ=255) ttCGTTCTAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACAAGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAt > 8:271653/1‑149 (MQ=255) cGTTCTAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAt < 7:7038/147‑1 (MQ=255) cGTTCTAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACAAGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAt < 7:271653/147‑1 (MQ=255) tAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCa > 1:309919/1‑149 (MQ=255) agagCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCAc > 7:98406/1‑149 (MQ=255) agCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGAAGTCGAT‑Gcgttgaaa < 4:145289/148‑8 (MQ=255) tGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTgc > 1:15477/1‑149 (MQ=255) tGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTgc > 2:133820/1‑149 (MQ=255) gCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACtt < 6:113874/128‑1 (MQ=255) gCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACtt > 5:113874/1‑128 (MQ=255) cGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTgcgc < 3:130669/149‑1 (MQ=255) ggCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCCGATGACCTATTCGAGGATCCATGAAATGATCGATC‑ATGACGAACGCAGCCGATTTCAACGCATCGACTTCGAT‑GCCCACAACGAGCTgcgc > 6:326796/1‑148 (MQ=255) gCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAg > 3:200816/1‑141 (MQ=255) gCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAg < 4:200816/141‑1 (MQ=255) gCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTg > 5:72002/1‑148 (MQ=255) cAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGa < 7:195505/148‑1 (MQ=255) cTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTg < 1:133820/149‑1 (MQ=255) aTGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGc < 8:181508/148‑1 (MQ=255) cacaTGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTg < 7:278940/148‑1 (MQ=255) acaTGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATA‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGc > 8:164432/1‑148 (MQ=255) aTGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCaa > 2:139737/1‑148 (MQ=255) aaCTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATCGATAACGAACGCCCCCGATTTAAACGCTACGACTACTAT‑CCCCCACCCTACATCCGCCAGGGGATGAAACCATGGCCTGAAAAAct > 1:253205/1‑149 (MQ=255) ctctAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAActgct < 1:183078/149‑1 (MQ=255) tctAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAActgctg > 6:181301/1‑149 (MQ=255) cTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCt > 8:49979/1‑149 (MQ=255) aGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAg > 7:40046/1‑101 (MQ=255) aGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAg < 8:40046/101‑1 (MQ=255) tCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGATGGCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAAc < 1:224294/134‑1 (MQ=255) gATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCaaa < 8:98406/149‑1 (MQ=255) gATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCaaa < 2:309919/149‑1 (MQ=255) ccTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTg > 6:278255/1‑140 (MQ=255) ccTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTg < 5:278255/140‑1 (MQ=255) cTATTCGAGTATCCATGCCATGAGCGATC‑ATGACGAACTCCGCCGCTTTCAACGCATCGACTTCGAT‑GGCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTTCTGGACCTTTCACTCAAAGGCTTTCTGATCGAAc < 5:326796/149‑1 (MQ=255) aTTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGc > 1:14815/1‑149 (MQ=255) ttCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCt < 8:133202/90‑1 (MQ=255) gAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAgggg > 8:292230/1‑85 (MQ=255) gAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAgggg < 7:292230/85‑1 (MQ=255) gAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCgg < 1:139737/149‑1 (MQ=255) gAGCATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTATCACTCAAAGGCTTGCTGATCGAACTCCCgg > 6:90249/1‑149 (MQ=255) catgccatgAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTg < 2:15477/149‑1 (MQ=255) atgccatgAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTgg < 7:164432/149‑1 (MQ=255) atgccatgAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTCAAACGCTGGCCTGCAAAACTGCTGGACCATGCACTCAACGGCTTGCTGATCGAAAGCCCCGCACACTgc > 8:292312/1‑148 (MQ=255) aGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACCCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGa < 5:285565/149‑1 (MQ=255) aTC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCt < 5:4314/148‑1 (MQ=255) aTC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCt > 6:4314/1‑148 (MQ=255) gAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGccc < 7:134873/119‑1 (MQ=255) gAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGccc > 8:134873/1‑119 (MQ=255) gAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGa > 5:66225/1‑149 (MQ=255) cgccgccgcTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGAttt < 8:17150/149‑1 (MQ=255) gccgcTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGACTGCAAAACTGCTGGACCTTTCACTCAATGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGAc < 5:90249/149‑1 (MQ=255) cAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCcatca < 2:14815/148‑1 (MQ=255) aaCGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCatcatc > 8:204429/1‑149 (MQ=255) gggCATCGAAGTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGAttt > 5:285338/3‑135 (MQ=255) gggCATCGAAGTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGAttt < 6:308123/133‑1 (MQ=255) gggCATCGAAGTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGAttt < 6:285338/133‑1 (MQ=255) gggCATCGAAGTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGAttt > 5:308123/3‑135 (MQ=255) ccGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCGAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCg < 7:292312/148‑1 (MQ=255) tCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTg > 2:17913/1‑149 (MQ=255) tCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTg > 1:238591/1‑149 (MQ=255) aCTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCGAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGa > 8:53754/1‑148 (MQ=255) | TTGCTTCGTTCTAGAGCCTCTGCGGCCAGATTGCAGTACCTGATGCACATGCAACTCTAACTATAGTCAGGATGACACCTGATGACCTATTCGAGGATCCATGCCATGAGCGATC‑ATGACGAACGCCGCCGCTTTCAACGCATCGACTTCGAT‑GCCCCCACCGAGCTGCGCCAGGGGCTGAAACGCTGGCCTGCAAAACTGCTGGACCTTTCACTCAAAGGCTTGCTGATCGAACGCCCGGCAGACTGGGATGCCGACCTGACTCAGGATTTCGACGCCATCATCGACCTGGA > NC_002947/5266684‑5266976 |
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: ‑ |