Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,887,465 | 0 | G→T | 5.4% | 0.8 ‑0.5 ‑0.1 | 55 | intergenic (‑121/‑38) | yobG/ECB_01797 | hypothetical protein/hypothetical protein |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (2/1): ref base (25/27): total (27/28) | ||||||||||
Fisher's exact test strand distribution p-value = 6.11e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.88e-01 |
TCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTG > REL606/1887431‑1887500
|
tCGGAATGAGTATCTTGCTTCAGCGCATATTGCCGc < 1:4690149/36‑1 (MQ=255)
tCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGc < 1:133686/36‑1 (MQ=255)
tCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGc < 1:3555087/36‑1 (MQ=255)
ggAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTa > 1:2375611/1‑36 (MQ=255)
ggAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTa < 1:2146516/36‑1 (MQ=255)
gcaTGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc < 1:2819168/34‑1 (MQ=255)
gAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc < 1:2895358/36‑1 (MQ=255)
gAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc > 1:957271/1‑36 (MQ=255)
gAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc < 1:4890631/36‑1 (MQ=255)
aaTGAGTAGCTTGCTTCAGCGCATATTGCCGCTACg < 1:5581277/36‑1 (MQ=255)
aaTGAGTAGCTTGCTTCAGCGCATATTGCCGCTACg > 1:909462/1‑36 (MQ=255)
aTGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGa > 1:2455569/1‑36 (MQ=255)
aTGAGGAGCTTGCTTCAGCGCATATTGCCGCTACGa > 1:556885/1‑36 (MQ=255)
tGAGTAGCTTGCTTCAGCGCATATTGCCTCTACGAt > 1:3639372/1‑36 (MQ=255)
tGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAt < 1:964170/36‑1 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt < 1:3873355/36‑1 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt < 1:3254669/36‑1 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt < 1:4368491/36‑1 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt < 1:3117778/36‑1 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt < 1:964796/36‑1 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt > 1:2314780/1‑36 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt > 1:1567971/1‑36 (MQ=255)
gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt < 1:1763877/36‑1 (MQ=255)
gTAGCTTGCTTCAGCGCATATTGCCTCTACGATTaa < 1:1474219/36‑1 (MQ=255)
gTAGCTTGCTTCAGCGCATATTGCCGCTACGATTaa < 1:4139745/36‑1 (MQ=255)
gTAGCTTGCTTCAGCGCATATTGCCGCTACGATTaa < 1:3451106/36‑1 (MQ=255)
gTAGCTTGCTTCAGCGCATATTGCCGCTACGATTaa > 1:2688278/1‑36 (MQ=255)
tAGCTTGCTTCAGCGCATATTGCCTCTACGATTAAg > 1:934551/1‑36 (MQ=255)
aGCTTGCTTCAGCGCATATTTCCGCTACGATTAAGc < 1:4636858/36‑1 (MQ=255)
gcttgcttCAGCGCATATTGCCGCTACGATTAAGCg > 1:5465423/1‑36 (MQ=255)
gcttgcttCAGCGCATATTGCCGCTACGATTAAGCg > 1:4254077/1‑36 (MQ=255)
gcttgcttCAGCGCATATTGCCGCTACGATTAAGCg > 1:4946310/1‑36 (MQ=255)
cttgcttCAGCGCATATTGCCGCTACGCTTAAGCGa > 1:1423096/1‑36 (MQ=255)
ttgcttCAGCGCATATTGCCGCTACGATTAAGCGaa > 1:1281292/1‑36 (MQ=255)
tgcttCAGCGCATATTGCCGCTACGATTAAGCGAAc < 1:2234931/36‑1 (MQ=255)
ttCAGCGCATATTGCCGCTACGATTAAGCGAACATa < 1:3442359/36‑1 (MQ=255)
ttCAGCGCATATGGCCGCTACGATTAAGCGAACATa < 1:5022686/36‑1 (MQ=255)
tCAGCGCATATTGCCGCTACGATTAAGCGAACATaa > 1:4953442/1‑36 (MQ=255)
cAGCGCATATTGCCGCTACGATTAAGCGAACATaaa < 1:526084/36‑1 (MQ=255)
gcgcATATTGCCGCTACGATTAAGCGAACATaaaaa > 1:4703504/1‑36 (MQ=255)
gcgaATATTGCCGCTACGATTAAGCGAACATaaaaa > 1:1088256/1‑36 (MQ=255)
cgcATATTGCCGCTACGATTAAGCGAACATAAAaag > 1:5352913/1‑36 (MQ=255)
gcATATTGCCGCTACGATTAAGCGAACATAAAaaga > 1:4825229/1‑36 (MQ=255)
gcATATTGCCGCTACGATTAAGCGAACATAAAaaga > 1:3804430/1‑36 (MQ=255)
atatTGCCGCTACGATTAAGCGAACATAAAaagaga > 1:4466694/1‑36 (MQ=255)
atatTGCCGCTACGATTAAGCGAACATAAAaagaga > 1:5428939/1‑36 (MQ=255)
atatTGCCGCTACGATTAAGCGAACATAAAaagaga < 1:3948888/36‑1 (MQ=255)
atatTGCCGCTACGATTAAGCGAACATAAAaagaga > 1:2070096/1‑36 (MQ=255)
tatTGCCGCTACGATTAAGCGAACATAAAaagagaa < 1:4438366/36‑1 (MQ=255)
tatTGCCGCTACGATTAAGCGAACATAAAaagagaa > 1:813416/1‑36 (MQ=255)
tatTGCCGCTACGATTAAGCGAACATAAAaagagaa < 1:2206669/36‑1 (MQ=255)
tatTGCCGCTACGATTAAGCGAACATAAAaagagaa > 1:1415270/1‑36 (MQ=255)
ttGCCGCTACGATTAAGCGAACATAAAaagagaaga > 1:5392211/1‑36 (MQ=255)
ttGCCGCTACGATTAAGCGAACATAAAaagagaaga < 1:3575979/36‑1 (MQ=255)
ttGCCGCTACGATTAAGCGAACATAAAaagagaaga < 1:805664/36‑1 (MQ=255)
gCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTg < 1:2218940/36‑1 (MQ=255)
|
TCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTG > REL606/1887431‑1887500
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 13 ≤
ATCG
< 22 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|