Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 22,700 | T→G | 17.7% | Q150H (CAA→CAC) | dck ← | deoxyadenosine/deoxycytidine kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 22,700 | 0 | T | G | 17.7% | 53.5 / 2.1 | 37 | Q150H (CAA→CAC) | dck | deoxyadenosine/deoxycytidine kinase |
Reads supporting (aligned to +/- strand): ref base T (12/18); new base G (5/2); total (17/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.12e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTCAATGCGAAGTTTTAAAACAGGGCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACATTGT > BS168_reference_genome/22583‑22834 | ttCAATGCGAAGTTTAAAAACATGGCAGGAATTAAAACGGCTTATCCAGTTCTGGTAGCGTGTATGCATTTCTTCCAAATAAGACCGGCGTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATg < 5:83879/139‑1 (MQ=255) cAATGCGAAGTTTTAAAACAGGGCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTAGTTGAATACGGGTTTAAAtttt > 3:75414/1‑139 (MQ=255) aTGCGAAGTTTTAAAACAGGGCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTGACTCCGGTTTAAAATGTtct > 2:38379/1‑139 (MQ=255) aGTTTTAAAACAGGGCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTTTCTGAAGCTCCATCTCACGGCCGCGTTGTTAAATACGGGTTAAAATGTTATCCAGa > 1:49373/1‑138 (MQ=255) gTTTTAAACCCGGGCAGGCATAAAACCCGCTTATCCAGTACTCGTCGCGTGAATGCATTACTACCCAAAAACACCGGCTGGGCTGAAGCGCCATCACACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGAtc < 1:308689/139‑1 (MQ=255) gTTTTAAAACAGGGCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGAtc < 7:171937/139‑1 (MQ=255) gggCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTTAAATGTTGTCCAGATATCCTTGTaaa > 5:347074/1‑138 (MQ=255) caTTAAAACGGCTTATCCAGTTCTCGTAGCGTGTAGCCATTTCATCCCAATAACACCGGCTTGTCTGACGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCat < 4:314344/139‑2 (MQ=255) caTTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGGTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCaa < 1:336901/139‑1 (MQ=255) aaaaCCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAacac > 7:211124/1‑138 (MQ=255) aaaaCCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGGTGTTCAATCCGGGTTAAAATGTTCTCCAGATGTCCATCTAAATAGATCAaccc > 6:132893/1‑136 (MQ=255) gtgtATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTGTCCAGATCTCCTTCTTAATATATCAACACAACTGGATGTGGGAAATACGGTGTCAt > 1:80574/1‑139 (MQ=255) gtgtATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGGTGTTCAATCCGGTTTTAAATGTTATGCAGCTCTCCTTCTAAATATATAAACACATCTGGGTGTGGGTAATAAGGTGTGAt > 1:234874/1‑139 (MQ=255) gtgtATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGGTGTTCAATCCGGGTTTAAATGTTGTGCAGCTATACGTGTTAATATATGAAAAAAACTGGATGTGGGGAATAAGGTGTCAt > 2:81940/1‑139 (MQ=255) tATGCATTTCATCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGGTTTGAAAACCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCAGGGTGTGGGTAAATCGGTGTGAtata > 5:308714/1‑136 (MQ=255) aTTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTg < 4:16265/139‑1 (MQ=255) ccAATAAGACCGGCGTTGCTGAAGCGCCATTTCACGGCCGCGGTGTTCAATCCGGTTTAAAATGTTCTCCAGCTCTCCATCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCAATACCATTGCTTCGaaa < 2:61076/139‑1 (MQ=255) cGGCTTGTCTGAAGCTCCCTCTCACGGCCGCGTTGTTCCATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTgtgt < 2:49373/139‑1 (MQ=255) cGGCTTGGCTGAAGCTCCATCTCACGGCCGCGTAGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTgtgt < 4:75414/139‑1 (MQ=255) gCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCaa < 8:44885/81‑1 (MQ=255) gCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCaa > 7:44885/1‑81 (MQ=255) tCTGAAGCTCCAGCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCtt < 4:95318/139‑1 (MQ=255) aGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAAt < 3:253900/139‑1 (MQ=255) gCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATc < 5:111704/139‑1 (MQ=255) gCTCCATCTCACGGCCGCGGTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATc < 8:211124/139‑1 (MQ=255) cATCTCACGGCCGCGTTGTTGAATCCGGTTTAAAATGTTCTCCCGATCTCCTTCGAAATAGATCAACACATCTGGATGTGGGAGAAACGGTGTCATTACCATTGCGTCGCAAAGGCGTGTGTAAATGTTATAATCAAcc > 5:293436/1‑139 (MQ=255) cATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAAcc > 1:204869/1‑139 (MQ=255) cATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAAc > 3:326055/1‑138 (MQ=255) cATCTCACGGCCGCGGTGTTGAATCCGGTTTAAAATGTTGTACCGATCTCCTTCTAAATATATAAACACAAATGGTTTTGGGGAATAAGGTGTCATTAACATTGCTTGGTAAAGGGTTGTGTATGTCTTATAATGaaa > 6:111704/1‑137 (MQ=255) cGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACa < 1:81940/139‑1 (MQ=255) cGGCCGCGTTGTTCAATACGGTTTTAAATGTTGTCCAGATCTCCTTATAAATAGATAAACACAAATGGATGTGGGTAATACGGTGTCATTAACATTTCTTCGAAAAGGCTTGTGTTAGTCTTATTATAAACCTTTGACa > 6:366980/1‑139 (MQ=255) cGGCCGCGTTGGTCAATCCGGTTTAAAATGTTATCCCGATCTCCTTCTAAATAGATCAACACAAATGGATGTGGGAAATACGGTGTCATTACCATTGCGTCGAAAAGGCATGTGTAAGTCTTTTAATCAACCTTTGACa > 6:329611/1‑139 (MQ=255) gCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACAtt < 6:308714/139‑1 (MQ=255) gCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACAtt < 5:132893/139‑1 (MQ=255) cgcgTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACATTGt < 4:326055/139‑1 (MQ=255) cgcgTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACATTGt < 5:366980/139‑1 (MQ=255) cgcgTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACATTGt < 2:80574/139‑1 (MQ=255) | TTCAATGCGAAGTTTTAAAACAGGGCAGGCATTAAAACCGCTTATCCAGTTCTCGTAGCGTGTATGCATTTCTTCCCAATAAGACCGGCTTGTCTGAAGCTCCATCTCACGGCCGCGTTGTTCAATCCGGTTTAAAATGTTCTCCAGATCTCCTTCTAAATAGATCAACACATCTGGATGTGGGAAATACGGTGTCATTACCATTGCTTCGAAAAGGCTTGTGTAAGTCTTATAATCAACCTTTGACATTGT > BS168_reference_genome/22583‑22834 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |