Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,689,625 | C→T | A213T (GCC→ACC) | bcsZ ← | endo‑1,4‑D‑glucanase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,689,625 | 0 | C | T | 88.0% | 54.7 / 1.4 | 25 | A213T (GCC→ACC) | bcsZ | endo‑1,4‑D‑glucanase |
Reads supporting (aligned to +/- strand): ref base C (1/2); new base T (12/10); total (13/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 |
TCATGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCC > NC_000913/3689502‑3689735 | tcatGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCaa > 1:479781/1‑136 (MQ=255) tcatGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCaa < 2:479781/136‑1 (MQ=255) aTGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGttt > 1:530144/1‑139 (MQ=255) gTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTcg > 1:917810/1‑139 (MQ=255) aaaCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCGGCCAGCCTTTGTCTTTCTCATAGCGCAACTAGTCTGGCGAAAATCATTTCGGGGTGGTTTCCAGCAATAATTGTTGATTGGTTTcgcg < 2:1042847/139‑1 (MQ=255) aCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGGCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCa < 2:548247/139‑1 (MQ=255) gcgGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGc > 2:575439/1‑139 (MQ=255) ctgctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGgcgcgc > 2:590651/1‑139 (MQ=255) aTCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGc < 1:590651/139‑1 (MQ=255) cAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt < 1:447050/101‑1 (MQ=255) cAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt > 2:447050/1‑101 (MQ=255) tGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTg > 1:378447/1‑139 (MQ=255) gTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGa < 2:462223/139‑1 (MQ=255) ttttttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGaa > 2:710031/1‑139 (MQ=255) cAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcggcg < 2:283980/139‑1 (MQ=255) ccTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCa > 1:741790/1‑139 (MQ=255) aGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCgg < 2:480928/95‑1 (MQ=255) aGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCgg > 1:480928/1‑95 (MQ=255) aCCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc < 2:378447/139‑1 (MQ=255) aCCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc > 2:554348/1‑139 (MQ=255) aCCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc < 2:917810/139‑1 (MQ=255) gCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt < 1:900449/42‑1 (MQ=255) gCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt > 2:900449/1‑42 (MQ=255) tttCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCa > 1:1035071/1‑44 (MQ=255) tttCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCa < 2:1035071/44‑1 (MQ=255) | TCATGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCC > NC_000913/3689502‑3689735 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |