Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 2,131,584 | A→C | 4.9% | intergenic (‑243/+57) | gatY ← / ← fbaB | D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,131,584 | 0 | A→C | 4.9% | 2.4 ‑2.3 ‑0.8 | 62 | intergenic (‑243/+57) | gatY/fbaB | D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase |
Reads supporting (aligned to +/- strand): new base (3/0): ref base (26/32): total (29/33) | ||||||||||
Fisher's exact test strand distribution p-value = 1.02e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.76e-01 |
GCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCC > REL606/2131550‑2131619
|
gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg < 1:6671208/36‑1 (MQ=255)
gCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg < 1:3764386/36‑1 (MQ=255)
cATGATAAGATGCNTCAGCGTCGCATCAGGCATAgg > 1:4428406/1‑36 (MQ=255)
tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt < 1:4369716/36‑1 (MQ=255)
tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt < 1:6846999/36‑1 (MQ=255)
tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt < 1:1843797/36‑1 (MQ=255)
tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt < 1:6604558/36‑1 (MQ=255)
tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt < 1:6123701/36‑1 (MQ=255)
aTAAGACGCTTCTGCGTCGCATCAGGCATAGGTTGc < 1:1498379/36‑1 (MQ=255)
aTAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGc > 1:5811566/1‑36 (MQ=255)
aaGACGCTTCAGCGTCGCATCAGGCATAGGTTGACg > 1:2283650/1‑36 (MQ=255)
aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg < 1:4356025/36‑1 (MQ=255)
aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg > 1:4516250/1‑36 (MQ=255)
aGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg < 1:6670155/36‑1 (MQ=255)
gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGa < 1:2109609/36‑1 (MQ=255)
tcGCGTCAGCGACGCATCAGGCATCGGCTGGGGGAt > 1:1664109/2‑36 (MQ=17)
aCGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt > 1:2984003/1‑36 (MQ=255)
aCGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt < 1:3955704/36‑1 (MQ=255)
cGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATg < 1:6665006/36‑1 (MQ=255)
cGCGTCAGCGACGCATCAGGCATCGGCTGGGGGAta > 1:6432677/1‑35 (MQ=17)
cGCGTCAGCGACGCATCAGGCATCGGCTGGGGGAta > 1:4501669/1‑35 (MQ=17)
cctCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg < 1:5001276/34‑1 (MQ=255)
cTTCAGCTTCGCATCAGGCATAGGTTGCCGGATgcg < 1:4814064/36‑1 (MQ=255)
cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg < 1:5941730/36‑1 (MQ=255)
cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATgcg > 1:6896467/1‑36 (MQ=255)
ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg > 1:5149809/1‑36 (MQ=255)
ttCAGCGTCGCATCAGGCATAGGTTGCCGGATgcgg > 1:797325/1‑36 (MQ=255)
ttCAGCGTCGCATCAGGCAGAGGTTGCCGGATgcgg > 1:384614/1‑36 (MQ=255)
tgcAGCGTCGCATCAGGCATGGGGTGCCGGATgcgg < 1:3180988/34‑1 (MQ=11)
tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc < 1:2509702/36‑1 (MQ=255)
tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc > 1:59192/1‑36 (MQ=255)
tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc < 1:5766060/36‑1 (MQ=255)
tCAGCGTCGCATCAGGCATAGGTTGCCGGATgcggc > 1:1734859/1‑36 (MQ=255)
tCAGCGTCGCATCAGGCATAGGTTGCCGGATgaggc > 1:5710470/1‑36 (MQ=255)
cAGCGTCGCATCAGGCATAGGTTGCCGGATgcggcg > 1:6654778/1‑36 (MQ=255)
aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt > 1:3998741/1‑36 (MQ=40)
aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt < 1:4773041/36‑1 (MQ=40)
aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt < 1:1871282/36‑1 (MQ=40)
aGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGt < 1:1609486/36‑1 (MQ=40)
gCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTa < 1:4648616/36‑1 (MQ=255)
gCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTa < 1:4823808/36‑1 (MQ=255)
gCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTa < 1:5713506/36‑1 (MQ=255)
gTCGCATCAGGCATAGGTTGCCGGATGCGGCGTaaa < 1:54143/36‑1 (MQ=255)
gTCGCATCAGGCATAGGTTGCCGGATGCGGCGTaaa > 1:4946421/1‑36 (MQ=255)
cGCATCAGGCATAGGTTGCCGGATGCGGCGTAAACg > 1:6434255/1‑36 (MQ=255)
cGCATCAGGCATACGTTGCCGGATGCGGCGTAAACg > 1:3250586/1‑36 (MQ=255)
cATCAGGCATAGGTTGCCGGATGCGGCGTAAACGcc < 1:712012/36‑1 (MQ=255)
cATCAGGCATAGGTTGCCGGATGCGGCGTAAACGcc > 1:4191562/1‑36 (MQ=255)
aTCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCt < 1:2811470/36‑1 (MQ=255)
aTCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCt < 1:6633207/36‑1 (MQ=255)
tCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCtt < 1:4205536/36‑1 (MQ=255)
cAGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTa > 1:5707181/1‑36 (MQ=255)
aGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTAt > 1:4776511/1‑36 (MQ=255)
aGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTAt > 1:4603756/1‑36 (MQ=255)
gCATAGGTTGCCGGATGCGGCGTAAACGCCTTATcc < 1:3357031/36‑1 (MQ=40)
aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt < 1:2048824/36‑1 (MQ=38)
aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt > 1:5736483/1‑36 (MQ=38)
aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt > 1:3424653/1‑36 (MQ=38)
aTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGt > 1:3620755/1‑36 (MQ=38)
aTAGGTTGCCGGATGCGGCGTAAACCCCTTATCCGt > 1:1481627/1‑36 (MQ=255)
aTAGGTTGCCGGAGGCGGTGTAAACGCCTTTTCCga > 1:2848833/1‑35 (MQ=11)
tAGATTGCCGGATGCGGCGTAAACGCCTTATCCGTc < 1:32392/36‑1 (MQ=12)
aGGTTGCCGGATGCGGCGTAAACGCCTTATCCGTcc > 1:1619652/1‑36 (MQ=35)
|
GCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCC > REL606/2131550‑2131619
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 15 ≤
ATCG
< 24 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|