| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AE015451 | 3,925,054 | G→T | 11.8% | V479V (GTC→GTA) | peaA ← | quinohaemoprotein amine dehydrogenase, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AE015451 | 3,925,054 | 0 | G | T | 11.8% | 101.4 / 7.6 | 34 | V479V (GTC→GTA) | peaA | quinohaemoprotein amine dehydrogenase, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (16/14); new base T (2/2); total (18/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01 | |||||||||||
GAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAACCGCCAGGTTCACGCCCTT > AE015451/3924841‑3925284 | gAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCAc < 1:578750/236‑1 (MQ=255) cGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGt < 1:306418/236‑1 (MQ=255) ccTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACgg < 1:414541/235‑1 (MQ=255) gttgttCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCa < 1:604971/236‑1 (MQ=255) gttgttCCAGCGCGGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTTACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCa < 1:14585/236‑1 (MQ=255) cGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCt > 1:89858/1‑211 (MQ=255) cGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTCTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAg < 1:418230/235‑1 (MQ=255) cGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGc > 2:595185/1‑204 (MQ=255) cGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCTGGTAGCCGa > 1:434546/1‑222 (MQ=255) aTGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTTACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCAt > 1:507985/1‑235 (MQ=255) tgCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGc > 2:520313/1‑219 (MQ=255) ccACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTc < 1:332315/236‑1 (MQ=255) aGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGcc > 2:31309/1‑182 (MQ=255) aGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCa < 1:553259/235‑1 (MQ=255) tACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAg < 1:520314/236‑1 (MQ=255) aCCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGa > 2:667408/1‑234 (MQ=255) ggCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGACCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAgg > 2:350552/1‑194 (MQ=255) gCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCAc > 1:539101/1‑235 (MQ=255) ttgGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTtct > 1:127285/1‑235 (MQ=255) gTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGACACGCCCCTGCACTTTCGGCACGGAGGCa > 1:335928/1‑205 (MQ=255) aCCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCtt > 2:170870/1‑155 (MQ=255) aCCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCt > 1:639732/1‑236 (MQ=255) cTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGc > 2:3621/1‑148 (MQ=255) ggCACGAACACACCGTCGGCCCGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCCCTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGgag > 1:318838/1‑213 (MQ=255) aCGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCa < 2:307672/175‑1 (MQ=255) gAACACACCGTCGGCCTGCATCGGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTGCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACttt < 2:37634/168‑1 (MQ=255) gTCGGCCTGCATCTGCCCGGCGAACTTTACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCa < 1:129946/86‑1 (MQ=255) gTCGGCCTGCATCTGCCCGGCGAACTTTACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCa > 2:129946/1‑86 (MQ=255) cTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAAcc < 1:667409/235‑1 (MQ=255) cATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATgcgc < 2:611356/181‑1 (MQ=255) aTCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGCGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAACCGc > 2:350676/1‑233 (MQ=255) aTCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAACCGcccg > 1:85665/1‑234 (MQ=255) aaCTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAACCGCCAGGTTCACGCCCtt < 1:721538/236‑1 (MQ=255) aattGACGTCTTCATCATCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAACCGCCAGGTTCACGCCCt < 2:312238/232‑1 (MQ=255) | GAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTGGCCACCGTCGGCCAGCGTGGCGATGACCTTCAGGTTACCGGCGTTGTTGGTCATCATCTTGCGCTCAGGGTTGGGCCCGGCACCGCCTGGCACGAACACACCGTCGGCCTGCATCTGCCCGGCGAACTTGACGTCTTCATCCTCGACAGCACGCTCGTTGAACGGCTCGACCTTCCAGCTGGCCGGCAGGTAGCCGATACGCAGCGGCTGGCCATTGGCATCCTTGCCCCAGGCTTCGGCCTCGAAACGCCCCTGCACTTTCGGCACGGAGGCACCGTTCTCGCCAATGCGCGCGATGGAGAACGCCGGCACGACCTTCACTTCGTCGACCTTGTCGTAAACCGCCAGGTTCACGCCCTT > AE015451/3924841‑3925284 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |