Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,836,062 | T→A | 11.0% | Q114L (CAG→CTG) | hycI ← | protease involved in processing C‑terminal end of HycE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,836,062 | 0 | T | A | 11.0% | 86.6 / 4.1 | 35 | Q114L (CAG→CTG) | hycI | protease involved in processing C‑terminal end of HycE |
Reads supporting (aligned to +/- strand): ref base T (15/16); new base A (2/2); total (17/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 = 4.61e-01 |
AATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCC > CP009273/2835923‑2836205 | aaTGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCACCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGa < 1:308002/149‑1 (MQ=255) gAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGcc > 2:332295/1‑144 (MQ=255) gAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGcc > 2:332199/1‑144 (MQ=255) ccACCGCTAACGGCGCGAAGCCGCCATTTCCTGCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGc < 2:235885/144‑1 (MQ=255) tAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGGGGGAAATACCTTCGCCAATCTcttctt > 1:49664/1‑149 (MQ=255) caccTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATcttctt < 1:35614/146‑1 (MQ=255) gAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAGTCACTTCGCCAATATCTTCTTTCaa < 2:235830/144‑1 (MQ=255) caGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTc < 1:164968/148‑1 (MQ=255) tccCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTaa < 2:321587/144‑1 (MQ=255) ccTTCCAGTCGTTGATAAACGGGTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGGAAATCACTTCGCCAATATGTTATTTCAACTGGTCGATAGGATACATCAACGACAt > 1:110722/1‑149 (MQ=255) ttCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACgg > 2:189736/1‑144 (MQ=255) ccAGTCGTTGATAAACGGTTTCGACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCa < 2:319147/144‑1 (MQ=255) cAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATAACTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTAt > 1:28126/1‑149 (MQ=255) gTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCCACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAAGATCTTCTTTCAACTGGTCGATCAGGTAATTCAACGGCagtatag > 2:136586/1‑138 (MQ=255) aTAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCATCTGGTCGATAAGGTAATTCAACGGCATGTTATggga > 2:207161/1‑143 (MQ=255) tAAACGTTTTCTACCTCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTa < 2:156135/144‑1 (MQ=255) aaCGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGtcatca < 1:237333/149‑1 (MQ=255) aCGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGGTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCCACTGGGCGATAAGGTACTTAAACGGCATGTTATGGTTGGTcatcat > 1:66006/1‑149 (MQ=255) aCGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACGGGTCGATAAGGTCATTTAACGGCATGTTATGGGTAGtc > 2:324424/1‑144 (MQ=255) gTTTGTACCGCATCTTTAATCGGCTGGGTCGTCGGGTAGTAAAGGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGCTCGATAAGGTAATTCATCGGCATGTTATGGGTAGtcatc < 2:32216/144‑1 (MQ=255) gTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGtcatc < 2:317395/144‑1 (MQ=255) ccGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACa > 2:245445/1‑144 (MQ=255) ctCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTTTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGa < 1:19547/148‑1 (MQ=255) tCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGat < 1:368659/149‑1 (MQ=255) tCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGat < 1:368561/149‑1 (MQ=255) tAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGACTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTAGATAAGGTAATTCAACGGCCTGTTATGGGTAGTCATCATACACATCTGGGCgg > 2:154817/1‑143 (MQ=255) gCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATatcat > 2:35649/1‑144 (MQ=255) tCATCGGGTAGTAAAAGCCAACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAGGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAATATCTCGCCGATATCATCGGTGTCgatg > 1:203363/1‑149 (MQ=255) cGGGTAGTAAGAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCgat < 2:279073/144‑1 (MQ=255) ggTAGTAAAAGTCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGGCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCgatga < 2:4139/144‑1 (MQ=255) tagtaAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGAt < 1:102344/149‑1 (MQ=255) taAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGAtctc > 1:180137/1‑149 (MQ=255) gCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGGAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCgg < 1:252792/149‑1 (MQ=255) cgacgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGc > 2:302972/1‑144 (MQ=255) ccGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATATCGGCTCTATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTa > 2:98018/1‑144 (MQ=255) tGATTGCCGAGGACAATCACTTCGCCAATATCTTCTCTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATCGGTAGTCTTCATAACCATCTCGGCGACATCATCCTGGTCGATGATGCGGATCTCGCCGGGGTTTAGcccc < 2:146581/144‑1 (MQ=255) | AATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCC > CP009273/2835923‑2836205 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |