Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,049,913 | C→T | 15.0% | intergenic (+46/‑55) | quiA → / → oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,049,913 | 0 | C | T | 15.0% | 118.1 / 12.8 | 54 | intergenic (+46/‑55) | quiA/oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base C (24/21); new base T (3/5); total (27/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01 |
CTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACC > NC_002947/4049768‑4050060 | ctcGCCGAAGACCGGCAAGTAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGt > 7:265815/1‑149 (MQ=255) tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTa > 5:180883/1‑149 (MQ=255) cGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAg > 6:80665/1‑149 (MQ=255) gCCGAAGACCTTCAAGCAGTAGGTGGTGATCTCTGGTGGCGGAGCCCGCCAGTCGACTGACCGTGGCTACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGA‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGACCCTTGttg < 4:88802/148‑3 (MQ=255) gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAGGGGCCTTGTAGGAGc < 3:78323/149‑1 (MQ=255) gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAGGGGCCTTGTAGGAGc < 5:58167/149‑1 (MQ=255) gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGc > 7:290521/1‑149 (MQ=255) gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTGCGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCgg < 1:75878/149‑1 (MQ=255) cGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCtt > 3:285059/1‑149 (MQ=255) gCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAGGGGCCTTGTAGGAGCGGGCTTGc < 5:164249/149‑1 (MQ=255) gATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAGGGGCCTTGTAGGAGCgg < 4:273827/127‑1 (MQ=255) gATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAGGGGCCTTGTAGGAGCgg > 3:273827/1‑127 (MQ=255) tCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaa < 8:265815/148‑1 (MQ=255) cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCCAGCCCGCACC‑GACAAGGCCCCGGGAGGCGCGGGGGGGGCCGCGACTGGCCCGCAAAg > 3:189669/1‑149 (MQ=255) cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCa > 6:71935/1‑146 (MQ=255) aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGa < 5:12809/149‑1 (MQ=255) ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAAc > 2:238905/1‑149 (MQ=255) cggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCga < 5:71935/149‑1 (MQ=255) ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCtt < 8:71263/113‑1 (MQ=255) ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCtt > 7:71263/1‑113 (MQ=255) ccAGTCGACTGACCGTGGCTACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGTAGCGGGCTTGccc < 2:5423/114‑1 (MQ=255) ccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGccc > 1:5423/1‑114 (MQ=255) aGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAACCCGCACCGTAAACGGCCCC‑GTAGGCGGCGGCGTGGCCCGGACTGCCCCCCAAAGACCCACCCAGATTGTGGAAt > 6:195358/1‑149 (MQ=255) gTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATg > 7:279865/1‑149 (MQ=255) aCCGTGGCTACTATGTTATTCCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGGTCC‑AACAAGCCCCCTGTCGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCg < 3:311185/148‑1 (MQ=255) aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCg < 1:238905/148‑1 (MQ=255) tgATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCCAGCCCGCGCC‑GACCAGGGCCCGGGAGGGGGGGGGGTGGCCGGGGAGGGCCCGGAAAGACCACCCGGGGGTCTGGAATGCCAAGCGCCATTGGCCGCAcccc > 4:295019/1‑149 (MQ=255) tgATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCgg < 3:306622/80‑1 (MQ=255) tgATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCgg > 4:306622/1‑80 (MQ=255) gATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAcccc < 6:180883/148‑1 (MQ=255) cTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCct > 7:75452/1‑149 (MQ=255) tctGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTAAAAGACAACCCGAGATTCTGTAATGCCCTCCGCAATTCGGATCACCCCCACCctgctgctg > 7:146474/1‑149 (MQ=255) ccGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGAc < 7:52944/90‑1 (MQ=255) ccGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggcc < 4:285059/148‑1 (MQ=255) ccGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCCCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggcc > 8:97547/1‑148 (MQ=255) ccGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCCCTCC‑AACAAGGCCCCAGTAGGAGCGGGCTTGCCCGCGAATGGCCCGCAAAGAc > 8:52944/1‑90 (MQ=25) aaaCCTCGT‑CGCGTGCTTCGGGGGCAAGCCGGCTCC‑TACAAGGCGGGTGTAGGATCGGGCTTGCGCGCTAATGCCCCTGAAAGACAACCCGAGATTCTGTAATGCGATCCGCAATTCGCATCACCCCCACGCTGCTGctggccctggc < 4:118075/148‑1 (MQ=255) aaCCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa < 5:132419/149‑1 (MQ=255) cgcgT‑CGCCTGCTTCGCGGGCAACCCCGCTCC‑AACAAGGCCCCTGTCGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGca < 5:195358/149‑1 (MQ=255) gcgT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac < 3:295019/149‑1 (MQ=255) gT‑CGCCGGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagca < 4:189669/149‑1 (MQ=38) cTGCTTCGCGGGCAACCCCCCACC‑AACAAAGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCACAAAAACAAACCCAGAATATGTAAATCCAAACGCAATTATCATAACCCCCCCCCAGCAGCTGGCCCAGGCcaccacaccagcc > 8:151295/1‑142 (MQ=255) cTGCTACGCGGCCACCCCCGCTCA‑AACACGCCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcc < 5:276875/149‑1 (MQ=255) gCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGCCCCTCAAAGCCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCGGGCCCTGGCCAGCACCCCTGCCCt > 5:224505/1‑149 (MQ=255) ttCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagcagg < 2:97554/142‑3 (MQ=255) cgcgGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc < 8:279865/148‑1 (MQ=255) aaGCCCGCTCC‑TACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac > 2:148047/1‑148 (MQ=255) ccGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct > 1:290964/1‑149 (MQ=255) gCTCC‑TACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa > 3:307629/1‑149 (MQ=255) gCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa > 1:221646/1‑149 (MQ=255) gCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa < 8:290521/149‑1 (MQ=255) c‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGACCTGGCCAGCACCACTGCCCTTGGCAGCGACGATCTCATGAc > 2:182283/1‑149 (MQ=255) cAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCAATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg > 8:187349/1‑149 (MQ=255) aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgc < 1:166950/149‑1 (MQ=255) ccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc < 2:290964/149‑1 (MQ=255) ccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATACGCAATTCGCATCACCCCCACCCTCCTGCTGGCCCTGGCCAGCACCACTGCCCTTGACAGCCACGAACTCATGACCCgcaccccc > 2:318953/1‑144 (MQ=255) | CTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGT‑CGCCTGCTTCGCGGGCAAGCCCGCTCC‑TACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACC > NC_002947/4049768‑4050060 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |