| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 978,122 | C→G | 56.0% | A341A (GCG→GCC) | dcp ← | dipeptidyl carboxypeptidase II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 978,122 | 0 | C | G | 56.0% | 3.1 / 23.7 | 25 | A341A (GCG→GCC) | dcp | dipeptidyl carboxypeptidase II |
| Reads supporting (aligned to +/- strand): ref base C (6/5); new base G (12/2); total (18/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.77e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 | |||||||||||
TTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTA > minE/978057‑978189 | tttAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGc > 1:2251063/1‑68 (MQ=255)tttAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGc > 1:2790798/1‑68 (MQ=255) aTTCGAAATAAGGCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATaa < 1:531140/68‑1 (MQ=255) aTTCGAAATAAGGCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATaa < 1:1675591/68‑1 (MQ=255) aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATaa < 1:1292093/68‑1 (MQ=255) aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATaa < 1:359857/68‑1 (MQ=255) cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATg > 1:1117210/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcg > 1:1575012/1‑68 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:1203146/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:3246045/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2801327/1‑69 (MQ=255) ggCTTGAGCTGCGCTTCATCGAGATCGAATTTCTCACGTCTTACCTGCTCGGCATAAAATGCCCAGTcc > 1:2575637/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) 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) ccTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCg < 1:2996445/69‑1 (MQ=255) ccGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTg < 1:90215/67‑1 (MQ=255) ccGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTa > 1:2957522/1‑69 (MQ=255) ccGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTa > 1:2150737/1‑69 (MQ=255) | TTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTA > minE/978057‑978189 |
| 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: ‑ |