Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 4,266,649 | C→A | 9.7% | intergenic (‑345/‑48) | acs ← / → nrfA | acetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 4,266,649 | 0 | C→A | 9.7% | 5.1 ‑2.8 0.0 | 53 | intergenic (‑345/‑48) | acs/nrfA | acetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome |
Reads supporting (aligned to +/- strand): new base (2/4): ref base (35/12): total (37/16) | ||||||||||
Fisher's exact test strand distribution p-value = 6.02e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
ATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTC > REL606/4266615‑4266683
|
aTGACAGTGTATGTGCGGTTACTGTCGTGAAAAACa < 1:782258/36‑1 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa > 1:5795566/1‑36 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa > 1:2389429/1‑36 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa > 1:7266178/1‑36 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa > 1:7194884/1‑36 (MQ=255)
tGACAGTGTAGGTGCGGTTACTTTCGTGAAAAACaa > 1:2129985/1‑36 (MQ=255)
tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACaa < 1:400565/36‑1 (MQ=255)
tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACaa > 1:6163886/1‑36 (MQ=255)
tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACaa < 1:6998740/36‑1 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt > 1:5463755/1‑36 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt > 1:7596362/1‑36 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt > 1:3732881/1‑36 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt < 1:3959198/36‑1 (MQ=255)
cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa > 1:2285525/1‑36 (MQ=255)
cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa < 1:4828138/36‑1 (MQ=255)
cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa > 1:7925611/1‑36 (MQ=255)
cAGTGTAAGTGCGGTTACTGTCGTGAAAAACAATaa < 1:5407866/36‑1 (MQ=255)
agggtAGGTGCGGTTACTGTCGTGTAAAAAAATaaa > 1:3269903/4‑36 (MQ=25)
aGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAa > 1:1397330‑M2/1‑35 (MQ=255)
gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa > 1:2262089‑M2/1‑34 (MQ=255)
gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa < 1:7278335‑M2/36‑3 (MQ=255)
tgtAGGTGCGGTTACTGTCGTGAAAATCAATataaa > 1:3913464/1‑36 (MQ=37)
tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa < 1:2036716‑M2/36‑4 (MQ=255)
tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa < 1:2729121‑M2/36‑4 (MQ=255)
tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa > 1:3354270‑M2/1‑33 (MQ=255)
gtAGGTGCGGTTACTGTTGTGAAAAACAATAAaaac > 1:6458527‑M2/1‑32 (MQ=255)
gtAGGCGCGCTTACTGTAGTGAAAAACAATAAAAAc < 1:4256303/36‑1 (MQ=25)
tAGGTGCGGTTACTGTCGTGAAAAACAATAAaaacc > 1:5475810‑M2/1‑31 (MQ=255)
tAGGTGCGGTTACTGTCGTGAAAAACAATAAaaacc < 1:1786231‑M2/36‑6 (MQ=255)
gggCGGTTACTGCAGTGAAAAACAATAAAAACCGcc < 1:6051002/34‑1 (MQ=37)
tGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCa > 1:4912596/1‑36 (MQ=255)
cGGTTACTGTCGTGAAAAACAATAAAAACCGCCAtt > 1:868948/1‑36 (MQ=255)
cGGTTACTGTCGTGAAAAAAAATAAAAACCGCCAtt > 1:61833/1‑36 (MQ=255)
ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg > 1:1215817/1‑36 (MQ=255)
ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg > 1:6793296/1‑36 (MQ=255)
gTTACTGTCGTAAAAAACAATAAAAACCGCCATTGc > 1:3563058/1‑36 (MQ=255)
ttACTGTCGTGAAAAACAATAAAAACCGCCATTGca > 1:1563161/1‑36 (MQ=255)
tCGTGAAAAACAAGAAAAACCGCCATTGCAACAATg > 1:2285148/1‑36 (MQ=255)
gTGAAAAACAATAAAAACCGCCATTGCAACAACGgc > 1:4217963/1‑36 (MQ=255)
aaaaaaaaAAAAAAAACCCGCCATTGCAACAATGgc < 1:5309197/33‑1 (MQ=16)
tGAAAAACAATAAAAACCGCCATTGCAACAATGgcg > 1:7850141/1‑36 (MQ=255)
aaaaaaaAAAAAAAACCCGCCATTGCAACAATGgcg < 1:122419/34‑1 (MQ=21)
gAAAAACAATAAAAACCGCCATTGCAACAATGgcgc > 1:7584869/1‑36 (MQ=255)
aaaaaaAAAAAAAAACCGCCATTGCAACAATGgcgc < 1:3108806/35‑1 (MQ=25)
aaaaaCAATAAAAACCGCCATTGCAACAATGGCGCa > 1:2059766/1‑36 (MQ=255)
aaaaCAATAAAAACCGCCATTGCAACAATGGCGCaa > 1:1845099/1‑36 (MQ=255)
aaaaCAAAAAaaaccgccgcagattggtatggcaaa > 1:4908588‑M2/1‑10 (MQ=255)
aaaCAATAAAAACCGCCCTTGCAACAATGCCGCaac > 1:6462075/1‑35 (MQ=25)
aaaAAAAAAAAACCGCCATTGCAACAATGGCGCAAt < 1:3095111/36‑1 (MQ=37)
aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt > 1:1741131/1‑36 (MQ=255)
aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt > 1:1596554/1‑36 (MQ=255)
aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt > 1:5047552/1‑36 (MQ=255)
aCAATAAAAACCGCCATTGCAACAATGGCGCAATTc > 1:819376/1‑36 (MQ=255)
|
ATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTC > REL606/4266615‑4266683
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 9 ≤
ATCG
< 22 ≤
ATCG
< 31 ≤
ATCG
< 37 ≤
ATCG
|