Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,999,913 | A→G | Q406R (CAG→CGG) | ubiB → | putative ubiquinone biosynthesis protein UbiB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,999,913 | 0 | A | G | 93.8% | 175.3 / 1.6 | 48 | Q406R (CAG→CGG) | ubiB | putative ubiquinone biosynthesis protein UbiB |
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (21/24); total (24/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.34e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 |
TACTGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACAGCG > REL606/3999844‑3999983 | tACTGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAAc < 1:992945/76‑1 (MQ=255) tACTGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGCGTTACTCCCGAAAAc > 1:717584/1‑76 (MQ=255) aCTGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAAcc > 1:1096843/1‑76 (MQ=255) aCTGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAAcc > 2:1046075/1‑76 (MQ=255) tGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCt < 1:1124157/76‑1 (MQ=255) gTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTg < 1:1188822/76‑1 (MQ=255) ttAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGc < 2:1100882/76‑1 (MQ=255) ttAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGc > 2:59638/1‑76 (MQ=255) aaTCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGctct > 1:62654/1‑76 (MQ=255) tCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTAt > 2:812066/1‑76 (MQ=255) cTGTTTAATACGGCGCGTCGCTTCAATAGGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATg < 1:658344/76‑1 (MQ=255) tttAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCg > 2:787868/1‑76 (MQ=255) aaTACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAg < 2:281110/76‑1 (MQ=255) tACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAggg < 1:1266732/76‑1 (MQ=255) tACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAggg > 1:287995/1‑76 (MQ=255) cGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGt > 2:1229160/1‑76 (MQ=255) gcgcGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAg < 2:1002140/76‑1 (MQ=255) cgcgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAgg < 1:265822/76‑1 (MQ=255) cgcgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAgg > 2:423771/1‑76 (MQ=255) cgcgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAgg < 2:560152/76‑1 (MQ=255) cgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGAc < 1:1240039/76‑1 (MQ=255) cgTCGCTTAAATAGGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGAc < 1:656882/76‑1 (MQ=255) ttCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGc < 2:81465/76‑1 (MQ=255) atatGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTa > 2:502455/1‑76 (MQ=255) atatGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTAGAAGGGGTAGGACGCCAGCTTTa > 2:629767/1‑76 (MQ=255) atGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATc > 2:70879/1‑76 (MQ=255) ggAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCg < 1:1208010/76‑1 (MQ=255) aaGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCa < 1:970034/76‑1 (MQ=255) tGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCt > 2:416193/1‑76 (MQ=255) cAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCg > 1:1212084/1‑76 (MQ=255) cGCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGAttt > 2:427468/1‑76 (MQ=255) gCAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTa < 1:6388/76‑1 (MQ=255) cAACTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTAt > 2:391775/1‑76 (MQ=255) aaCTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATg > 1:28313/1‑76 (MQ=255) aaCTGGTGTTACTCCAGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATg > 2:621467/1‑76 (MQ=255) aCTGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATgg < 1:323880/76‑1 (MQ=255) ttGGTGTTGCTCCGGAAGGCCCTGCTCTATGTCGAGGGGGTAGGGGGCCAGCTTTATCCGCAGCTCGATTTATGGa < 1:599042/75‑1 (MQ=255) tGGTGTTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGaa < 1:519087/76‑1 (MQ=255) gtgtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGaaaa < 2:1190652/76‑1 (MQ=255) tgtTACTCGAGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAc < 1:231887/76‑1 (MQ=255) tgtTACTCCGGAAATCCCTGCGCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAc < 1:719469/76‑1 (MQ=255) tgtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGCTTTATGGAAAAc > 1:286918/1‑76 (MQ=255) tgtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAc < 1:183512/76‑1 (MQ=255) tgtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAc < 1:958527/76‑1 (MQ=255) tgtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAc < 1:895396/76‑1 (MQ=255) tgtTACTCCAGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAc > 2:48335/1‑76 (MQ=255) gtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAAcn < 1:1097769/76‑2 (MQ=255) gtTACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACa < 1:1015954/76‑1 (MQ=255) ttACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACAg > 2:1161918/1‑76 (MQ=255) tACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACAGc > 1:347761/1‑76 (MQ=255) tACTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACAGc > 1:431224/1‑76 (MQ=255) aCTCCGGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACAGCg > 2:1008390/1‑76 (MQ=255) | TACTGTTAAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTATGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAGCTCGATTTATGGAAAACAGCG > REL606/3999844‑3999983 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |