| 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 | 98.2% | 192.7 / ‑3.4 | 57 | coding (215/1071 nt) | dhaK | dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (23/33); total (23/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01 | |||||||||||
GTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑G‑CCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGC > NC_000913/1250489‑1250755 | gTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑G‑ccc > 2:1474910/1‑139 (MQ=255) gCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑G‑cccccc > 1:992148/1‑138 (MQ=255) gTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATccc < 1:1474910/139‑1 (MQ=255) aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa < 2:646471/139‑1 (MQ=255) aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa < 2:1223675/139‑1 (MQ=255) aGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGAGTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa < 2:1308531/139‑1 (MQ=255) aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc < 1:60409/139‑1 (MQ=255) aaTATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc < 1:1234792/139‑1 (MQ=255) tatCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGa < 1:906089/139‑1 (MQ=255) tCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc < 2:451167/135‑1 (MQ=255) tCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc > 1:451167/1‑135 (MQ=255) tCGCCGGTGTAATTTTTGATAAGCAACAGTACACCTTCGCCGCCATCAACTTGAATGGCGCATTCAAAGATTTTATCCGGCGTCGGAGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGata > 2:753767/1‑139 (MQ=255) cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc > 2:113349/1‑134 (MQ=255) cGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc < 1:113349/134‑1 (MQ=255) cGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAAc < 2:697069/139‑1 (MQ=255) gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑G‑ccccc < 1:1373070/111‑1 (MQ=255) gtgtAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑G‑ccccc > 2:1373070/1‑111 (MQ=255) aaTTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGt < 2:992148/139‑1 (MQ=255) ttGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc < 1:1399128/139‑1 (MQ=255) gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc > 2:853931/1‑137 (MQ=255) gATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATc < 1:853931/137‑1 (MQ=255) tCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCgt > 2:906439/1‑139 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCgtgt < 2:1424830/127‑1 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCgtgt > 1:1424830/1‑127 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc > 1:45635/1‑134 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc < 2:45635/134‑1 (MQ=255) tCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCa > 1:59872/1‑139 (MQ=255) cgccgcCATCAACTTGCGTGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGACAATTTCGCCCGGAGAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAc > 2:879196/1‑139 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa > 2:1290949/1‑94 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa < 1:1290949/94‑1 (MQ=255) ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCAccc < 1:906439/139‑1 (MQ=255) aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGc > 2:290005/1‑80 (MQ=255) aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGc < 1:290005/80‑1 (MQ=255) aCTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAgg < 2:59872/139‑1 (MQ=255) gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa > 1:966561/1‑81 (MQ=255) gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGa < 2:966561/81‑1 (MQ=255) gCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGa > 1:584370/1‑139 (MQ=255) cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagc < 2:786123/126‑1 (MQ=255) cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagc > 1:786123/1‑126 (MQ=255) gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGAC‑CCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc < 2:646700/139‑1 (MQ=255) gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc < 2:122933/76‑1 (MQ=255) gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc > 1:122933/1‑76 (MQ=255) aTTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTg < 1:879196/139‑1 (MQ=255) gATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAgg > 1:1064459/1‑139 (MQ=255) ttttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc > 2:1528396/1‑103 (MQ=255) ttttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc < 1:1528396/103‑1 (MQ=255) tATCGGGCGTCGGTGA‑GTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACca > 1:948665/1‑71 (MQ=255) tATCGGGCGTCGGTGA‑GTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACca < 2:948665/71‑1 (MQ=255) gCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTc > 1:918658/1‑139 (MQ=255) gTCGGTGAGGTGAAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg < 2:1484652/139‑1 (MQ=255) ggtgaAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc < 2:1321323/46‑1 (MQ=255) ggtgaAAATTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGAcc > 1:1321323/1‑46 (MQ=255) aaaaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAAc < 2:597429/50‑1 (MQ=255) aaaaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAAc > 1:597429/1‑50 (MQ=255) aaaaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc < 2:1530839/80‑1 (MQ=255) aaaaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc > 1:1530839/1‑80 (MQ=255) aaaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCg < 2:918658/139‑1 (MQ=255) aaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGa > 1:1020966/1‑42 (MQ=255) aaTTTCGCCCGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGa < 2:1020966/42‑1 (MQ=255) cGGACAG‑GCCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATg < 2:1064459/139‑1 (MQ=255) ggACAG‑GCCCCCCGAAACCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCTGCCACCCGACAGCAGGGCGACTTTTCCTGCATCAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATgg < 1:942500/139‑2 (MQ=255) | GTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATCAACAGTACACCTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAG‑G‑CCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGC > NC_000913/1250489‑1250755 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |