Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 978,110 | C→T | 42.4% | R345R (CGG→CGA) ‡ | dcp ← | dipeptidyl carboxypeptidase II |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 978,110 | 0 | C | T | 42.4% | 7.0 / 21.2 | 33 | R345R (CGG→CGA) ‡ | dcp | dipeptidyl carboxypeptidase II |
Reads supporting (aligned to +/- strand): ref base C (4/15); new base T (12/2); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.61e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGAT > minE/978042‑978178 | 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: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:2567613/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:359857/68‑1 (MQ=255) aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATaa < 1:1292093/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:1249564/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:90394/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:534255/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:3397235/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1203146/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2575637/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:3246045/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1633559/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1607320/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2801327/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1554795/1‑69 (MQ=255) ccgTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGAt < 1:2039223/69‑1 (MQ=255) | TCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGAT > minE/978042‑978178 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |