Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,936,634 | A→T | 7.9% | K119N (AAA→AAT) | rbsB → | D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,936,634 | 0 | A | T | 7.9% | 186.1 / 4.8 | 62 | K119N (AAA→AAT) | rbsB | D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor |
Reads supporting (aligned to +/- strand): ref base A (27/30); new base T (3/2); total (30/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
GCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGA > NC_000913/3936494‑3936780 | gCGACCGTGCAGGACTTACCCGTTCGCGGCACAAAAATTCTGCTGATTTACCCGCCCGCCTCCGACGCAGTGGGTAATGGTGTGAAGATGGCTAACCAGGCTATCATACCGGGTATCTCTCGTGGGCTCCAGGCAAGGAATGGGGAAg > 6:118513/1‑148 (MQ=255) gCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGa < 3:56646/146‑1 (MQ=255) cGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTTGGTAATGCTGTGAAGATTGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAgt < 5:118513/148‑1 (MQ=255) aaCGTGCAGGACTTAACCGATCGCGGAACAAAAATACTGCTGATTAACCCGACCGACTCCGACGCAGAGGGAAAAGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCGAGGCAACGAAAGGTGAAgtgg > 5:335375/1‑148 (MQ=255) ggACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCcac < 2:110810/148‑1 (MQ=255) gTTCGCGGCACAAAAATTCTGCTGATTAACCCGCCCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCt > 4:274493/1‑147 (MQ=255) tCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGa > 4:350892/1‑147 (MQ=255) gcgGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATaa < 2:302766/148‑1 (MQ=255) gcgGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATa < 3:33036/147‑1 (MQ=255) aaaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGa < 1:142117/114‑1 (MQ=255) aaaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGa > 2:142117/1‑114 (MQ=255) aaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTgg > 1:168888/1‑148 (MQ=255) aaaaTTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTgg > 6:61800/1‑148 (MQ=255) tAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGTACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt > 3:290326/1‑148 (MQ=255) aCCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg < 5:61800/147‑1 (MQ=255) aCCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGGGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg < 3:350892/147‑1 (MQ=255) cccGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGt < 3:274493/148‑1 (MQ=255) ccgaccgaCTCCTACGCAGTGGGTAATGCTGTTAAGATGGCTAACCAGGGGAACATCCCGGTATTCTCTCTTGACCGCCAGGCTACGAAAGGTTAAGTGGTTAGCCACATTGCTTCTGATAACGTACTGGGCGGCATAATCGCTGGTg < 6:335375/148‑1 (MQ=255) ccgaccgaCTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTg > 6:65898/1‑118 (MQ=255) ccgaccgaCTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTg < 5:65898/118‑1 (MQ=255) gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAg < 2:112099/84‑1 (MQ=255) gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAg > 1:112099/1‑84 (MQ=255) gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg > 7:21067/1‑148 (MQ=255) gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg > 5:250624/1‑148 (MQ=255) gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg > 4:194690/1‑148 (MQ=255) tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGt < 8:21067/148‑1 (MQ=255) tGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa < 6:250624/148‑1 (MQ=255) cTGTGAAGATGGCTAACCAGGCGTACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg < 4:290326/147‑1 (MQ=255) tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt > 4:337388/1‑148 (MQ=255) tGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt > 1:309979/1‑148 (MQ=255) gCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa > 3:69819/1‑148 (MQ=255) cTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa < 5:62875/109‑1 (MQ=255) cTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa > 6:62875/1‑109 (MQ=255) aCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAATGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCg < 6:385063/147‑1 (MQ=255) aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCt > 1:155550/1‑148 (MQ=255) cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCAtt > 6:205746/1‑148 (MQ=255) tCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc > 8:233937/1‑137 (MQ=255) tCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc < 7:233937/137‑1 (MQ=255) cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTAc > 2:262720/1‑147 (MQ=255) ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa > 6:179756/1‑104 (MQ=255) ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa < 5:179756/104‑1 (MQ=255) ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa > 1:346532/1‑122 (MQ=255) ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa < 2:346532/122‑1 (MQ=255) ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGt > 8:109498/1‑141 (MQ=255) ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGt < 7:109498/141‑1 (MQ=255) ttGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATcc < 3:337388/146‑1 (MQ=255) aCCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc < 8:248659/128‑1 (MQ=255) aCCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGc > 7:248659/1‑128 (MQ=255) gCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCg < 5:205746/148‑1 (MQ=255) aGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGa < 5:206106/148‑1 (MQ=255) aGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGa < 4:69819/148‑1 (MQ=255) ggCAACGAATGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg > 8:279740/1‑95 (MQ=255) ggCAACGAATGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg < 7:279740/95‑1 (MQ=255) ggCAACGAATGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGAGAAGAAAGCGGGTGAAg > 8:279739/1‑95 (MQ=255) ggCAACGAATGGTGAAGTGGTGAGACACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg < 7:279739/95‑1 (MQ=255) ggCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCAATGGTCGTACATCCGCAGCCCGTGaa > 6:206106/1‑148 (MQ=255) gCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAAc < 2:168888/148‑1 (MQ=255) gCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAAc < 3:194690/148‑1 (MQ=255) cGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTgg > 4:150578/1‑148 (MQ=255) gAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGaa < 7:102635/142‑1 (MQ=255) gAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGaa > 8:102635/1‑142 (MQ=255) gAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGc < 2:309979/148‑1 (MQ=255) aaGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGa < 3:150578/148‑1 (MQ=255) aaGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGa < 2:155550/148‑1 (MQ=255) | GCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTGCAAGGCATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGA > NC_000913/3936494‑3936780 |
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: ‑ |