Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,049,910 | C→G | 23.4% | intergenic (+43/‑58) | 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,910 | 0 | C | G | 23.4% | 60.9 / 15.7 | 34 | intergenic (+43/‑58) | quiA/oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base C (10/16); new base G (5/3); total (15/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAG > NC_002947/4049762‑4050056 | cTACGTATCGCCGCAGACCGGCAAGCAGTACGTGGTGATCACAGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGGCGCCTGCTTAGCGGGCCACCCCGCGCCTACAACGc > 6:362627/1‑149 (MQ=255) cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAGAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGgcc < 5:172085/148‑1 (MQ=255) tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTg > 5:266696/1‑149 (MQ=255) gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa < 1:132163/147‑1 (MQ=255) gAAGACCGTCAAGCAGTACGTGGTGATCACTGCTGGCGGCTCCCGCAAGTCGACTGACCGTGGCGACTATGTGATTTGTTACGCTCTGCCGTATAACCTCGTCGCCTGATTCTCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc < 6:355149/149‑1 (MQ=37) gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc < 6:266696/149‑1 (MQ=255) gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg < 2:316134/149‑1 (MQ=255) cGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGGCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc < 1:342008/149‑1 (MQ=255) gtggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGCCGACAAGGGCCCAGTAGGAGCGGGCGTGGCCGCGaaaggc > 2:119318/1‑145 (MQ=255) gggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGGGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGccc < 5:291346/148‑1 (MQ=255) tCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGATTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa < 1:119318/149‑1 (MQ=255) tCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCAGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa < 4:221869/149‑1 (MQ=255) ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑c < 5:362627/149‑1 (MQ=255) aCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc < 6:146950/97‑1 (MQ=255) aCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc > 5:146950/1‑97 (MQ=255) aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg < 6:199090/95‑1 (MQ=255) aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg > 5:199090/1‑95 (MQ=255) gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTg < 5:179934/128‑1 (MQ=255) gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTg > 6:179934/1‑128 (MQ=255) gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGTGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGGAATGCCATCCGCAATTCGCATCAccccc > 6:184474/1‑149 (MQ=255) gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAccc < 1:59759/148‑1 (MQ=255) ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggcc < 5:165836/148‑1 (MQ=255) cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGTGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGccc < 5:184474/147‑1 (MQ=255) cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCCGCTGGCCAGGGCCAgcagcagggtgg < 3:71324/147‑7 (MQ=38) acgcgGGCAAGCCCGCTCCTACAAGGCCCCTTGAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc < 7:280605/148‑1 (MQ=255) aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac > 1:149698/1‑148 (MQ=255) ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct > 7:76395/1‑149 (MQ=255) ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCCCCACTGCCCTTGCCAGCGACGAtct > 2:251748/1‑149 (MQ=255) ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct > 1:284867/1‑149 (MQ=255) cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct > 8:251769/1‑148 (MQ=255) cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc < 8:76395/149‑1 (MQ=255) ccTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa > 4:296783/1‑149 (MQ=255) aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc > 5:87852/1‑127 (MQ=255) aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc < 6:87852/127‑1 (MQ=255) aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCAGCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCCCTGCCCTTGCCAGCGCCGATCTCAGGACCCgc > 4:247806/1‑149 (MQ=255) gCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcag > 5:9219/1‑148 (MQ=255) | CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAG > NC_002947/4049762‑4050056 |
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: ‑ |