Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 988,180 | 0 | G→. | 4.0% | 0.8 0.0 ‑0.1 | 50 | coding (237/747 nt) | kdsB | 3‑deoxy‑manno‑octulosonate cytidylyltransferase |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (2/0): ref base (34/14): total (36/14) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e+00 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.17e-01 |
ACGCGCGCCGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGA > REL606/988145‑988215
|
aCGCGCGCCGATCATCAGTCAGGAACAGAACGTCTg > 1:1759419/1‑36 (MQ=255)
cgcgcgCCGATCATCAGTCAGGAACAGAACGTCTgg < 1:2682451/36‑1 (MQ=255)
gcgcgcCGATCATCAGTCAGGAACAGAACGTCTGGc > 1:2558471/1‑36 (MQ=255)
cgcgcCGATCATCAGTCAGGAACAGAACGTCTGGCg > 1:1724914/1‑36 (MQ=255)
gcCGATCATCAGTCAGGAACAGAACGTCTGGCGGaa < 1:423380/36‑1 (MQ=255)
cGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGt < 1:4107024/36‑1 (MQ=255)
cGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGt < 1:5499232/36‑1 (MQ=255)
gATCATCAGTCAGGAACAGAACGTCTGGCGGAAGtt > 1:5861426/1‑36 (MQ=255)
gATCATCAGTCAGGAACAGAACGTCTGGCGGAAGtt < 1:1359642/36‑1 (MQ=255)
atcatcAGTCAGGAACAGACCGTCTGGCGGAAGTTg < 1:3840765/36‑1 (MQ=255)
atcatcAGTCAGGAACAGAACGTCTGGCGGAAGTTg > 1:66887/1‑36 (MQ=255)
atcatcAGTCAGGAAAAGAACGTATGGCGGAAGTTg > 1:3569643/1‑36 (MQ=37)
tcatcaGTCAGGAACAGAACGTCTGGCGGAAGTTGt > 1:2601322/1‑36 (MQ=255)
catcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTc < 1:5210107/36‑1 (MQ=255)
catcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTc < 1:1465179/36‑1 (MQ=255)
catcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTc > 1:1056368/1‑36 (MQ=255)
catcaGTCAGAAACAGAACGTCTGGCGGAAGTTGTc < 1:2970938/36‑1 (MQ=255)
atcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTCg > 1:5427402/1‑36 (MQ=255)
atcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTCg > 1:1574694/1‑36 (MQ=255)
taagtcagGAACAGAACGTCTGGCGGGAGTTGTCGa > 1:4383253/3‑36 (MQ=37)
cagtcagGAACAGAACGTCTGGCGGAAGTTGTCGaa < 1:2929734/36‑1 (MQ=255)
agtcagGAACAGAACGTCTGGAGGAAGTTGTCGaaa > 1:2408547/1‑36 (MQ=255)
gtaagGAACAGAACGTCT.GCGGAAGTTGTCGaaaaa > 1:3160544/4‑36 (MQ=37)
cagGAACAGAACGTCTGGCGGAAGTTGTCGAAAAAt > 1:539511/1‑36 (MQ=255)
cagAAACAGAACGTCTGGCGGAAGTTGTCGAAAAAt < 1:2561545/36‑1 (MQ=255)
agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg > 1:4305993/1‑36 (MQ=255)
agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg > 1:3055619/1‑36 (MQ=255)
agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg > 1:1486589/1‑36 (MQ=255)
agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg > 1:1155624/1‑36 (MQ=255)
ggAACAGAACTTCTGGCGGAAGTTGTCGAAAAATgc < 1:5588551/36‑1 (MQ=255)
ggAACAGAACGTCTGGCGGAAGTTGTCGAAAAATgc > 1:634907/1‑36 (MQ=255)
gAACAGAACGTCTGGCGGAAGTTGTCGAAAAATgcg > 1:4633207/1‑36 (MQ=255)
gAACAGAACGTCTGGCGGAAGTTGTCGAAAAATgcg > 1:832661/1‑36 (MQ=255)
gAACAGAACGTCTGGCGGAAGTTGCCGAAAAATggg > 1:1320810/1‑34 (MQ=37)
aCAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCa > 1:1819073/1‑36 (MQ=255)
aGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCCtt > 1:1044367/1‑36 (MQ=255)
aGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCAtt > 1:3173130/1‑36 (MQ=255)
aGAACGTCTGGCGGAAGTTGTCGAAAAATGAGAAtt > 1:5177819/1‑36 (MQ=37)
aGAACGTCTGGAGGAAGTTGTCGAAAAATGCGCAtt > 1:1963048/1‑36 (MQ=255)
aGAACGTCT.GCGGAAGATGTCGACAAATGCGCATTc > 1:292966/1‑36 (MQ=21)
gAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTc > 1:4481000/1‑36 (MQ=255)
gAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTc > 1:4575860/1‑36 (MQ=255)
gAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTc > 1:1939066/1‑36 (MQ=255)
aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCCTTCa > 1:130747/1‑36 (MQ=255)
aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCa > 1:4832614/1‑36 (MQ=255)
aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCa < 1:4797747/36‑1 (MQ=255)
aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCa > 1:4096314/1‑36 (MQ=255)
aCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAg > 1:3983975/1‑36 (MQ=255)
cGTCTGGCGGGAGTTGTCGAAAAATGCGCATTCAGc < 1:679638/36‑1 (MQ=255)
gTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGcg > 1:4576926/1‑36 (MQ=255)
gTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGcg < 1:1005327/36‑1 (MQ=255)
tGGCGGAAGTTGTCGAAAAATGCGCATTCAGcgacg > 1:5524317/1‑36 (MQ=255)
ggCGGAAGTTGTCGAAAAATGCGCATTCAGcgacga > 1:5052528/1‑36 (MQ=255)
|
ACGCGCGCCGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGA > REL606/988145‑988215
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 18 ≤
ATCG
< 27 ≤
ATCG
< 37 ≤
ATCG
|