Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,893,550 | 1 | .→G | 2.9% | 4.0 0.0 ‑0.5 | 70 | intergenic (+5/‑51) | kup/insJ‑5 | potassium transporter/IS150 hypothetical protein |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Reads supporting (aligned to +/- strand): new base (1/1): ref base (33/35): total (34/36) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e+00 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.97e-01 |
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT > REL606/3893515‑3893585
|
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa < 1:6025886/36‑1 (MQ=255)
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa < 1:44577/36‑1 (MQ=255)
tATCGAACTGGGCACTCAGGTCGAAATCTAACGCCa > 1:4258227/1‑36 (MQ=255)
aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g > 1:903041/1‑36 (MQ=255)
aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g < 1:1813694/36‑1 (MQ=255)
cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GAc < 1:4619890/36‑1 (MQ=255)
cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GAc > 1:4329564/1‑36 (MQ=255)
gAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACt > 1:3479372/1‑36 (MQ=255)
gAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACt > 1:6693877/1‑36 (MQ=255)
aCTGGGTACTCAGGTCTAAATCTAACGCCA..GACTGt < 1:2643515/36‑1 (MQ=255)
aCTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGt > 1:5503759/1‑36 (MQ=255)
aCTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGt < 1:2202490/36‑1 (MQ=255)
cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa < 1:6799242/36‑1 (MQ=255)
cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa < 1:3520982/36‑1 (MQ=255)
cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa > 1:6115126/1‑36 (MQ=255)
cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa > 1:1339050/1‑36 (MQ=255)
tGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTAc < 1:3359013/36‑1 (MQ=255)
gTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGc > 1:2286584/1‑36 (MQ=255)
aCTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAc < 1:5353453/36‑1 (MQ=255)
aCTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAc < 1:1854745/36‑1 (MQ=255)
cTCAGGTCGAAATCTAACGCCA..GACTGTACTTCAcc < 1:6620667/36‑1 (MQ=255)
cTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAcc > 1:7043161/1‑36 (MQ=255)
cTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAcc > 1:3279528/1‑36 (MQ=255)
cTCAGGTAGAAATCTAACGCCA..GACTGTACTGCAcc < 1:1493422/36‑1 (MQ=255)
tCAGGTCGAAATCTAACGCCA..GACTGTACTGCAccc < 1:4145591/36‑1 (MQ=255)
tCAGGTCGAAATCTAACGCCA..GACTGTACTGCAccc < 1:2701724/36‑1 (MQ=255)
cAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCa < 1:6565609/36‑1 (MQ=255)
cAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCa < 1:5218593/36‑1 (MQ=255)
aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt < 1:4210210/36‑1 (MQ=255)
aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt < 1:5184194/36‑1 (MQ=255)
aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt > 1:7305295/1‑36 (MQ=255)
gTCGAAATCTAACGCCA..GACTGTACTGCACCCAttt < 1:2571732/36‑1 (MQ=255)
gTCGAAATCTAACGCCA..GACTGTACTGCACCCAttt > 1:2339688/1‑36 (MQ=255)
gTCGAAATCTAACGCCA..GACTGTACTGCACCCAttt < 1:423694/36‑1 (MQ=255)
tCGAAATCTAACGCCA..GACTGTACTGCACCCAtttt > 1:1086055/1‑36 (MQ=255)
tCGAAATCTAACGCCA..GACTGTACTGCACCCAtttt < 1:2429068/36‑1 (MQ=255)
cGAAATCTAACGCCAG.GACTGTACTGCACCCAtttt > 1:6198436/1‑36 (MQ=25)
cGAAATCTAACGCCA..GACTGTACTGCACCCATTttg < 1:3455184/36‑1 (MQ=255)
cGAAATCTAACGCCA..GACTGTACTGCACCCATTttg < 1:4975839/36‑1 (MQ=255)
gAAATCTAACGCCA..GACTGTACTGCACCCATTttgt < 1:1123347/36‑1 (MQ=255)
gAAATCTAACGCCA..GACTGTACTGCACCCATTttgt > 1:6075048/1‑36 (MQ=255)
aaaTCTAACGCCA..GACTGTACTGCACCCATTttgtt > 1:1538032/1‑36 (MQ=255)
aaTCTAACGCCA..GACTGTACTGCACCCATTttgttg < 1:51526/36‑1 (MQ=255)
aaTCTAACGCCA..GACTGCACTGCACCCATTttgttg > 1:7300889/1‑36 (MQ=255)
aTCTAACGCCAGGGACTGTACTGCACCCATTttgtt < 1:6984964/36‑1 (MQ=17)
aTCTAACGCCA..GACTGTACTGCACCCATTTTGTTgg > 1:5676672/1‑36 (MQ=255)
tCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGa > 1:5291489/1‑36 (MQ=255)
tCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGa < 1:3253005/36‑1 (MQ=255)
tCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGa > 1:4419196/1‑36 (MQ=255)
cTAACGCCA..GACTGTACTGCACCCATTTTGTTGGAc < 1:4358178/36‑1 (MQ=255)
aaCGCCA..GACTGTACTGCACCCATTTTTTTGGACGa > 1:2136078/1‑36 (MQ=255)
aaCGCCA..GACTGTACTGCACCCATTTTGTTGGACGa < 1:5912004/36‑1 (MQ=255)
aCGCCA..GACTGTACTGCACCCATTTTGTTGGACGAt > 1:32132/1‑36 (MQ=255)
aCGCCA..GACTGTACTGCACCCATTTTGTTGGACGAt > 1:4678343/1‑36 (MQ=255)
aCGCCA..GACTGTACTGCACCCATTTTGTTGGACGAt < 1:2988156/36‑1 (MQ=255)
cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg > 1:5065330/1‑36 (MQ=255)
cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg > 1:3168875/1‑36 (MQ=255)
cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg > 1:6500074/1‑36 (MQ=255)
cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg < 1:1085652/36‑1 (MQ=255)
gCCA..GACTGTACTGCACCCATTTTGTTGGACGATGa < 1:1079770/36‑1 (MQ=255)
ccA..GACTGTCCTGCACCCATTGTGTTGGACGATGaa < 1:3930709/36‑1 (MQ=17)
ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:5872723/1‑36 (MQ=255)
ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa < 1:2778370/36‑1 (MQ=255)
ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:452964/1‑36 (MQ=255)
ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:4008109/1‑36 (MQ=255)
ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa < 1:315520/36‑1 (MQ=255)
ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:681215/1‑36 (MQ=255)
cA..GACTGTACTGCCCCCATTTTGTTGGACGATGaaa < 1:5518838/36‑1 (MQ=255)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa < 1:1868054/36‑1 (MQ=40)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:2105405/1‑36 (MQ=40)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:2823119/1‑36 (MQ=40)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa < 1:6644060/36‑1 (MQ=40)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:542311/1‑36 (MQ=40)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:5546571/1‑36 (MQ=40)
cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:5960908/1‑36 (MQ=40)
a..GACTGTACTGCACCCATTTTGTTGGGCGATGAAAt > 1:7389806/1‑36 (MQ=255)
a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt < 1:3422945/36‑1 (MQ=38)
a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt > 1:3878815/1‑36 (MQ=38)
a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt < 1:6047420/36‑1 (MQ=38)
|
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT > REL606/3893515‑3893585
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 16 ≤
ATCG
< 26 ≤
ATCG
< 32 ≤
ATCG
< 37 ≤
ATCG
|