Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 2,086,059 | T→C | 20.4% | R389R (CGT→CGC) | sbcB → | exonuclease I |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,086,059 | 0 | T | C | 20.4% | 74.8 / 16.5 | 40 | R389R (CGT→CGC) | sbcB | exonuclease I |
Reads supporting (aligned to +/- strand): ref base T (15/16); new base C (2/6); total (18/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01 |
CTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGG > W3110S.gb/2085991‑2086125 | cTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGCGCCGCGCaa > 1:1054946/1‑71 (MQ=255) cTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGCaa > 1:1653548/1‑71 (MQ=255) ataACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTa > 1:4060979/1‑71 (MQ=255) aCGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAAttt > 1:2629364/1‑67 (MQ=255) aCGGCTCTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAAtt > 1:4784593/1‑66 (MQ=255) aCGGCTCTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTAACg > 1:5044829/1‑71 (MQ=255) cGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGCAATTTACCgg < 1:5114616/71‑1 (MQ=255) ttttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTg > 1:3303012/1‑71 (MQ=255) tttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGCAATTTACCGGCACTg < 1:3848385/70‑1 (MQ=255) tGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCa < 1:4099401/71‑1 (MQ=255) gACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTAc > 1:1369875/1‑55 (MQ=255) aCGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTAcc > 1:758223/1‑55 (MQ=255) gCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACttt < 1:4123677/71‑1 (MQ=255) cAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGCAATTTACCGGCACTGGATATCACtttt < 1:3181341/71‑1 (MQ=255) cGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGGAATTTACCggg > 1:3531393/1‑49 (MQ=255) gTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGa > 1:4822458/1‑70 (MQ=255) gTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACAGGCACTGGata > 1:2436755/1‑57 (MQ=255) tGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATa < 1:444716/71‑1 (MQ=255) tGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATa < 1:161570/71‑1 (MQ=255) gcagcaATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTttgtt < 1:2885204/66‑1 (MQ=255) aTGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGcc > 1:2850889/1‑41 (MQ=255) aTGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACt < 1:4943795/54‑1 (MQ=255) aTGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAt > 1:3255199/1‑71 (MQ=255) gAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAtt > 1:1204536/1‑70 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTg < 1:1698855/70‑1 (MQ=255) aaTTGTGCTGGAAACCGAGCCGCGTAATTTACCGGc > 1:2191587/1‑36 (MQ=255) ggAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGa < 1:1148057/48‑1 (MQ=255) ggAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGa < 1:4438880/48‑1 (MQ=255) ggAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGa < 1:3473181/48‑1 (MQ=255) ggAAACCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATa < 1:1890840/50‑1 (MQ=255) ggAAACCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATa < 1:132481/50‑1 (MQ=255) ggAAACCGAGCC‑CGTAATTTACCGGCACTGGATATCACTTTTGTTGa < 1:1626728/47‑1 (MQ=38) aaaCCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTc < 1:2761007/71‑1 (MQ=255) aaaCCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTc < 1:2827138/71‑1 (MQ=255) ccGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGaaa < 1:5077884/57‑1 (MQ=255) ccGAGCCGCGTAATTTACCGGCACTGGATAACACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAAt < 1:1438721/71‑1 (MQ=255) gCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATc > 1:573902/1‑71 (MQ=255) ccGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAtt > 1:2786351/1‑48 (MQ=255) tccgTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCgg < 1:3701986/69‑1 (MQ=255) tccgTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCgg < 1:2308862/69‑1 (MQ=255) | CTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGG > W3110S.gb/2085991‑2086125 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |