Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 3,453,694 | +G | 6.1% | coding (2142/2553 nt) | mrcA → | fused penicillin‑binding protein 1a: murein transglycosylase/murein transpeptidase |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,453,694 | 1 | .→G | 6.1% | 3.7 ‑1.6 ‑0.7 | 33 | coding (2142/2553 nt) | mrcA | fused penicillin‑binding protein 1a: murein transglycosylase/murein transpeptidase |
Reads supporting (aligned to +/- strand): new base (2/0): ref base (13/18): total (15/18) | ||||||||||
Fisher's exact test strand distribution p-value = 1.99e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.09e-01 |
CGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAGATG > REL606/3453659‑3453728
|
cGTGCAGGTCGTGATTTGCGGCGTCGCGATATcggc > 1:4597351/1‑36 (MQ=255)
cGTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc > 1:4029498/1‑36 (MQ=255)
cGTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc > 1:3773878/1‑36 (MQ=255)
gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g > 1:6843295/1‑36 (MQ=255)
gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g < 1:2502913/36‑1 (MQ=255)
gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g > 1:2123351/1‑36 (MQ=255)
gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g < 1:1464704/36‑1 (MQ=255)
gTGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.g > 1:1544550/1‑36 (MQ=255)
tGCAGGTCGTGATTTGCAGCGTCGCGATATcggc.gg > 1:2751056/1‑36 (MQ=255)
gCAGGTCGTGATTTGCAGCGTCGCGATATCGGC.ggg > 1:1611453/1‑36 (MQ=255)
ggTCGTGATTTGCAGCGTCGCGATATCGGC.ggggga > 1:6504121/1‑33 (MQ=255)
gTCTTGATTTGCACCGCCGCGATATCGGC.GGGaaaa < 1:5680035/36‑1 (MQ=25)
gTCGTGATTTGCAGCGTCGCGATATCGGC.GGGgaaa > 1:2455643/1‑36 (MQ=255)
tCGTGATTTGCAGCGTCGCGATATCGGC.GGGGAAAc > 1:8140820/1‑36 (MQ=255)
tCGTGATTTGCAGCGTCGCGATATCGGC.GGGAAAAc < 1:8336320/36‑1 (MQ=255)
cgggATTGGCAGTGTCGCGATATCTGC.GGGAAAAcc < 1:6329713/33‑1 (MQ=19)
cGTGATTTGCATCGTCTCGATATCGGC.GGGAAAAcc < 1:3309908/36‑1 (MQ=37)
cGTGATTTGCAGCGTCGCGATATCGGC.GGGGAAAcc > 1:4383416/1‑36 (MQ=255)
gTGATTTGCAGCGTCGCGATATCGGCGGGGAAAAcc > 1:5549352/1‑36 (MQ=37)
gTGATTTGCAGCGTCGCGATATCGGC.GGGAAAACCg < 1:3265508/36‑1 (MQ=255)
tGATTTGCAGCGTCGCGATATCGGCGGGGAAAACCg > 1:809379/1‑36 (MQ=37)
tGATTTGCAGCGTCGCGATATCGGC.GGGGAAACCgg > 1:3081895/1‑36 (MQ=255)
aTTTGCAGCGTCGCGATATCGGC.GGGGAAACCCGTa > 1:6559921/1‑36 (MQ=37)
aTTTGCAGCGTCGCGATATCGGC.GGGAAAACCGGTa < 1:4723115/36‑1 (MQ=255)
gCAGCGTCGCGATATCGGC.GGGAAAACCGGTACCAc < 1:2649419/36‑1 (MQ=255)
cAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACt < 1:5227/36‑1 (MQ=255)
cAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACt < 1:5018652/36‑1 (MQ=255)
gCGTCGCGATATCGGC.GGGGAAACCGGTACCACTaa > 1:8443617/1‑36 (MQ=255)
cGTCGCGATATCGGC.GGGAAAACCGGTACCACTAAc > 1:718048/1‑36 (MQ=255)
cGTCGCGATATCGGC.GGGAAAACCGGTACCACTAAc < 1:7515203/36‑1 (MQ=255)
cGTCGCGATATCGGC.GGGAAAACCGGGACCCCTCAc > 1:3600675/1‑36 (MQ=25)
gTCGCGATATCGGC.GGGGAAACCGGTACCACTCACa > 1:6927227/1‑36 (MQ=37)
gTCGCGATATCGGC.GGGAAAACCGGTACCACTAACa < 1:5770365/36‑1 (MQ=255)
gTCGCGATATCGGC.GGGAAAACCGGTACCACTAACa < 1:8133561/36‑1 (MQ=255)
gTCGCGATATCGGC.GGGAAAAACGGTACCACTAAcc > 1:8195788/1‑35 (MQ=255)
tCGCGATATCGGC.GGGAAAGCCGGTACCACTAACAg < 1:2325150/36‑1 (MQ=255)
tCGCGATATCGGC.GGGAAAACCGGTACCACTAACAg < 1:1470525/36‑1 (MQ=255)
cgATATCGGC.GGGAAAACCGGTACCACTAACAGTTc < 1:4826498/36‑1 (MQ=255)
cgATATCGGC.GGGAAAACCGGTACCACTAACAGTTc < 1:4634778/36‑1 (MQ=255)
gATATCGGC.GGGAAAACCGGTACCACTAACAGTTCg < 1:561442/36‑1 (MQ=255)
atatCGGC.GGGAAAAGCGGTACCACTAACAGTTCGa < 1:4389707/36‑1 (MQ=255)
atatCGGC.GGGAAAACCGGGACCACTAACAGTTCgt > 1:1158358/1‑35 (MQ=255)
tagcggc.ggGAAAACCGGTACCACTAACAGTTCGaa < 1:7719420/33‑1 (MQ=255)
tCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAg > 1:3680526/1‑36 (MQ=255)
tCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAg > 1:957362/1‑36 (MQ=255)
cggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGa < 1:7289805/36‑1 (MQ=255)
ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt > 1:543395/1‑36 (MQ=255)
ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt < 1:4065396/36‑1 (MQ=255)
ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt < 1:1756027/36‑1 (MQ=255)
ggc.ggGAAAACCGGTACCACTAACAGTTCGAAAGAt > 1:8441554/1‑36 (MQ=255)
ggc.ggGAAAACCGGTACCACCAACAGTTCGAAAGAt > 1:4312749/1‑36 (MQ=255)
gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg < 1:6772383/36‑1 (MQ=255)
gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg > 1:4938494/1‑36 (MQ=255)
gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg > 1:7345659/1‑36 (MQ=255)
gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg < 1:3956419/36‑1 (MQ=255)
gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg < 1:1800639/36‑1 (MQ=255)
gc.ggGAAAACCGGTACCACTAACAGTTCGAAAGATg > 1:867785/1‑36 (MQ=255)
|
CGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC.GGGAAAACCGGTACCACTAACAGTTCGAAAGATG > REL606/3453659‑3453728
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 20 ≤
ATCG
< 29 ≤
ATCG
< 37 ≤
ATCG
|