Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,730,188 | 0 | G→A | 10.2% | 3.1 ‑10.6 0.0 | 97 | 2098 | hypF | carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases |
Rejected: Prediction has biased read strand distribution supporting polymorphism. | ||||||||||
Reads supporting (aligned to +/- strand): new base (1/11): ref base (62/23): total (63/34) | ||||||||||
Fisher's exact test strand distribution p-value = 2.61e-05 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
ATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATAC > REL606/2730153‑2730223
|
aTAATGCGCCAGACGTGCACGCAGCAAACGGTTATg > 1:2694348/1‑36 (MQ=255)
tAATGCGCCAGACGTGCACGCAGCAAACGGTTATGa < 1:7205247/36‑1 (MQ=255)
tAATGCGCCAGACGTGCACGCAGCAAACGGTTATGa > 1:3243117/1‑36 (MQ=255)
tAATGCGCCAGACGTGCACGCAGCAAACGGTTATGa > 1:645827/1‑36 (MQ=255)
aaTGCGCCAGACGTGCACGCAGCAAACGGTTATGaa < 1:7624395/36‑1 (MQ=255)
aaTGCGCCAGACGTGCACGCAGCAAACGGTTATGaa < 1:5232107/36‑1 (MQ=255)
aaTGCGCCAGACGTGCACGAAGCAAACGGTTATGaa < 1:5565360/36‑1 (MQ=255)
aTGCGCCAGACGTGCACGGAGCAAACGGTTATGAAt > 1:5440429/1‑36 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt < 1:3083495/36‑1 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt > 1:7986335/1‑36 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt > 1:6046320/1‑36 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt > 1:6148075/1‑36 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt > 1:2356353/1‑36 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt > 1:2264760/1‑36 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt < 1:4103923/36‑1 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt < 1:4406359/36‑1 (MQ=255)
aTGCGCCAGACGTGCACGCAGCAAACGGGTATGAAt > 1:6396434/1‑36 (MQ=255)
tGCGCCAGACGTGCACGCAGCAAACGGTTATGAATa > 1:3168689/1‑36 (MQ=255)
gcgcCAGACGTGCACGCAGCAAACGGTTATGAATaa < 1:5454741/36‑1 (MQ=255)
gcgcCAGACGTGCACGCAGCAAACGGTTATGAATaa > 1:33621/1‑36 (MQ=255)
gcCAGACGTGCACGCAGCAAACGGTTATGAATAAcc > 1:5825479/1‑36 (MQ=255)
ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc < 1:3195962/36‑1 (MQ=255)
ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc < 1:5810501/36‑1 (MQ=255)
ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc > 1:3446615/1‑36 (MQ=255)
ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc > 1:5836494/1‑36 (MQ=255)
cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc > 1:3091385/1‑36 (MQ=255)
cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc < 1:2396399/36‑1 (MQ=255)
cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc < 1:2368406/36‑1 (MQ=255)
cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc < 1:6933777/36‑1 (MQ=255)
cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc < 1:1203135/36‑1 (MQ=255)
cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc > 1:5168893/1‑36 (MQ=255)
cAGACGTGCACGAAGCAAACGGTTATGAATAAcccc < 1:5833441/36‑1 (MQ=255)
aGACGTGCACGCAGCAAACGGTTATGAATAACCccg < 1:7129924/36‑1 (MQ=255)
aGACGTGCACGCAGCAAACGGTTATGAATAACCccg > 1:6986060/1‑36 (MQ=255)
aGACGTGCACGCAGCAAACGGTTATGAATAACCccg < 1:6003955/36‑1 (MQ=255)
aGACGTGCACGCAGCAAACGGTTATGAATAACCccg > 1:8080330/1‑36 (MQ=255)
aGACGTGCACGCAGCAAACGGTTATGAATAACCccg < 1:2245692/36‑1 (MQ=255)
aGACGTGCACGCAGCAAACGGTTATGAATAACCccg < 1:6620244/36‑1 (MQ=255)
gACGTGCACGCAGCAAACGGTTATGAATAACCccgc > 1:2256144/1‑36 (MQ=255)
gACGTGCACGCAGCAAACGGTTATGAATAACCccgc < 1:2132675/36‑1 (MQ=255)
aCGTGCACGCAGCAAACGGTTATGAATAACCccgcc < 1:7280489/36‑1 (MQ=255)
aCGTGCACGCAGCAAACGGGTATGAATAACCccgcc > 1:6391047/1‑36 (MQ=255)
aCGTGCACGAAGCAAACGGTTATGAATAACCccgcc < 1:7602496/36‑1 (MQ=255)
cGTGCACGCAGCAAACGGTTATGAATAACCccgccg > 1:2761372/1‑36 (MQ=255)
cGTGCACGCAGAAAACGGTTATGAATAACCccgccg < 1:1401010/36‑1 (MQ=255)
gTGCACGCAGCAAACGGTTATGAATAACCCcgccgc > 1:8106939/1‑36 (MQ=255)
tGCACGCAGCAAACGGTTATGAATAACCCCGCCGCt > 1:4805776/1‑36 (MQ=255)
gCACGCAGCAAACGGTTATGAATAACCCCGCCGCTa > 1:7262257/1‑36 (MQ=255)
gCACGCAGCAAACGGTTATAAATAACCCCGCCGCTa > 1:8184735/1‑36 (MQ=255)
cACGCAGCAAACGGTTATGAATAACCCCGCCGCTaa > 1:6544376/1‑36 (MQ=255)
cACGCAGCAAACGGTTATGAATAACCCCGCCGCTaa > 1:7300164/1‑36 (MQ=255)
aCGCAGCAAACGGTTATGAATAACCCCGCCGCTaaa > 1:7161363/1‑36 (MQ=255)
gcagcaAACGGTTATGAATAACCCCGCCGCTaaaca > 1:3877326/1‑34 (MQ=255)
gcagcaAACGGTTATGAATAACCCCGCCGCTAAATa > 1:8363715/1‑36 (MQ=255)
cagcaAACGGTTATGAATAACCCCGCCGCTAAATAc > 1:1813910/1‑36 (MQ=255)
cagcaAACGGTTATGAATAACCCCGCCGCTAAATAc > 1:2662436/1‑36 (MQ=255)
cagcaAACGGTTATGAATAACCCCGCCGCTAAATAc > 1:4114079/1‑36 (MQ=255)
agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc > 1:8153725/1‑36 (MQ=255)
agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc > 1:2870555/1‑36 (MQ=255)
agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc > 1:3900174/1‑36 (MQ=255)
agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc > 1:4067528/1‑36 (MQ=255)
gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa > 1:4485348/1‑36 (MQ=255)
gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa > 1:2406566/1‑36 (MQ=255)
gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa > 1:6155646/1‑36 (MQ=255)
gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa > 1:5430528/1‑36 (MQ=255)
gaaaaCGGTTATAAAAAACCCCGCCGCTAAATACCa < 1:1720970/34‑1 (MQ=25)
aaaaaGGGTTATAAAAAACCCCGCCGCTAAATACCa < 1:1878477/34‑1 (MQ=21)
caAACGGTTATGAATAACCCCGCCGCTAAATACCAg > 1:2931677/1‑36 (MQ=255)
caAACGGTTATGAATAACCCCGCCGCTAAATACCAg > 1:5909481/1‑36 (MQ=255)
caAACGGTTATGAATAACCCCGCCGCTAAATACCAg > 1:5107312/1‑36 (MQ=255)
aaaaCGTTTATAAATAACCCCGCCGCTAAATACCAg < 1:2150293/35‑1 (MQ=37)
aaaaCGTTTATAAAAAACCCCGCCGCTAAATACCAg < 1:5840656/35‑1 (MQ=21)
aaaCGGTTATGAATAACCCCGCCGCTAAATACCAGc > 1:2781100/1‑36 (MQ=255)
aaaCGGTTATGAATAACCCCGCCGCTAAATACCAGc > 1:4831770/1‑36 (MQ=255)
aaaCGGTTATAAATAACCCCGCCGCTAAATACCAGc < 1:1190200/36‑1 (MQ=255)
aaCGGTTATGAATAACCCCGCCGCTAAATACCAGCg > 1:2517254/1‑36 (MQ=255)
aaCGGTTATGAATAACCCCGCCGCTAAATACCAGCg > 1:7903269/1‑36 (MQ=255)
aCGGTTATGAATAACCCCGCCGCTAAATACCAGCGt > 1:1743421/1‑36 (MQ=255)
aCGGTTATGAATAACCCCGCCGCTAAATACCAGCGt > 1:7815693/1‑36 (MQ=255)
cgtttATAAAAAACCCCGCCGCTAAATACCAGCGTa < 1:7807256/33‑1 (MQ=21)
cGGTTATAAAAAACCCCGCCGCTAAATACCAGCGTa < 1:2236492/36‑1 (MQ=37)
ggTTATGAATAACCCCGCCGCTAAATACCAGCGTAg > 1:8327851/1‑36 (MQ=255)
ggTTATGAATAACCCCGCCGCTAAATACCAGCGTAg > 1:6508528/1‑36 (MQ=255)
ggTTATGAATAACCCCGCCGCTAAATACCAGCGTAg > 1:7860099/1‑36 (MQ=255)
ggTTATGAATAACCCCGCCGCTAAATACCAGCATAg > 1:128444/1‑36 (MQ=255)
ggTTAAAAATAACCCCGCCGCTAAATACCAGCGTAg < 1:3955380/36‑1 (MQ=37)
gTTATGAATAACCCCGCCGCTAAATACCAGCGTAGt > 1:5888046/1‑36 (MQ=255)
ttATAAATAACCCCGCCGCTAAATACCAGCGTAGTg < 1:4985828/36‑1 (MQ=255)
ttATAAAAAACCCCGCCGCTAAATACCAGCGTAGTg < 1:2752921/36‑1 (MQ=37)
tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa > 1:7896713/1‑36 (MQ=255)
tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa > 1:5480410/1‑36 (MQ=255)
tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa > 1:6509416/1‑36 (MQ=255)
tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa > 1:6531673/1‑36 (MQ=255)
tATAAATAACCCCGCCGCTAAATACCAGCGTAGTGa < 1:5239886/36‑1 (MQ=255)
aTGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt > 1:7500001/1‑36 (MQ=255)
aTGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt > 1:5366843/1‑36 (MQ=255)
aTGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt > 1:6472944/1‑36 (MQ=255)
tGAATAACCCCGCCGCTAAATACCAGCGTAGTGATa > 1:5666058/1‑36 (MQ=255)
gAATAACCCCGCCGCTAAATACCAGCGTAGTGATAc < 1:3420138/36‑1 (MQ=255)
gAATAACCCCGCCGCTAAATACCAGCGTAGTGATAc > 1:3881958/1‑36 (MQ=255)
|
ATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATAC > REL606/2730153‑2730223
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 13 ≤
ATCG
< 24 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|