Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,602,297 | C→A | 21.3% | D68Y (GAT→TAT) | ECB_01505 ← | conserved hypothetical protein |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,602,297 | 0 | C→A | 21.3% | 4.9 ‑4.6 ‑0.0 | 32 | 202 | ECB_01505 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): new base (1/6): ref base (18/7): total (19/13) | ||||||||||
Fisher's exact test strand distribution p-value = 1.02e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01 |
CCATCAGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCGCCCCACCACGGCATGTTAATCCGCTTGC > REL606/1602264‑1602332
|
ccATCAGGCTTGCCACATCCCAAACA.CGCGGATCaa < 1:3785803/36‑1 (MQ=255)
ccATCAGGCTTGCCACATCCCAAACA.CGCGGATCaa < 1:7004702/36‑1 (MQ=255)
cAGGCTTGCCACATCCCAAACA.CGCGGATCAAAAcc > 1:5452131/1‑36 (MQ=255)
cAGGCTTGCCACATCCCAAACA.CGCGGATAAAAAcc > 1:5225868/1‑36 (MQ=255)
aGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCg > 1:4332580/1‑36 (MQ=255)
aGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCg < 1:3894868/36‑1 (MQ=255)
aGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCg < 1:1577844/36‑1 (MQ=255)
ggCTTGCCACATCCCAAACA.CGCGGATCAAAACCGc < 1:4938078/36‑1 (MQ=255)
cTTGCCACATCCCAAACA.CGCGGATCAAAACCGccc > 1:6999892/1‑36 (MQ=255)
ttGCCACATCCCAAACA.CGCGGATCAAAACCGcccc > 1:7068346/1‑36 (MQ=255)
tGCCACATCCCAAACA.CGCGGATCAAAACCGCCcca > 1:3708941/1‑36 (MQ=255)
gCCACATCCCAAACA.CGCGGATCAAAACCGCCacac > 1:7750290/1‑36 (MQ=255)
gCCACACCCCAAACA.CGCGGATCAAAACCGCCccac < 1:2654754/36‑1 (MQ=255)
gCCACACCCCAAACA.CGCGAATCAAAACCGCCccac < 1:397183/36‑1 (MQ=18)
cacaTCCCAAACA.CGCGGATCAAAACCGCCccacca > 1:7931161/1‑36 (MQ=255)
cacaTCCCAAACA.CGCGGATCAAAACCCCCccacca > 1:2122270/1‑36 (MQ=255)
caTCCCAAACA.CGCGGATCAAAACCGCCCCACCACg > 1:4714562/1‑36 (MQ=255)
caTCCCAAACA.CGCGGATCAAAACCGCCCAACCACg > 1:3237960/1‑36 (MQ=255)
cacccAAAAACCGCGGAAAAAAACCGCCCCACCACg < 1:3567183/34‑1 (MQ=0)
aTCCCAAACA.CGCGGATCAAAACCGCCCCACCACgg > 1:7496004/1‑36 (MQ=255)
tCCCAAACA.CGCGGATCAAAACCTTCCCACCACGGc > 1:731508/1‑36 (MQ=18)
tCCCAAACA.CGCGGATCAAAACCGCCCCACCACGGc > 1:7539352/1‑36 (MQ=255)
ccaaaaCA.CCCGAATAAAAACCGCCCCACCACGGCa < 1:6356746/33‑1 (MQ=9)
caaaaCA.CGCGAATAAAAACCGCCCCACCACGGCAt < 1:2714389/34‑1 (MQ=17)
aaCA.CGCGGATCAAAACCGCCACAACACGGCATGtt > 1:5573608/1‑36 (MQ=18)
aca.cGCGGATAAAAACCGCCCCACCACGGCATGTTa < 1:6687436/36‑1 (MQ=255)
ca.cGCGAATAAAAACCGCCCCACCACGGCATGTTaa < 1:4806137/36‑1 (MQ=18)
cgcgGATCAAAACCGCCCCACCACGGCATGTTAATc > 1:5384374/1‑36 (MQ=255)
cgcgGATAAAAACCGCCCCACCACGGCATGTTAATc < 1:4064734/36‑1 (MQ=255)
cgGATCAAAACCGCCCCACCACGGCATGTTAATCCg > 1:3085316/1‑36 (MQ=255)
gATCAAAACCGCCCCACCACGGCATGTTAATCCGCt > 1:240397/1‑36 (MQ=255)
aTCAAAACCGCCCCACCACGGCATGTTAATCCGCtt > 1:225556/1‑36 (MQ=255)
acAAACCCGCCCCACCACGGCATGTTAATCCGCTTg < 1:3159940/35‑1 (MQ=255)
cAAAACCGCCCCACCACGGCATGTTAATCCGCTTGc > 1:1344541/1‑36 (MQ=255)
|
CCATCAGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCGCCCCACCACGGCATGTTAATCCGCTTGC > REL606/1602264‑1602332
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 9 ≤
ATCG
< 19 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|