| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,304,649 | C→T | 12.9% | intergenic (+256/+459) | eco → / ← mqo | ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,304,649 | 0 | C | T | 12.9% | 66.0 / 4.0 | 31 | intergenic (+256/+459) | eco/mqo | ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain |
| Reads supporting (aligned to +/- strand): ref base C (15/12); new base T (2/2); total (17/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01 | |||||||||||
AAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGG > NC_000913/2304516‑2304713 | aaGTGCGGCACCCCGTAGGTCAGACAAGGCGGTGACGGCGCATCCGAGATCCACCGCCCGAGCCTGTTGCCTGCTGCGACGCGGGCGGGGCTTATCAGGGCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa < 5:215200/148‑1 (MQ=14)aaGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc > 3:129457/1‑149 (MQ=255) gTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTGATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc > 5:199581/1‑149 (MQ=255) ggCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACGGCAGTGAAGTGCTCCACCACGAAGGACCGATAAGGCGCTTAagcaga > 2:119682/1‑145 (MQ=11) ccacccGAACGTCAGACTAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATc < 5:281080/146‑1 (MQ=255) cACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGACTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAcggcggg > 6:17701/1‑137 (MQ=255) ccccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcaccgcctac > 1:390472/1‑144 (MQ=255) aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATATGACGGTGACGCCGCATCCGTCCTCTCACg > 5:167787/1‑148 (MQ=11) ccAAGGCTGTCACTCCGCATCCGCCATCCACCGCCCGCGCCGGGTGCCTGCTGCGCCGCTGGCGGGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGccc < 5:442164/148‑1 (MQ=11) aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc > 1:211535/1‑109 (MQ=35) aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc < 2:211535/109‑1 (MQ=35) aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATc > 5:508622/1‑132 (MQ=21) aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATc < 6:508622/132‑1 (MQ=34) tggTCACGCCGCATGCGTCAGCCATCGGCCGAGCGGGTTGGTTGCTGCTCCGCTGTCGCGTCTTATCAGGCCCACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg < 3:302334/148‑1 (MQ=11) ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg > 6:516968/1‑105 (MQ=35) ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg < 5:516968/105‑1 (MQ=35) tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt > 1:297636/1‑149 (MQ=33) cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa < 3:421148/102‑1 (MQ=32) cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa > 4:421148/1‑102 (MQ=32) gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc < 4:413485/113‑1 (MQ=32) gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc > 3:413485/1‑113 (MQ=32) cgcCTCCGACATCCAACGCCCGAGCCGGTTGCCTGATTCGACTCTGGAGCGTCTTATCAGGCCTACACCGCTGTGAAGTTCTCCACCCCGTACGTCGTATAAGGCGGTTTCTCCGCATCCGACATCTAACGCCCAAGCCg < 4:424960/140‑1 (MQ=11) cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCAGGCGCGTCTAATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGAAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCTAACGCCCAAGCCg > 3:424960/1‑140 (MQ=255) tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATc > 4:274240/1‑122 (MQ=32) ccGACATCCAACGCCCGTTGAGGTTGCCTGATGCGACGATGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGGGCGCCACCCCGTAGGTCGGAGAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa < 3:316095/149‑1 (MQ=11) gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgc > 2:132423/1‑106 (MQ=17) gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgc < 1:132423/106‑1 (MQ=17) aCATCCCACGCCCGAGCCGGTTGCCTTCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc < 1:119682/149‑1 (MQ=12) tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg > 1:192201/1‑149 (MQ=32) cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGc < 2:125448/92‑1 (MQ=14) cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGc > 1:125448/1‑92 (MQ=14) aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc < 4:114897/106‑1 (MQ=32) aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc > 3:114897/1‑106 (MQ=32) | AAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGG > NC_000913/2304516‑2304713 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |