Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,326,566 | A→T | D67E (GAT→GAA) | mreD ← | cell wall structural complex MreBCD transmembrane component MreD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,326,566 | 0 | A | T | 98.7% | 619.2 / ‑0.0 | 152 | D67E (GAT→GAA) | mreD | cell wall structural complex MreBCD transmembrane component MreD |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base T (101/49); total (102/50) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.51e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01 |
Only 100 of 153 total aligned reads displayed.
ATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAAC > REL606/3326493‑3326638 | aTTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTc > 2:4384868/1‑76 (MQ=255) tttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTcc > 2:4221630/1‑76 (MQ=255) ttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCa > 1:3538800/1‑76 (MQ=255) tcggcgcCACCAGGAACGCAATGAGGCACATCGCAAATCCGGGTCCCCAAGGCGTGGAGCCGCTGATCAGTTCCAg < 2:2797052/73‑1 (MQ=255) tCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAg > 1:3864669/1‑76 (MQ=255) tCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAg > 2:6450068/1‑76 (MQ=255) cAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGt > 2:2343803/1‑76 (MQ=255) aGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTa > 1:5056147/1‑76 (MQ=255) aGCGCCACCAGGAAAGAAAGAATGCTCATCGCCAATGCGCGTACGCCAAGCGGGGAGCCGCTGATCAGTTCCAGTa < 2:5441482/76‑1 (MQ=255) ccccaccaGGTAAGCAATGATGCTCACCGCCAAGCCGCGCACGCCAAGCGGCGAGCCGCTGATCAGTTCCAGTATg < 1:3214774/74‑1 (MQ=255) gcCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATgg > 1:4177636/1‑76 (MQ=255) cccacaaGGTAAGCAAGGATGCCCACCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATgg < 1:1578188/75‑1 (MQ=255) ccaccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGc > 2:4425415/1‑76 (MQ=255) cacgaGGAAAGAAATGATCCTCACCGCAAATACGGGTACCCAAGGCGGGGAGCCGCTGATCAGTTCCAGTATGGCa < 2:4512353/76‑1 (MQ=255) caccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCa > 2:4707386/1‑76 (MQ=255) caccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCa > 1:5673387/1‑76 (MQ=255) caccaGGAAAGCAATGTTGCCCATCGCCAATCCGCGTACCCAAAGCGGGGAGCCGCTGATCAGTTCCAGTATGGCa < 2:857731/76‑1 (MQ=255) accaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCAc > 2:2434254/1‑76 (MQ=255) accaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCAc > 1:6200149/1‑76 (MQ=255) ccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCAcc > 2:6941105/1‑76 (MQ=255) aGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCa > 1:5013436/1‑76 (MQ=255) gTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATc > 1:1881698/1‑76 (MQ=255) tAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGAACCCATca > 2:2185357/1‑76 (MQ=255) tAAGAAATGAGCCGCACCCCCAATGCGCGTACGCCAGGCGGCGAGCCGCTGATCAGTTCCAGTATGGCACCCATca < 1:6899585/76‑1 (MQ=255) gaaGAAATGATCCTCATCGCCAATCCCCGTACGCCAAGCGGGGAGCCGCTGATCAGTTCCAGTATGGCACCCATca < 2:6904446/75‑1 (MQ=255) cAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACaaa > 1:4543525/1‑76 (MQ=255) atggatgCTCAGCGCCAATCCGCGTCCGCCAGGCGTGGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACaaaa < 2:6255459/73‑1 (MQ=255) aaTGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACaaaa > 2:1815449/1‑76 (MQ=255) atgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGGGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAAc > 2:6530002/1‑76 (MQ=255) atgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAAc > 1:6036042/1‑76 (MQ=255) tgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACGCATCACAAAAcc > 1:1745030/1‑76 (MQ=255) tgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAAcc > 1:3867841/1‑76 (MQ=255) gCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCtgtg > 2:2340783/1‑76 (MQ=255) cTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGc > 1:4928700/1‑76 (MQ=255) tCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGcc > 2:5295158/1‑76 (MQ=255) aTCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCca > 2:1726224/1‑76 (MQ=255) aTCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCca > 1:4128792/1‑76 (MQ=255) tCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcac > 2:3590340/1‑76 (MQ=255) tCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcac > 2:623419/1‑76 (MQ=255) cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcaca > 1:2527186/1‑76 (MQ=255) cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcaca > 2:2999454/1‑76 (MQ=255) ccaaaTCCGCGTACCCAAGGGGGCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAt < 2:4250300/73‑1 (MQ=255) caaaTGCGCGTCCCCAAAGCGGCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAtt < 2:3956066/74‑1 (MQ=255) cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt > 2:6559136/1‑76 (MQ=255) cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt > 2:2055312/1‑76 (MQ=255) cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt > 2:2424531/1‑76 (MQ=255) cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt > 1:4862987/1‑76 (MQ=255) aaTCCGCGGACGCCAGGGGTGGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa < 1:3372694/76‑1 (MQ=255) aaTCCGCGCACGCAAGGGGTGGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa < 1:3540784/76‑1 (MQ=255) aaTACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa > 2:578616/1‑76 (MQ=255) cgcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTAcgcg > 2:6942559/1‑76 (MQ=255) cgcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTAcgcg > 2:2995297/1‑76 (MQ=255) gcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGa > 2:1237633/1‑76 (MQ=255) gcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGCATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGa > 1:742285/1‑76 (MQ=255) gTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATg > 2:55982/1‑76 (MQ=255) gTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATg > 2:102439/1‑76 (MQ=255) tACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGa > 2:4606898/1‑76 (MQ=255) tACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGa > 1:6595556/1‑76 (MQ=255) gCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGc > 1:1008598/1‑76 (MQ=255) gCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGc > 1:6105901/1‑76 (MQ=255) ccAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCa > 2:4991549/1‑76 (MQ=255) cAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCaa > 1:3495850/1‑76 (MQ=255) aaGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAg > 2:2661428/1‑76 (MQ=255) aaGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAg > 1:237431/1‑76 (MQ=255) gCGTGGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGc < 2:4658567/76‑1 (MQ=255) ggggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa < 2:6743943/73‑1 (MQ=255) ggggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa < 2:3229007/73‑1 (MQ=255) gcggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa < 1:1213936/73‑1 (MQ=255) gggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAg < 2:318313/74‑1 (MQ=255) ggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAgg < 1:2138038/75‑1 (MQ=255) cGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAgg > 2:5302780/1‑76 (MQ=255) gAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGa > 1:1308729/1‑76 (MQ=255) gAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGa > 1:2237889/1‑76 (MQ=255) gAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGa > 2:6890122/1‑76 (MQ=255) aGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGAt < 1:5154273/76‑1 (MQ=255) aGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGAt < 1:4654358/76‑1 (MQ=255) aGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGAt < 1:2340908/76‑1 (MQ=255) gCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATc < 1:1698556/76‑1 (MQ=255) ccGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATcc < 1:6450068/76‑1 (MQ=255) ccGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATcc > 1:6817648/1‑76 (MQ=255) cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa > 2:5731917/1‑76 (MQ=255) cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa < 2:5089947/76‑1 (MQ=255) cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa < 2:158732/76‑1 (MQ=255) cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa > 2:3965859/1‑76 (MQ=255) gCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCaa > 2:1390448/1‑76 (MQ=255) cTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAAt > 2:5981403/1‑76 (MQ=255) cTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAAt > 2:4635371/1‑76 (MQ=255) cTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAAt > 1:3744409/1‑76 (MQ=255) tGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATa > 2:4788580/1‑76 (MQ=255) tGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATa > 1:1573432/1‑76 (MQ=255) tGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATa < 1:777985/76‑1 (MQ=255) gATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATac > 1:624395/1‑76 (MQ=255) gATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGACCCAATac > 1:6895011/1‑76 (MQ=255) aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca < 1:3739052/76‑1 (MQ=255) aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca < 1:573065/76‑1 (MQ=255) aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca < 1:1772623/76‑1 (MQ=255) aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca < 1:6886139/76‑1 (MQ=255) tCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATacaa < 1:192531/76‑1 (MQ=255) cAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATacaac < 2:1035005/76‑1 (MQ=255) cAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATacaac < 1:1059633/76‑1 (MQ=255) | ATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAAC > REL606/3326493‑3326638 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |