| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 978,112 | G→T | 42.4% | R345R (CGG→AGG) ‡ | dcp ← | dipeptidyl carboxypeptidase II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 978,112 | 0 | G | T | 42.4% | 15.3 / 26.8 | 33 | R345R (CGG→AGG) ‡ | dcp | dipeptidyl carboxypeptidase II |
| Reads supporting (aligned to +/- strand): ref base G (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 | |||||||||||
ACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAA > minE/978051‑978180 | acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:933006/69‑1 (MQ=255)acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2985986/69‑1 (MQ=255)acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:3294473/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:1691452/69‑1 (MQ=255)acCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:1322632/69‑1 (MQ=255)acCGTGTTTAATTCAAAATATGGCTTGAGCTACGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGt < 1:2645641/69‑1 (MQ=255) tttAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGc > 1:2790798/1‑68 (MQ=255) tttAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGc > 1:2251063/1‑68 (MQ=255) tAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTTCCGCCGTACCTg > 1:2611139/1‑60 (MQ=255) aTTCGAAATAAGGCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATaa < 1:531140/68‑1 (MQ=255) aTTCGAAATAAGGCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATaa < 1:1675591/68‑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:3397235/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:534255/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:90394/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:3246045/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2575637/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2801327/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1633559/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1607320/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1554795/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1249564/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1203146/1‑69 (MQ=255) ccgTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGAt < 1:2039223/69‑1 (MQ=255) cgTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATa < 1:2386968/69‑1 (MQ=255) gTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCt > 1:930553/1‑43 (MQ=255) gTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATaa > 1:2922242/1‑69 (MQ=255) gTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATaa > 1:3363488/1‑69 (MQ=255) | ACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAA > minE/978051‑978180 |
| 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: ‑ |