| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,250,629 | (C)5→6 | coding (210/1071 nt) | dhaK ← | dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,250,624 | 1 | . | C | 97.9% | 164.1 / ‑3.4 | 47 | coding (215/1071 nt) | dhaK | dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit |
| Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (22/24); total (23/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01 | |||||||||||
CGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCA > NC_000913/1250487‑1250756 | cGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑c < 1:277708/139‑1 (MQ=255) cTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑cc‑ccccg < 2:1362836/139‑3 (MQ=255) tGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑cc‑ccccga > 2:115970/1‑136 (MQ=255) tGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGACCGGACAGG‑CCACCCGa > 2:1313843/1‑139 (MQ=255) gTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATccc < 1:1382894/139‑1 (MQ=255) aaGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTg < 2:837032/139‑1 (MQ=255) aaGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTg > 1:837032/1‑139 (MQ=255) gAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAc < 1:292735/139‑1 (MQ=255) cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGa < 1:592158/132‑1 (MQ=255) cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGa > 2:592158/1‑132 (MQ=255) ggTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAAcc < 2:513696/139‑1 (MQ=255) gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACca > 1:399489/1‑139 (MQ=255) gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑cc‑ccc < 2:95211/89‑1 (MQ=255) gTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑cc‑ccc > 1:95211/1‑89 (MQ=255) ccTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc > 2:782471/1‑102 (MQ=255) ccTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc < 1:782471/102‑1 (MQ=255) cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAgg > 2:512064/1‑78 (MQ=255) cTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAgg < 1:512064/78‑1 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc < 2:1141761/99‑1 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc > 1:1141761/1‑99 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa > 2:226618/1‑139 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa > 1:296285/1‑139 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc < 1:1174365/96‑1 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc < 2:573261/96‑1 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc > 2:1174365/1‑96 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc > 1:573261/1‑96 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc < 1:209553/134‑1 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc > 2:209553/1‑134 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc > 2:628103/1‑139 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc > 2:357527/1‑139 (MQ=255) gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGa < 1:190515/81‑1 (MQ=255) gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGa > 2:190515/1‑81 (MQ=255) gcgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAcc < 2:1063192/120‑1 (MQ=255) gcgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAcc > 1:1063192/1‑120 (MQ=255) cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCg < 1:478514/108‑1 (MQ=255) cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCg > 2:478514/1‑108 (MQ=255) gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc < 2:399489/139‑1 (MQ=255) ttttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAgggg > 2:384317/1‑139 (MQ=255) ttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCa > 2:1431242/1‑139 (MQ=255) tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCAt > 1:224797/1‑139 (MQ=255) gggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGc < 2:461739/139‑1 (MQ=255) ggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc > 2:894575/1‑76 (MQ=255) ggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc < 1:894575/76‑1 (MQ=255) gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc > 2:1312431/1‑56 (MQ=255) gCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGAcc < 1:1312431/56‑1 (MQ=255) cGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCg > 2:445494/1‑139 (MQ=255) gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACATTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg < 2:296285/139‑1 (MQ=255) gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg < 1:628103/139‑1 (MQ=255) tgaggtgaAAATTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGAc < 1:445494/139‑1 (MQ=255) aaaTTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg < 1:1431242/139‑1 (MQ=255) aTTTCGCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGa < 1:1313843/139‑1 (MQ=255) gCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACttt < 1:436320/94‑1 (MQ=255) gCCCGGACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACttt > 2:436320/1‑94 (MQ=255) gACAGGCCC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCa < 1:357527/139‑1 (MQ=255) | CGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGG‑CC‑CCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCA > NC_000913/1250487‑1250756 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |