Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,756,588 | 0 | C→T | 11.0% | 4.8 ‑9.0 0.0 | 77 | 299 | ygbJ | predicted dehydrogenase, with NAD(P)‑binding Rossmann‑fold domain |
Rejected: Prediction has biased read strand distribution supporting polymorphism. | ||||||||||
Reads supporting (aligned to +/- strand): new base (10/0): ref base (24/43): total (34/43) | ||||||||||
Fisher's exact test strand distribution p-value = 1.20e-04 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
Only 100 of 105 total aligned reads displayed.
ACAACATCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGCAAGAAATTG > REL606/2756553‑2756623
|
acaacaTCTGAAACCCGGTACGGCAGTAA.TGGTttc < 1:2966140/36‑1 (MQ=255)
acaacaTCTGAAACCCGGTACGGCAGTAA.TGGTttc > 1:1619198/1‑36 (MQ=255)
acaacaTCTGAAACCCGGTACGGCAGGGA.TGGTttc > 1:4060988/1‑36 (MQ=37)
acaacaTCTGAAACCCGGTACGGCAGGCA.TTGTttc > 1:4870364/1‑36 (MQ=25)
caacaTCTGAAACCCGGTACGGCAGTAA.TGGTttct < 1:2583155/36‑1 (MQ=255)
caacaTCTGAAACCCGGTACGGCAGTAA.TGGTttct < 1:5480769/36‑1 (MQ=255)
aacaTCTGAACCCCGGTACGGCAGTAA.TGGTttctt < 1:3374522/36‑1 (MQ=255)
aacaTCTGAAACCCGGTACGGCAGTTA.TTGTttttt > 1:2491616/1‑36 (MQ=25)
aacaTCTGAAACCCGGTACGGCAGTAA.TGGTttctt < 1:291067/36‑1 (MQ=255)
aacaTCTGAAACCCGGTACGGCAGTAA.TGGTttctt < 1:1567986/36‑1 (MQ=255)
aacaTCTGAAACCCGGTACGGCAGCAA.TTGTttttt > 1:928114/1‑36 (MQ=25)
acaTCTGAAGCCCGTTACGGCCGTAA.TGGTttcttc < 1:1573884/36‑1 (MQ=25)
acaTCTGAAACCCGGTTCGGCAGTAA.TGGTttcttc < 1:5442449/36‑1 (MQ=255)
acaTCTGAAACCCGGTACGGCAGTTA.CGGTttcctt > 1:4460162/1‑33 (MQ=25)
acaTCTGAAACCCGGTACGGCAGTAA.TGGTttcttc < 1:5883705/36‑1 (MQ=255)
acaTCTGAAACCCGGTACGGCAGGAA.TGGGttcttc > 1:1332120/1‑36 (MQ=37)
caTCTGAACCCCGGTACGGCAGTAA.TGGTTTCTTcc < 1:3575310/36‑1 (MQ=255)
caTCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTcc < 1:3304065/36‑1 (MQ=255)
caTCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTcc < 1:1432108/36‑1 (MQ=255)
caTCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTcc < 1:5859887/36‑1 (MQ=255)
caTCTGAAACCCGGTACGGCAGTAA.TGGGttctttc > 1:4199695/1‑34 (MQ=37)
caTCTGAAACCCGGCACGGCAGTAA.TGGTTTTTTcc > 1:4451687/1‑36 (MQ=37)
aTCTGAAGCCCGGTACGGCAGTAA.TGGTTTCTTCCa < 1:508648/36‑1 (MQ=255)
aTCTGAAACCCGGTACGGCAGGAA.TTGTttttttcc > 1:2861463/1‑33 (MQ=16)
aTCTGAAACCCGGTACGGCAGGAA.TGGTttttttcc > 1:4045622/1‑33 (MQ=25)
tCTGAAACCCGGTACGGCAGTAA.TGGTTTTTTcccc > 1:5247699/1‑34 (MQ=38)
tCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCAc > 1:1427369/1‑36 (MQ=255)
tCTGAAACCCGGTACGGCAGGAAGTGGTTTTTTccc > 1:1297591/1‑35 (MQ=16)
tCTGAAACCCGGTACGGCAGGAA.AGGTtttttctcc > 1:4791044/1‑33 (MQ=16)
cTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACt > 1:4968746/1‑36 (MQ=255)
cTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACt < 1:1897044/36‑1 (MQ=255)
tGAAACCCGGTACGGCAGTAA.TGGTTTTTTCCCCTa > 1:749265/1‑36 (MQ=37)
tGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTa > 1:5869778/1‑36 (MQ=255)
tGAAACCCGGTACGGCAGGAA.TGGGTTTTTTCACtt > 1:1306062/1‑35 (MQ=16)
gAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTAt < 1:3043542/36‑1 (MQ=255)
gAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTAt < 1:4079393/36‑1 (MQ=255)
gAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTAt < 1:2513019/36‑1 (MQ=255)
caaCCAGGTCCGGCAGTAA.TGGTTTCTTCCACTATc < 1:5374225/35‑1 (MQ=37)
aaaCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATc < 1:2305724/36‑1 (MQ=255)
aaaCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATc < 1:53801/36‑1 (MQ=255)
aaaCCCGGTACGGCAGGAA.TGGTTTTTTTCTCTATc > 1:5054793/1‑36 (MQ=19)
aaCCCGGTACGGCAGTAA.TGGTTTTTTTCACTAttc > 1:670746/1‑34 (MQ=37)
aaCCCGGTACGGCAGTAA.TGGTTTCTTCCCCTATCg > 1:4541233/1‑36 (MQ=255)
aCCCGGTACGGCAGTAA.TGGTTTCTTCCTCTAtttc > 1:3108847/1‑33 (MQ=37)
aCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGc < 1:2555216/36‑1 (MQ=255)
aCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGc < 1:3912907/36‑1 (MQ=255)
aCCCGGTACGGCAGTAA.TGGGTTCTTACACTATCgg > 1:4130291/1‑35 (MQ=37)
cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt < 1:3067282/36‑1 (MQ=255)
cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt < 1:344390/36‑1 (MQ=255)
cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt < 1:1894954/36‑1 (MQ=255)
cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt < 1:5558458/36‑1 (MQ=255)
cccGGTACGGCAGGTA.TGGGTTTTTTCACCATCCCt > 1:4211464/1‑36 (MQ=0)
ccGGTACGTCAGTAA.TGGTTTCTTCCACTATCGCTa < 1:3719914/36‑1 (MQ=255)
ggTACGGCAGTTA.TGGTTTTTTTCCCTATCGCTTGt > 1:5539446/1‑36 (MQ=16)
gTACGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTg < 1:3650622/36‑1 (MQ=255)
aCGGCAGTAA.TGGTTTCTTCCACTATCTCTAGTGCt > 1:1861362/1‑36 (MQ=255)
aCGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCt < 1:681700/36‑1 (MQ=255)
cGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTg < 1:4016012/36‑1 (MQ=255)
cGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTg < 1:4451598/36‑1 (MQ=255)
cGGCAGGAA.TGGTTTCTTCCACTATCGCTAGTGCtc > 1:5365123/1‑35 (MQ=255)
gtcAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGa < 1:1591543/34‑1 (MQ=255)
gtcAGTAA.TGGTTTCTTCCACCATCGCTAGTGCTGa < 1:5419059/34‑1 (MQ=25)
ggCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGa < 1:5084423/36‑1 (MQ=255)
ggCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGa > 1:1151245/1‑36 (MQ=255)
gCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAt < 1:79208/36‑1 (MQ=255)
gCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAt > 1:3346845/1‑36 (MQ=255)
gCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAt > 1:5905124/1‑36 (MQ=255)
cAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATg < 1:2741762/36‑1 (MQ=255)
cAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATg < 1:598804/36‑1 (MQ=255)
cAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATg < 1:3965434/36‑1 (MQ=255)
aGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAttc > 1:1087336/1‑34 (MQ=255)
aGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgg > 1:4741808/1‑35 (MQ=255)
aGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgc < 1:4590110/36‑1 (MQ=255)
tAA.TGGTTTCTTCCACTCTCGCTAGTGCTGATgcgc < 1:246287/36‑1 (MQ=255)
tAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgcgc < 1:2884296/36‑1 (MQ=255)
tAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgcgc < 1:2637939/36‑1 (MQ=255)
tAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgcgc < 1:5647171/36‑1 (MQ=255)
aa.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGCa < 1:1606146/36‑1 (MQ=255)
aa.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGCa > 1:4288536/1‑36 (MQ=255)
aa.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGCa < 1:5084613/36‑1 (MQ=255)
a.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGCaa > 1:2138063/1‑36 (MQ=255)
a.TGGTTTCTTCCACGATCGCTAGTGCTGATGCGCaa > 1:3110139/1‑36 (MQ=255)
tGGTTTCTTCCACTATCGCTAGTGCTTATTCGCAAg > 1:2989184/1‑36 (MQ=37)
ggTTTCTTCCACTATCGCTATTGCTGATGCGCAAGa > 1:2431111/1‑36 (MQ=255)
ggTTTCTTCCACTATCGCTAGTGCTGATGCGCAAGa < 1:4656226/36‑1 (MQ=255)
ggTTTCTTCCAATATCGCTAGCGCTGATGCGAAAGa > 1:1153020/1‑36 (MQ=25)
gTTTCTTCCACTATCGCTAGTGCTGATGCGCAAGaa < 1:1057263/36‑1 (MQ=255)
gTTTCTTCCACTATCGCTAGTGCTGATGCGCAAGaa > 1:4702837/1‑36 (MQ=255)
gTTTCTTCCACTATCGCTAGTGCTGATGCGCAAGaa < 1:4439239/36‑1 (MQ=255)
tttCTTCCACTATCGCTAGTGCTGATGCGCAAGaaa < 1:2098967/36‑1 (MQ=255)
tttCTTCCACTATCGCTAGTGCTGATGCGCAAGaaa < 1:5064823/36‑1 (MQ=255)
tttCTTCCACTATCGCTAGTGATTATGTGCAAGaaa > 1:5031159/1‑36 (MQ=25)
ttcttcCACTATCGCTAGTGCTGATTCGCCAGAAAt > 1:3703759/1‑36 (MQ=37)
ttcttcCACTATCGCTAGTGCTGATGCGCAAGAAAt < 1:1337643/36‑1 (MQ=255)
ttcttcCACTATCGCTAGTGCTGATGCGCAAGAAAt > 1:2849521/1‑36 (MQ=255)
tcttcCACTATCGCTAGTGCTGATGCGCAAGAAAtt < 1:2369827/36‑1 (MQ=255)
tcttcCACTATCGCTAGTGCTGATGAGCAAGAAAtt > 1:4318213/1‑36 (MQ=255)
cttcCACTATCGCTAGTGCTGATGCGCAAGAAATTg < 1:2408456/36‑1 (MQ=255)
cttcCACTATCGCTAGTGCTGATGCGCAAGAAATTg < 1:2682147/36‑1 (MQ=255)
cttcCACTATCGCTAGTGCTGATGCGCAAGAAATTg > 1:3135416/1‑36 (MQ=255)
|
ACAACATCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGCAAGAAATTG > REL606/2756553‑2756623
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 7 ≤
ATCG
< 14 ≤
ATCG
< 27 ≤
ATCG
< 37 ≤
ATCG
|