| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,785,101 | A→C | 50.0% | A80A (GCT→GCG) | pepE ← | (alpha)‑aspartyl dipeptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,785,101 | 0 | A | C | 50.0% | 0.2 / 42.7 | 40 | A80A (GCT→GCG) | pepE | (alpha)‑aspartyl dipeptidase |
| Reads supporting (aligned to +/- strand): ref base A (11/9); new base C (13/7); total (24/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 | |||||||||||
CCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACA > minE/2785038‑2785158 | cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:1107389/69‑1 (MQ=37)cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:385098/69‑1 (MQ=37)cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:3171608/69‑1 (MQ=37)cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:299822/69‑1 (MQ=37)cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:2910625/69‑1 (MQ=37)cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:1124274/69‑1 (MQ=37)cccACGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTACCGCCGCCGACAATCACAATTTCCGCA‑tt < 1:1566189/69‑1 (MQ=37) tttCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACggg > 1:3082237/1‑69 (MQ=255) tCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACggg < 1:3175265/67‑1 (MQ=255) tCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACggg < 1:7653/67‑1 (MQ=255) ggAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGc < 1:1426796/69‑1 (MQ=255) gTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTTCAATCGCGGCAACGGGATCGACAACGCTATg < 1:1931485/69‑1 (MQ=255) gTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTTCAATCGCGGCAACGGGATCGACAACGCTATg < 1:2660407/69‑1 (MQ=255) ttCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATa > 1:2966455/1‑69 (MQ=255) cccGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATAcc > 1:991799/1‑69 (MQ=255) cccGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATAcc > 1:432350/1‑69 (MQ=255) cccGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATAcc > 1:2677328/1‑69 (MQ=255) cccGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATAc > 1:211033/1‑68 (MQ=255) gccgccgACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATTCCGGt < 1:2052653/69‑1 (MQ=255) gccgccgACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGt < 1:1107066/69‑1 (MQ=255) gccgccgACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGt < 1:844598/69‑1 (MQ=255) gccgccgACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACCCTATGAATACCGGt < 1:3110152/69‑1 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATc > 1:3192768/1‑41 (MQ=37) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:3168205/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:503066/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:389711/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:3325535/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:1949231/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:2508907/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:3126565/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:3075907/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:2450975/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:2839614/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAACGGGATCGACAACg > 1:2733842/1‑48 (MQ=255) aCAATCACAATTTCCGCA‑TTTTCAATCGCGGCAAAGGGATCGACAACg > 1:2431656/1‑48 (MQ=255) cACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATAcc > 1:3238512/1‑54 (MQ=255) cACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAaca > 1:1306543/1‑67 (MQ=255) cACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAaca > 1:2072882/1‑67 (MQ=255) cACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAaca > 1:2068848/1‑67 (MQ=255) cACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAaca > 1:1790843/1‑67 (MQ=255) | CCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCA‑TTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACA > minE/2785038‑2785158 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |