Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,955,845 | T→G | 7.2% | intergenic (+71/+133) | elbB → / ← PP_5221 | putative hydrolase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,955,845 | 0 | T | G | 7.2% | 155.6 / 2.2 | 52 | intergenic (+71/+133) | elbB/PP_5221 | putative hydrolase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (26/22); new base G (2/2); total (28/24) | |||||||||||
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.87e-01 |
GAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGTCAGCCTTTC > NC_002947/5955697‑5955987 | gAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGt < 4:161055/149‑1 (MQ=255) cAAGTCGATCAGCGAGGCGGCGGGCGGGATCTAGAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtg > 2:57792/1‑149 (MQ=255) gTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTgg > 6:374597/1‑149 (MQ=255) aTCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGc > 5:68797/1‑148 (MQ=255) tCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACCAGGAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCAACAGGTCCAGCGCAGGCTTCAAGGCCAAGCGGTACCTGTGGGGCa > 7:7168/1‑148 (MQ=255) cGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGATTCTTCGCAGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTgc > 4:333458/1‑149 (MQ=255) gAGGCGGCGGGCGCGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCCCAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTCCCTTAGCGGCAACTgcg > 8:165038/1‑149 (MQ=38) cgggATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTgg < 6:309210/128‑1 (MQ=25) cgggATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTgg > 5:309210/1‑128 (MQ=25) gggATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCtgtg < 6:330071/149‑1 (MQ=255) ggATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGCCGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCCGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTgg > 8:323845/1‑149 (MQ=255) tACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtgtg > 8:204495/1‑120 (MQ=255) tACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtgtg < 7:204495/120‑1 (MQ=255) aGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTgc < 5:374597/149‑1 (MQ=255) cTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTgg‑‑‑g > 8:47322/1‑139 (MQ=255) cTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTgg‑‑‑g < 7:47322/139‑1 (MQ=255) gtggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGAt > 7:149075/1‑149 (MQ=255) ccGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCtt > 5:350974/1‑149 (MQ=255) aGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTg > 1:225616/1‑147 (MQ=255) gTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTa < 2:225616/148‑1 (MQ=255) tGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAgg < 7:168678/149‑1 (MQ=255) tGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGc > 5:216261/1‑149 (MQ=255) aCGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGAGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTAc < 6:133693/149‑1 (MQ=255) aGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGC‑‑‑GGCAACTCCGCGATCGCTTGTAGGAGCGGGATTACCCGCGaa < 1:57792/149‑1 (MQ=255) ttACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGAt < 1:168461/92‑1 (MQ=255) ttACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGAt > 2:168461/1‑92 (MQ=255) aCTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCgg < 7:349456/147‑1 (MQ=255) cTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCa < 5:296512/148‑1 (MQ=255) gCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTg > 5:252965/1‑149 (MQ=255) gCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCAGAGAACACCGGCAGAGCCGGTg > 2:64286/1‑149 (MQ=255) gCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGAGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTg > 6:335455/1‑149 (MQ=255) cttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGcc < 6:350974/147‑1 (MQ=255) tcGCGGGCATGCCCGCTGCCACAGGTCCAGCGCTGGCTTCTAGGCCACGCGGTACCTGTGTGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTTGGAGCGGGTTTACCCGCGAACACCGGCTCTGCCGGTgttcgc < 8:7168/148‑6 (MQ=255) gcgGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCAc > 6:329968/1‑149 (MQ=255) gggCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGAGGGGAGGGATGTGGACGAAAGATTTTATGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCa > 1:122493/1‑149 (MQ=25) aTGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCgcgc > 1:166351/1‑148 (MQ=255) tGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCa < 7:323845/149‑1 (MQ=255) cGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCt < 8:149075/149‑1 (MQ=255) tCCCTCTGGTCCAGGGCAGGCTTCAAGGCCACGCGGGACCTGTGGGGCAACTGCGCTATATCCTGTGG‑‑‑GGCAACTGCGCGATCGCGTGTAGGAGGGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTggg < 7:165038/149‑1 (MQ=37) cccACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGc > 8:283539/1‑149 (MQ=25) acaGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTg < 2:26266/123‑1 (MQ=255) acaGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTg > 1:26266/1‑123 (MQ=255) acaGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGt < 6:216261/149‑1 (MQ=255) cAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTa > 6:137788/1‑148 (MQ=255) aGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGTGCTGTCGCCTACa < 5:329968/149‑1 (MQ=255) ttCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGt > 8:82725/1‑149 (MQ=255) tCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCAATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGTAATCTGt > 8:329853/1‑148 (MQ=255) ggCCACGCGGTACCTGTGGGGAAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGTCAGcc > 2:89433/1‑149 (MQ=255) gCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGTCAGCCt < 2:122493/149‑1 (MQ=255) cACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGt > 2:25898/1‑70 (MQ=255) cACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGt < 1:25898/70‑1 (MQ=255) cACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGTCAGCCttt < 6:143593/149‑1 (MQ=255) aCGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGTCAGCCTTTc < 6:68797/149‑1 (MQ=255) | GAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGG‑‑‑GGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTGTCAGCCTTTC > NC_002947/5955697‑5955987 |
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: ‑ |