Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,168,136 | C→G | 22.1% | R90G (CGA→GGA) ‡ | pabC → | 4‑amino‑4‑deoxychorismate lyase |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,168,136 | 0 | C→G | 22.1% | 3.1 ‑2.5 ‑0.0 | 18 | 268 | pabC | 4‑amino‑4‑deoxychorismate lyase |
Reads supporting (aligned to +/- strand): new base (4/0): ref base (5/9): total (9/9) | ||||||||||
Fisher's exact test strand distribution p-value = 8.24e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.31e-01 |
TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG > REL606/1168101‑1168171
|
tACTGAAAGTCGTGATCAGTCGCGGTAGTggcgggc < 1:810967/36‑1 (MQ=255)
tACTGAAAGTCGTGATCAGTCGCGGTAGTggcgggc < 1:3801357/36‑1 (MQ=255)
aaGTCGTGATCAGTCGCGGTAGTGGCGGGGGGgggg > 1:1954491/1‑35 (MQ=37)
aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGgggg > 1:1789389/1‑35 (MQ=255)
aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt < 1:33180/36‑1 (MQ=255)
aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg > 1:189200/1‑34 (MQ=25)
aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg > 1:3499161/1‑34 (MQ=25)
aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa < 1:5016397/36‑1 (MQ=255)
aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGcaca < 1:2215117/36‑1 (MQ=255)
tCAGTCGCGGTAGTGGCGGGCGAGGGTACAGcactt < 1:2727851/36‑3 (MQ=255)
cgcgGTAGTGGCGGGCTGGTGTACAGCACATTGAAc < 1:667143/36‑1 (MQ=25)
cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAAc < 1:5130302/36‑1 (MQ=255)
cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAAc < 1:1737914/36‑1 (MQ=255)
gcgGTAGTGGCGGGGGGGGGGGGAGCACAATGaaaa > 1:3816148/1‑34 (MQ=9)
cgGTAGTGGCGGGCGGGGGGACAGGACAATGAAAAg > 1:261073/1‑36 (MQ=16)
aGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGa < 1:477321/36‑1 (MQ=255)
gTGGCGGGCGTGGGTACAGCCCATTGAACAGCGGAc > 1:1161171/1‑36 (MQ=37)
tGGCGGGCGAGGGTACAGCACATTGATCAGCGGAcc < 1:3865600/36‑1 (MQ=255)
tGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc > 1:5363078/1‑36 (MQ=255)
tGGCGGGCGAGGGTACAGCACATGGAACAGCGGAcc > 1:4334752/1‑36 (MQ=255)
gcgggcgAGGGTACAGCACATTGAACAGCGGACCAg < 1:4865052/36‑1 (MQ=255)
gcgggcgAGGGTACAGCACATTGAACAGCGGAACAg > 1:2175603/1‑36 (MQ=255)
gggcgAGGGTACAGCACATTGAACAGCGGCCCAGcc < 1:3714672/36‑2 (MQ=255)
gggcgAGGGTACAACACATTGAACAGCGGACCAGCa > 1:3108442/1‑36 (MQ=255)
gcgAGGGTACAGCACATTGAACAGCGGACCAGCAAc > 1:549268/1‑36 (MQ=255)
gcgAGGGTACAGCACATTGAACAGCGGACCAGCAAc > 1:5791935/1‑36 (MQ=255)
cgAGGGTACAGCACATTGAACAGCGGACCAGCAAcg > 1:5844205/1‑36 (MQ=255)
|
TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG > REL606/1168101‑1168171
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 7 ≤
ATCG
< 16 ≤
ATCG
< 27 ≤
ATCG
< 37 ≤
ATCG
|