Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,585,813 | 0 | C→A | 7.7% | 0.2 ‑5.9 ‑0.0 | 57 | 752 | yfhR | predicted peptidase |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/5): ref base (38/14): total (38/19) | ||||||||||
Fisher's exact test strand distribution p-value = 2.78e-03 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.47e-01 |
CCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATG > REL606/2585778‑2585848
|
cccATGGCAGCACAGCGAAAAGTTGTATAGCCTGGc > 1:4147985/1‑36 (MQ=255)
ccATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCa < 1:4147770/36‑1 (MQ=255)
cATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaa > 1:1916426/1‑36 (MQ=255)
cATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaa > 1:1002072/1‑36 (MQ=255)
tGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaaaa < 1:5345658/36‑1 (MQ=255)
ggCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAag > 1:381347/1‑36 (MQ=255)
ggCAGCACAGCGAAAAGCTGTATAGCCTGGCAAAag > 1:4168452/1‑36 (MQ=255)
gcagcaCAGCGAAAAGTTGTATAGCCTGGCAAAaga > 1:5496500/1‑36 (MQ=255)
cagcaCAGCGAAAAGTTGTATAGCCTGGCAAAagag < 1:4794460/36‑1 (MQ=255)
cagcaCAGCGAAAAGTTGTATAGCCTGGCAAAagag > 1:5579036/1‑36 (MQ=255)
agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGc > 1:963518/1‑36 (MQ=255)
agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGc > 1:1631632/1‑36 (MQ=255)
agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGc > 1:5232711/1‑36 (MQ=255)
gcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGcc > 1:4400398/1‑36 (MQ=255)
cacaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCa < 1:5678951/36‑1 (MQ=255)
acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa > 1:4926619/1‑36 (MQ=255)
acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa > 1:5383683/1‑36 (MQ=255)
acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa > 1:5665726/1‑36 (MQ=255)
aGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaaaa < 1:2783004/36‑1 (MQ=255)
aGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaaaa > 1:4308913/1‑36 (MQ=255)
aGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaaaa > 1:573803/1‑36 (MQ=255)
ggaaaaaGTTGAATAGCCGGGAAAAAGAGCCAAAAc < 1:1619556/33‑1 (MQ=16)
gaaaaaaTTTGTATAGCCTGGAAAAAGAGCCAAAAc < 1:2360972/33‑1 (MQ=25)
gCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAAc > 1:435355/1‑36 (MQ=255)
cGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACg > 1:2875280/1‑36 (MQ=255)
gAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACgg > 1:1072666/1‑36 (MQ=255)
aaaGTTTTATAGCCTGGAAAAAGAGCCAAAACGTtt < 1:3796968/36‑1 (MQ=25)
aaaGTTGTATAGCCTGGCAAAAGAGCCAAAACGGtt < 1:5876419/36‑1 (MQ=255)
aaGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTg < 1:4798199/36‑1 (MQ=255)
aGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGa > 1:2983583/1‑36 (MQ=255)
gTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGAt > 1:4759317/1‑36 (MQ=255)
gTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGAt > 1:1677778/1‑36 (MQ=255)
gTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGAt > 1:5205346/1‑36 (MQ=255)
ttGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATc > 1:3065649/1‑36 (MQ=255)
tGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATcc > 1:3883175/1‑36 (MQ=255)
tGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATcc > 1:2532401/1‑36 (MQ=255)
gTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCt < 1:118421/36‑1 (MQ=255)
gTATACCCTGGAAAAAGAGCCAAAACGGTTGATCCt < 1:4297992/36‑1 (MQ=37)
tataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTa < 1:3397795/36‑1 (MQ=255)
tataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTa < 1:5937429/36‑1 (MQ=255)
ataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTaa < 1:3963234/36‑1 (MQ=255)
taGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAAt > 1:537524/1‑36 (MQ=255)
taGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAAt > 1:3396892/1‑36 (MQ=255)
aGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATc > 1:4118772/1‑36 (MQ=255)
aGCATGGCAAAAGAGCCAAAACGGTTGATCCTAATc < 1:5698641/36‑1 (MQ=255)
gacTGGAAAAAGAGCCATAACGGTTGATCCTAATcc < 1:513104/34‑1 (MQ=37)
gCCTGGCAAAAGAGCCAAAACGGTTGATCCTCATcc > 1:251650/1‑36 (MQ=255)
gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc > 1:4869809/1‑36 (MQ=255)
gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc > 1:3731765/1‑36 (MQ=255)
gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc > 1:603769/1‑36 (MQ=255)
gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc > 1:1684602/1‑36 (MQ=255)
gCCTGGCAAAAGAGCCAAAACGGTTGATCATAATcc > 1:1499712/1‑36 (MQ=255)
cTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCg < 1:2882061/36‑1 (MQ=255)
cTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCg > 1:828606/1‑36 (MQ=255)
ggCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGa > 1:5451149/1‑36 (MQ=255)
gCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGAt > 1:2597708/1‑36 (MQ=255)
gCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGAt < 1:3989561/36‑1 (MQ=255)
gCAAAAGAGCCAAAACGGTTGATCCTAATCCAGGAt > 1:1748543/1‑36 (MQ=255)
cAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATg > 1:1917220/1‑36 (MQ=255)
cAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATg < 1:599095/36‑1 (MQ=255)
cAAAAGAGCCAAAACGGTTGATCCTAATCCAGGATg > 1:530656/1‑36 (MQ=255)
|
CCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATG > REL606/2585778‑2585848
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 20 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|