Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,136,629 | 0 | G→C | 11.8% | 3.1 ‑6.8 0.0 | 74 | 319 | uppP | undecaprenyl pyrophosphate phosphatase |
Rejected: Prediction has biased read strand distribution supporting polymorphism. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/10): ref base (35/29): total (35/39) | ||||||||||
Fisher's exact test strand distribution p-value = 1.14e-03 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
GTTAAACAATGACTTAATCGTGTCGTGGAACAACAGCCCCAGTACCACCGCCGGAATCATCCCCAGCAAA > REL606/3136594‑3136663
|
gTTAAACAATGACTTAATCGTGTCGTGGAACAACAg > 1:4098714/1‑36 (MQ=255)
gTTAAACAATGACTTAATCGTGTCGTGGAACAACAg > 1:5989578/1‑36 (MQ=255)
gTTAAACAATGACTTAATCGTGTCGTGGAACAACAg > 1:6638980/1‑36 (MQ=255)
ttAAACAATGACTTAATCGTGTCGTGGAACAACAGc > 1:5852923/1‑36 (MQ=255)
ttAAACAATGACTTAATCGTGTCGTGGAACAACAGc < 1:3923779/36‑1 (MQ=255)
ttAAACAATGACTTAATCGTGTCGTGGAACAACAGc < 1:7809208/36‑1 (MQ=255)
tAAACAATGACTTAATCGTGTCGTGGAACAACAGcc > 1:4870969/1‑36 (MQ=255)
tAAACAATGACTTAATCGTGTCGTGGAACAACAGcc < 1:3760578/36‑1 (MQ=255)
aaaCAATGACTTAATCGTGTCGTGGAACAACAGccc < 1:2956495/36‑1 (MQ=255)
aaaCAATGACTTAATCGTGTCGTGGAACAACAGccc < 1:7738110/36‑1 (MQ=255)
aaCAATGACTTAATCGTGTCGTGGAACAACAGcccc > 1:5230280/1‑36 (MQ=255)
aaCAATGACTTAATCGTGTCGTGGAACAACAGcccc < 1:7050217/36‑1 (MQ=255)
aCAATGACTTAATCGTGTCGTGGAACAACAGCCCCa < 1:4133417/36‑1 (MQ=255)
aCAATGACTTAATCGTGTCGTGGAACAACAGCCCCa > 1:3252734/1‑36 (MQ=255)
aCAATGACGTAATCGTGTCGTGGAACAACAGCCCCa > 1:7145747/1‑36 (MQ=255)
cAATGACTTAATCGTGTCGTGGAAGAACAGCCCCAg < 1:1500824/36‑1 (MQ=255)
cAATGACTTAATCGTGTCGTGGAACAACAGCCCCAg > 1:5374720/1‑36 (MQ=255)
cAATGACTTAATCGTGTCGTGGAACAACAGCCCCAg > 1:868522/1‑36 (MQ=255)
aaTGACTTAATCGTGTCGTGGAACAACAGCCCCAGt > 1:1288461/1‑36 (MQ=255)
aaTGACTTAATCGTGTCGTGGAACAACAGCCCCAGt < 1:5803891/36‑1 (MQ=255)
aaTGACTTAATCGTGTCGTGGAACAACAGCCCCAGt < 1:1458474/36‑1 (MQ=255)
aTGACTTAATCGTGTCGTGGAACAACAGCCCCAGTa < 1:1085255/36‑1 (MQ=255)
aTGACTTAATCGTGTCGTGGAACAACAGCCCCAGTa > 1:4323293/1‑36 (MQ=255)
tggcTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:5498252/33‑1 (MQ=255)
tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:3672270/36‑1 (MQ=255)
tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:5265487/36‑1 (MQ=255)
tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:1526084/36‑1 (MQ=255)
tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:5347728/36‑1 (MQ=255)
tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:2166016/36‑1 (MQ=255)
tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac < 1:6892232/36‑1 (MQ=255)
aCTTAATCGTGTCGTGGAACAACAGCCCCAGTacca > 1:7406899/1‑36 (MQ=255)
aCTTAATCGTGTCGTGGAACAACAGCCCCAGTacca > 1:1509211/1‑36 (MQ=255)
ttAATCGTGTCGTGGAACAACAGCCCCAGTaccacc < 1:7805924/36‑1 (MQ=255)
ttAATCGTGTCGTGGAACAACAGCCCCAGTaccacc < 1:4679966/36‑1 (MQ=255)
taaaTCGTGTCGTGGAACAACAGCCCCAGTaccacc < 1:4602271/34‑1 (MQ=255)
tAATCGTGTCGTGTAACAACAGCCCCAGTCCCAccg < 1:5215277/36‑1 (MQ=37)
tAATCGTGTCGTGGAACAACAGCCCCAGTACCAccg < 1:5286297/36‑1 (MQ=255)
aaTCGTGTCGTGTAACAACAGCCCCAGTACCAccgc < 1:5919715/36‑1 (MQ=255)
aaTCGTGTCGTGGAACAACAGCCCCAGTACCAccgc > 1:6303706/1‑36 (MQ=255)
aaTCGTGTCGTGGAACAACAGCCCCAGTACCAccgc > 1:3984553/1‑36 (MQ=255)
aTCGTGTCGTGGAACAACAGCCCCAGTACCAccgcc > 1:4503467/1‑36 (MQ=255)
aTCGTGTCGTGGAACAACAGCCCCAGTACCAccgcc > 1:1115682/1‑36 (MQ=255)
aTCGTGTCGTGGAAAAACAGCCCCAGTACCAccgcc < 1:1688634/36‑1 (MQ=255)
aTCGTGTCGTGAAAAAACAGCCCCAGTACCAccgcc < 1:991633/36‑1 (MQ=37)
tcgtgtcgtgGAACAACAGCCCCAGTACCAccgccg > 1:5116978/1‑36 (MQ=255)
tcgtgtcgtgGAACAACAGCCCCAGTACCAccgccg > 1:1797960/1‑36 (MQ=255)
cgtgtcgtgGAACAACAGCCCCAGTACCACCGCCgg > 1:2828808/1‑36 (MQ=255)
gtgtcgtgGAACAACAGCCCCAGTACCACCGCCGGa > 1:7220419/1‑36 (MQ=255)
gtgtcgtgGAACAACAGCCCCAGTACCACCGCCGGa > 1:6809480/1‑36 (MQ=255)
tgtcgtgGAACAACAGCCCCAGTACCACCGCCGGaa > 1:6988092/1‑36 (MQ=255)
tgtcgtgGAACAACAGCCCCAGTACCACCGCCGGaa > 1:3453596/1‑36 (MQ=255)
gtcgtgGAACAACAGCCCCAGTACCACCGCCGGAat > 1:5800058/1‑36 (MQ=255)
gtcgggGAACAACAGCCCCAGTACCACCGCCGGAat < 1:678882/36‑1 (MQ=255)
cgtgGAACAACAGCCCCAGTACCACCGCCGGAatca > 1:2992174/1‑36 (MQ=255)
cgtgGAACAACAGCCCCAGTACCACCGCCGGAatca > 1:6194926/1‑36 (MQ=255)
cgtgGAACAACAGCCCCAGTACCACCGCCGGAatca > 1:4131103/1‑36 (MQ=255)
cgtgAAACAACACCCCCAGTACCACCGCCGGAatca < 1:2915895/36‑1 (MQ=38)
gtgGAACAACAGCCCCAGTACCACCGCCGGAatcat > 1:2723366/1‑36 (MQ=255)
ggAACAACAGCCCCAGTACCACCGCCGGAATCATcc > 1:4098557/1‑36 (MQ=255)
ggAACAACAGCCCCAGTACCACCGCCGGAATCATcc > 1:6494079/1‑36 (MQ=255)
aaacaacAGCCCCAGTACCACCGCCGGAATCATccc < 1:4699676/35‑1 (MQ=255)
aacaccACCCCAACTCCCACCGCCGGAATCATcccc < 1:4768733/36‑1 (MQ=16)
aacaacAGCCCCAGTACCACCGCCGGAATCATcccc > 1:2232078/1‑36 (MQ=255)
aacaacACCCCCAGTACCACCGCCGGAATCATcccc < 1:2314648/36‑1 (MQ=255)
aacaacACCCCCAGTACCACCGCCGGAATCATcccc < 1:4757652/36‑1 (MQ=255)
aacaacACCCCCAAGGCCACCGCCGGAATCATcccc < 1:541416/36‑1 (MQ=19)
aaaaacaCCCCCAGCCCCCCCGCCGGAATCATcccc < 1:4087094/33‑1 (MQ=12)
ccaacaCCCCCATTACCACCGCCGGAATCATCCCCa < 1:1876333/35‑1 (MQ=37)
caaacaGCCCCAGTACCACCGCCGGAATCATCCCCa < 1:274305/34‑1 (MQ=255)
acaacaGCCCCAGTACCACCGCCGGAATCATCCCCa > 1:7038/1‑36 (MQ=255)
acaacaCCCCAAGTACCACCGCCGGAATCATCCCCa < 1:1400367/36‑1 (MQ=38)
aaaacaCCCCCATTACCACCGCCGGAATCATCCCCa < 1:6988538/34‑1 (MQ=37)
aaaacaCCCCCACTGCCCCCGCCGGAATCATCCCCa < 1:1803113/34‑1 (MQ=16)
caacaGCCCCAGTACCACAGCCGGAATCATCCCCAg > 1:725395/1‑36 (MQ=255)
caGCCCCAGTACCACCGCCGGAATCATCCCCAGCaa > 1:730022/1‑36 (MQ=255)
caGCCCCAGGACCACCGCCGGAATCATCCCCAGCaa > 1:5453561/1‑36 (MQ=255)
aGCCCCAGTACCACCGCCGGAATCATCCCCAGCaaa > 1:2402792/1‑36 (MQ=255)
|
GTTAAACAATGACTTAATCGTGTCGTGGAACAACAGCCCCAGTACCACCGCCGGAATCATCCCCAGCAAA > REL606/3136594‑3136663
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 7 ≤
ATCG
< 20 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|