Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,061,993 | T→C | 5.4% | Y277H (TAT→CAT) | yihM → | putative sugar phosphate isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,061,993 | 0 | T | C | 5.4% | 261.5 / 3.9 | 74 | Y277H (TAT→CAT) | yihM | putative sugar phosphate isomerase |
Reads supporting (aligned to +/- strand): ref base T (33/37); new base C (2/2); total (35/39) | |||||||||||
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.01e-01 |
GCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCAAACCATCTTC > NC_000913/4061851‑4062134 | gCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc < 2:364871/148‑1 (MQ=255) gCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc > 1:402086/1‑148 (MQ=255) ccGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc > 1:318122/1‑148 (MQ=255) gTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAg < 3:939286/148‑1 (MQ=255) ttCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGa < 1:916491/148‑1 (MQ=255) aaGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAg < 3:988272/148‑1 (MQ=255) aGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGt < 2:318122/148‑1 (MQ=255) gTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGaa < 4:346967/148‑1 (MQ=255) aCGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc > 4:539802/1‑130 (MQ=255) aCGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc < 3:539802/130‑1 (MQ=255) ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg > 3:598976/1‑122 (MQ=255) ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg > 1:601355/1‑122 (MQ=255) ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg < 2:601355/127‑6 (MQ=255) ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg < 4:598976/127‑6 (MQ=255) aTCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaaa < 2:402086/148‑1 (MQ=255) gaGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa < 2:579516/113‑1 (MQ=255) gaGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa > 1:579516/1‑113 (MQ=255) ggTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCa > 3:954827/1‑112 (MQ=255) ggTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCa < 4:954827/112‑1 (MQ=255) atatGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGa < 3:632979/92‑1 (MQ=255) atatGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGa > 4:632979/1‑92 (MQ=255) ggCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGt < 1:1058529/130‑1 (MQ=255) ggCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGt > 2:1058529/1‑130 (MQ=255) cTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaa > 2:381437/1‑105 (MQ=255) cTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaa < 1:381437/105‑1 (MQ=255) ggAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc > 1:476166/1‑76 (MQ=255) ggAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc < 2:476166/76‑1 (MQ=255) aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACAcc > 1:283520/1‑91 (MQ=255) aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACAcc < 2:283520/91‑1 (MQ=255) aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaaga > 3:20862/1‑133 (MQ=255) aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaaga < 4:20862/133‑1 (MQ=255) aagaGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAaca < 2:1017759/86‑1 (MQ=255) aagaGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAaca > 1:1017759/1‑86 (MQ=255) aagaGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTaa > 2:1021577/1‑148 (MQ=255) agagGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAAt < 1:1021577/148‑1 (MQ=255) tGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTa > 4:56459/1‑100 (MQ=255) tGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTa < 3:56459/100‑1 (MQ=255) tACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTg > 4:1076863/1‑81 (MQ=255) tACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTg < 3:1076863/81‑1 (MQ=255) aTGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc < 4:386674/93‑1 (MQ=255) aTGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc > 3:386674/1‑93 (MQ=255) ttCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTaa < 1:744303/122‑1 (MQ=255) ttCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTaa > 2:744303/1‑122 (MQ=255) ccGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGtacta > 1:204211/1‑148 (MQ=255) cGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGtactac > 4:318926/1‑148 (MQ=255) tttGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGAt < 1:1020011/108‑1 (MQ=255) tttGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGAt > 2:1020011/1‑108 (MQ=255) tGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTAcaacaa < 2:204211/148‑1 (MQ=255) gACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAAt < 3:882464/111‑1 (MQ=255) gACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAAt > 4:882464/1‑111 (MQ=255) gACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCa > 1:1082890/1‑148 (MQ=255) aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCAtg < 1:1055916/123‑1 (MQ=255) aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCAtg > 2:1055916/1‑123 (MQ=255) aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTAcaacaa < 1:501793/143‑1 (MQ=255) aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTAcaacaa > 2:501793/1‑143 (MQ=255) aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCaa < 3:318926/148‑1 (MQ=255) gatgatAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAag < 2:1082890/148‑1 (MQ=255) tgatAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAagag > 2:338845/1‑148 (MQ=255) atAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATc < 3:1140631/102‑1 (MQ=255) atAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATc > 4:1140631/1‑102 (MQ=255) aaTCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCa > 1:892522/1‑110 (MQ=255) aaTCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCa < 2:892522/110‑1 (MQ=255) aTCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCaaa > 1:767449/1‑148 (MQ=255) tCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAagag < 2:322839/142‑1 (MQ=255) tCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAAGTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAagag > 1:322839/1‑142 (MQ=255) ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaaagaa > 4:562758/1‑84 (MQ=255) ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaaagaa < 3:562758/84‑1 (MQ=255) ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCACTTAATCAGATAAATCAtg < 2:863786/109‑1 (MQ=255) ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCAtg > 1:863786/1‑109 (MQ=255) tGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaa > 2:569300/1‑77 (MQ=255) tGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaa < 1:569300/77‑1 (MQ=255) aTCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCAAACCAtcttc < 1:338845/148‑1 (MQ=255) tCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGta > 2:277984/1‑113 (MQ=255) tCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGta < 1:277984/113‑1 (MQ=255) | GCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCAAACCATCTTC > NC_000913/4061851‑4062134 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |