Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,584,604 | 0 | T→A | 7.0% | 2.1 ‑6.4 ‑0.0 | 70 | 431 | yhiF | predicted DNA‑binding ranscriptional regulator |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/5): ref base (49/16): total (49/21) | ||||||||||
Fisher's exact test strand distribution p-value = 1.68e-03 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.98e-01 |
AACGGAAGATATTCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTACTGCCATAAACACAATATCATGA > REL606/3584569‑3584639
|
aaCGGAAGATATTCTGGAGAAACTGAAAATATCGCt < 1:3472504/36‑1 (MQ=255)
aCGGAAGATATTCTGGAGAAACTGAAAATATCGCTa > 1:2675998/1‑36 (MQ=255)
aCGGAAGATATTCTGGAGAAACTGAAAATATCGCTa > 1:6604344/1‑36 (MQ=255)
aCGGAAGAGATTCTGGAGAAACTGAAAATATCGCTa < 1:6086686/36‑1 (MQ=255)
gtaaGATATTCTGGAGAAACTGAAAATATCGCTaaa < 1:5509537/34‑1 (MQ=255)
ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa > 1:2937580/1‑36 (MQ=255)
ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa > 1:3943199/1‑36 (MQ=255)
ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa > 1:4551059/1‑36 (MQ=255)
ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa > 1:2078957/1‑36 (MQ=255)
ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa < 1:1791646/36‑1 (MQ=255)
gAAGATATTCTGGAGAAACTGAAAATATCGCTaaaa > 1:3886160/1‑36 (MQ=255)
gAAGATATTCTGGAGAAACTGAAAATATCGCTaaaa > 1:6957131/1‑36 (MQ=255)
gAAGATATTCTGGAGAAACTGAAAATATCGCTaaaa < 1:2584657/36‑1 (MQ=255)
aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa > 1:359825/1‑36 (MQ=255)
aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa > 1:147542/1‑36 (MQ=255)
aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa > 1:4943866/1‑36 (MQ=255)
aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa > 1:4287149/1‑36 (MQ=255)
aGATATTCTGGAGAAACTGAAAATATCGCTAAAAAc < 1:981933/36‑1 (MQ=255)
aGATATTCTGGAGAAACTGAAAATATCGCTAAAAAc < 1:4799346/36‑1 (MQ=255)
aGATATTCTGGAGAAACTGAAAATATCGCTAAAAAc > 1:1526290/1‑36 (MQ=255)
gATATTCTGGAGAAACTGAAAATATCGCTAAAAACg > 1:7660349/1‑36 (MQ=255)
gATATTCTGGAGAAACTGAAAATATCGCTAAAAACg > 1:6132742/1‑36 (MQ=255)
atatTCTGGAGAAACTGAAAATATCGCTAAAAACGt > 1:6616577/1‑36 (MQ=255)
atatTCTGGAGAAACTGAAAATATCGCTAAAAACGt < 1:3865640/36‑1 (MQ=255)
atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc < 1:7587091/36‑1 (MQ=255)
atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc < 1:3984101/36‑1 (MQ=255)
atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc < 1:827901/36‑1 (MQ=255)
atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc > 1:2311943/1‑36 (MQ=255)
ttCTGGAGAAACTGAAAATATCGCTAAAAACGTTCt > 1:7533433/1‑36 (MQ=255)
ttCTGGAGAAACTGAAAATATCGCTAAAAACGTTCt < 1:1428112/36‑1 (MQ=255)
tCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTa > 1:1161809/1‑36 (MQ=255)
tCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTa < 1:4931725/36‑1 (MQ=255)
cTGGAGAAACTGAAAATATCGCTAAAAACGTTCTAc < 1:1683378/36‑1 (MQ=255)
cTGGAGAAACTGAAAATATCGCTAAAAACGGTCTAc > 1:2991814/1‑36 (MQ=255)
tGGAGAAACTGAAAATATCGCTAAAAACGTTCTACt > 1:6322484/1‑36 (MQ=255)
tGGAGAAACTGAAAATATCGCTAAAAACGTTCTACt < 1:595714/36‑1 (MQ=255)
tGGAGAAACTGAAAATATCGCTAAAAACGTTCTACt < 1:2467638/36‑1 (MQ=255)
ggAGAAACTGAAAATATCGCTAAAAACGTTCTACTg < 1:4181751/36‑1 (MQ=255)
ggAGAAACTGAAAATATCGCTAAAAACGTTCTACTg > 1:5306294/1‑36 (MQ=255)
agaAACTGAAAATATCGCTAAAAACGTTCTACTGcc > 1:5868849/1‑36 (MQ=255)
gaAACTGAAAATATCGCTAAAAACGTTCTACTGCCa > 1:372292/1‑36 (MQ=255)
gaAACTGAAAATATCGCTAAAAACGTTCTACTGCCa > 1:6396713/1‑36 (MQ=255)
aaaCTGAAAATATCGCTAAAAACGTTCTATTGCCAt > 1:1040020/1‑36 (MQ=255)
aaaCTGAAAATATCGCTAAAAACGTTCTACTGCCAt > 1:3241532/1‑36 (MQ=255)
aaaCTGAAAATATCGCTAAAAACGTTCTACTGCCAt > 1:7178199/1‑36 (MQ=255)
aaaCTGAAAATATCGCTAAAAACGTTCTACTGCCAt > 1:1435211/1‑36 (MQ=255)
aCTGAAAATATCGCTAAAAACGTTCTACTGCCATaa > 1:1044063/1‑36 (MQ=255)
aCTGAAAATATCGCTAAAAACGTTCTACTGCCATaa > 1:1897102/1‑36 (MQ=255)
aCTGAAAATATCGCTAAAAACGTTCTACTGACATaa > 1:2800067/1‑36 (MQ=255)
tGAAAATATCGCTAAAAACGTTCTACTGCCATAAac > 1:5822177/1‑36 (MQ=255)
tGAAAATATCGCTAAAAACGTTCTACTGCCATAAac > 1:6740144/1‑36 (MQ=255)
aaaaaaCATCCCAAAAACCGTTCAACTGCCATAAac < 1:2476702/34‑1 (MQ=12)
gAAAATATCGCTAAAAACGTTCTACTGCCATaaaaa > 1:3723184/1‑34 (MQ=255)
gAAAATATCGCAAAAAACGTTCTACTGCCATAAaca < 1:2731919/36‑1 (MQ=255)
aaaaTATCGCTAAAAACGTTCTACTGCCATAAacac > 1:5201676/1‑36 (MQ=255)
aaaaTATCGCTAAAAACGTTCTACTGCCATAAacac > 1:3478030/1‑36 (MQ=255)
aaaaTATCGCTAAAAACGTTCTACTGCCATAAacac > 1:2617281/1‑36 (MQ=255)
aaaaTATCCAAAAAAAGATTCCACTGCCATAAacac < 1:3973101/36‑1 (MQ=9)
aaaTATCGCTAAAAACGTTCTACTTCCATAAAcaca > 1:6780863/1‑36 (MQ=255)
aaaTATCGCTAAAAACGTTCTACTGCCATAAAcaca > 1:7171389/1‑36 (MQ=255)
aaaTATCGAAAAAAACGTTCTACTGCCATAAAcaca < 1:7616051/36‑1 (MQ=38)
aaTATCGAAAAAAACGTTCTACTGCCATAAACACaa < 1:5408286/36‑1 (MQ=38)
atatCGCTAAAAACGTTCTACTGCCAGAAACACAat > 1:7014364/1‑36 (MQ=255)
atCGCTAAAAACGTTCTACTGCCATAAACACAatat > 1:1590894/1‑36 (MQ=255)
atCGCTAAAAACGTTCTACTGCCATAAACACAatat > 1:1042337/1‑36 (MQ=255)
cGCTAAAAACGTTCTACTGCCATAAACACAATATCa > 1:6764257/1‑36 (MQ=255)
cGCTAAAAACGTTCTACTGCCATAAACACAATATCa > 1:2821032/1‑36 (MQ=255)
agCTAAAAAAGTTCTACTGAAATAAACACAATATCa > 1:196330/2‑36 (MQ=25)
gCTAAAAACGTTCTACTGCCATAAACACAATATCat > 1:6467182/1‑36 (MQ=255)
cTAAAAACGTTCTACTGCCATAAACACAATATCatg > 1:6090050/1‑36 (MQ=255)
cTAAAAACGTTCTACTGCCATAAACACAATATCatg > 1:7162272/1‑36 (MQ=255)
tAAAAACGTTCTACTGCCATAAACACAATATCatga > 1:2722985/1‑36 (MQ=255)
tAAAAACGTTCTACTGCCATAAACACAATATCatga > 1:5496776/1‑36 (MQ=255)
|
AACGGAAGATATTCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTACTGCCATAAACACAATATCATGA > REL606/3584569‑3584639
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 11 ≤
ATCG
< 24 ≤
ATCG
< 32 ≤
ATCG
< 37 ≤
ATCG
|