| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 978,092 | A→G | 35.8% | L351L (CTT→CTC) | dcp ← | dipeptidyl carboxypeptidase II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 978,092 | 0 | A | G | 35.8% | 21.3 / 20.3 | 39 | L351L (CTT→CTC) | dcp | dipeptidyl carboxypeptidase II |
| Reads supporting (aligned to +/- strand): ref base A (5/20); new base G (12/2); total (17/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 | |||||||||||
TACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCC > minE/978035‑978140 | tACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAAttt > 1:1009210/1‑69 (MQ=255) tCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCccgc < 1:1132394/69‑1 (MQ=255) tCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCccgc < 1:3276647/69‑1 (MQ=255) tCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCccgc < 1:3255404/69‑1 (MQ=255) tCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCccgc < 1:1355139/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2567613/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:933006/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:3294473/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2985986/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2719899/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2506737/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2320811/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:1956657/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:1852563/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:1322632/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:1691452/69‑1 (MQ=255) acCGTGTTTAATTCAAAATATGGCTTGAGCTACGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2645641/69‑1 (MQ=255) tttAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGc > 1:2251063/1‑68 (MQ=255) tttAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGc > 1:2790798/1‑68 (MQ=255) tAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTTCCGCCGTACCTg > 1:2611139/1‑60 (MQ=255) aTTCGAAATAAGGCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATaa < 1:1675591/68‑1 (MQ=255) aTTCGAAATAAGGCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATaa < 1:531140/68‑1 (MQ=255) aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCccgc < 1:2887541/50‑1 (MQ=255) aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATaa < 1:1292093/68‑1 (MQ=255) aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATaa < 1:359857/68‑1 (MQ=255) cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATg > 1:1117210/1‑69 (MQ=255) tatGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:667854/51‑1 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcg > 1:1575012/1‑68 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1554795/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2575637/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2801327/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1607320/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:3246045/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1249564/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1203146/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:3397235/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1633559/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:534255/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:90394/1‑69 (MQ=255) | TACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCC > minE/978035‑978140 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |