Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,853,162 | T→A | 21.8% | I169N (ATT→AAT) | alkB → | alpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,853,162 | 0 | T | A | 21.8% | 39.7 / 8.8 | 23 | I169N (ATT→AAT) | alkB | alpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase |
Reads supporting (aligned to +/- strand): ref base T (11/7); new base A (2/3); total (13/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.18e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01 |
GTCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTGAACAGGG > NC_002947/3853014‑3853310 | gTCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAtt < 2:189096/150‑1 (MQ=255) tCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTg > 2:221741/1‑150 (MQ=255) aaCCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTgcg > 2:142861/1‑150 (MQ=255) aTTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa > 1:395591/1‑150 (MQ=255) tACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGc < 1:143970/150‑1 (MQ=255) aCCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGcc < 2:574859/150‑1 (MQ=255) cTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTc > 1:144909/1‑150 (MQ=255) gAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTg < 1:15329/150‑1 (MQ=255) gCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTggg < 1:142861/150‑1 (MQ=255) ccTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGACTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTgggg > 2:390149/1‑150 (MQ=255) aCCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGAcc < 1:221741/150‑1 (MQ=255) ggATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCa > 2:477469/1‑150 (MQ=255) gTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGc > 1:40423/1‑150 (MQ=255) gATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGAt > 2:406339/1‑150 (MQ=255) gtgtTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGaa < 2:395591/150‑1 (MQ=255) cTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcggcgg > 2:243848/1‑150 (MQ=255) ggcggcCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAAc < 2:40423/150‑1 (MQ=255) cggcCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAAAACCAATGCGCTTCCCAAACGGCCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCt < 1:381650/150‑1 (MQ=255) aCAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGCCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCgg > 2:500903/1‑150 (MQ=255) tggCTCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTg > 1:8160/2‑150 (MQ=255) tCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTg > 1:587178/1‑150 (MQ=255) tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTCGTCTGGGGCGGCGAAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCTGGTGTGCATCCGCGCCTggg < 2:469109/108‑1 (MQ=255) tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTCGTCTGGGGCGGCGAAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCTGGTGTGCATCCGCGCCTggg > 1:469109/1‑108 (MQ=255) aTTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTGAACAggg < 2:144909/150‑1 (MQ=255) | GTCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTGAACAGGG > NC_002947/3853014‑3853310 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |