Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,049,957 | C→T | 25.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,957 | 0 | C | T | 25.0% | 45.9 / 15.5 | 32 | intergenic (+46/‑55) | quiA/oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base C (11/13); new base T (3/5); total (14/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01 |
ACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCT > NC_002947/4049808‑4050106 | aCGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGgccc < 2:222090/150‑1 (MQ=255) aCGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGACAGCCCGCGCCTACCCGGcccc > 2:1544098/1‑150 (MQ=255) gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa < 2:1549927/150‑1 (MQ=255) gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa < 2:1468228/150‑1 (MQ=255) ccGAAGACCGCCGAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGGCGGTCGACTGACCGGGGCGGCTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAg < 2:534783/150‑1 (MQ=38) cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc < 2:1749990/150‑1 (MQ=255) aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg < 2:1255889/150‑1 (MQ=255) aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg < 2:1227034/150‑1 (MQ=255) aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGAATATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg < 2:1591194/150‑1 (MQ=255) aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCCGCAGTCTACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg < 1:409109/150‑1 (MQ=255) aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg < 2:581610/150‑1 (MQ=255) aGACCGGCAAGCAGTACGTGGTGATCAATGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGCCCCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg < 2:1011196/150‑1 (MQ=255) gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc > 1:1202510/1‑150 (MQ=255) ggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAg < 2:248760/111‑1 (MQ=38) ggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAg > 1:248760/1‑111 (MQ=38) tACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgct < 2:624594/150‑1 (MQ=255) cGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgc < 1:1544098/150‑1 (MQ=255) gcgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcacc > 1:1344283/1‑150 (MQ=255) gggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcg > 1:1796342/1‑150 (MQ=255) aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga < 1:1098570/150‑1 (MQ=255) cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCCTCAAAGCCAACCCGAGATTCTGTAATGCCATCCGCAATTCCCATCACCCCCACCCTGCTGCTGGCCCTGGCAACCACCACTCCCCTTCCCAGcgacgacc > 2:1746384/1‑148 (MQ=255) ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc > 2:1628198/1‑150 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAAcc < 2:1836899/58‑1 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAAcc > 1:1836899/1‑58 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGCTCTCAt > 2:94539/1‑150 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt > 1:522202/1‑150 (MQ=255) gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCAAGCGACGATCTCag > 2:110508/1‑149 (MQ=255) cTCCTACAGGGGCCTTGTAGGAGCGGGCTTG‑CCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa > 1:483871/1‑150 (MQ=255) cTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc > 2:372512/1‑150 (MQ=255) cTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccag < 2:287184/110‑3 (MQ=255) cTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccag > 1:287184/1‑108 (MQ=255) aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg < 2:1344283/150‑1 (MQ=255) aGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcag < 1:94539/150‑1 (MQ=255) gCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCCTGACCCgcaaca > 2:145911/1‑150 (MQ=255) cccTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcagca > 2:1498024/1‑148 (MQ=255) cccTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcagca < 1:1498024/148‑1 (MQ=255) cTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCt > 2:839210/1‑150 (MQ=255) | ACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATG‑CCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCT > NC_002947/4049808‑4050106 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |