Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,119,132 | +G | 7.4% | coding (117/456 nt) | csgA → | cryptic curlin major subunit |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,119,132 | 1 | .→G | 7.4% | 3.1 ‑0.7 ‑0.2 | 27 | coding (117/456 nt) | csgA | cryptic curlin major subunit |
Reads supporting (aligned to +/- strand): new base (2/0): ref base (15/10): total (17/10) | ||||||||||
Fisher's exact test strand distribution p-value = 5.16e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.19e-01 |
CGGCGGTAACCACGGTGGTGGCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTACG > REL606/1119099‑1119166
|
cggcggTAACCACGGTGGTGGCGGTAATAATAGC.gg < 1:312173/36‑1 (MQ=255)
gcggTAACCACGGTGGTGGCGGTAATAATAGC.gggc > 1:1197190/1‑34 (MQ=255)
cggTAACCACGGTGGTGGCGGGAATAATAGC.GG.ccc > 1:2823194/1‑36 (MQ=255)
gTAACCACGGTGGTGGCGGTAATAATAGC.GG.gccca > 1:7933847/1‑34 (MQ=38)
aaCCACGGTGGTGGCGGTAATAATAGA.GG.GCCCAAt > 1:1431114/1‑36 (MQ=25)
aCCACGGTGGTGGCGGTAATAATAGC.GG.CCCAAAtt < 1:4503182/36‑1 (MQ=255)
ccACGGTGGTGGCGGTAATAATAGC.GG.CCCCAATTc > 1:58611/1‑36 (MQ=255)
cACGGTGGTGGCGGTGATAATAGT.GGCCCCAAAttt > 1:5338334/1‑35 (MQ=16)
cACGGTGGTGGCGGGAATAATATC.GG.CCCCAAtttt > 1:6373282/1‑34 (MQ=21)
cACGGTGGTGGCGGAAATAATAGCGGG.CCCAAAtct > 1:8231317/1‑34 (MQ=21)
aCGGTGGTGGCGGTAATAATAGC.GG.CCCAAATTCTg < 1:3597947/36‑1 (MQ=255)
aCGGTGGTGGCGGTAATAATAACGGG.CCCAAAtttt > 1:3715826/1‑34 (MQ=21)
aCGGTGGTGGCGGATATAATAGC.GG.CCCGCATTcgt > 1:201586/1‑34 (MQ=16)
ggtggtGGCGGTAATAATAGC.GG.CCCAAATTCTGAg < 1:6362414/36‑1 (MQ=255)
tggtggCGGTAATAATAGC.GG.CCCAAATTTTGAGCt > 1:8291967/1‑36 (MQ=255)
ggtggCGGTAATAATAGC.GG.CCCCAATTCTGAGCTg > 1:1269216/1‑36 (MQ=255)
gggCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGaa < 1:4615320/35‑1 (MQ=255)
gCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACa > 1:4047429/1‑36 (MQ=255)
gCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACa < 1:2581573/36‑1 (MQ=255)
cGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACAt > 1:1311964/1‑36 (MQ=255)
ggTAATAATAGC.GG.CCCAAATTCTGAGCTGAACAtt < 1:6043328/36‑1 (MQ=255)
ggTAATAATAGC.GG.CCCAAATTCTGAGCTGAACAtt < 1:6225061/36‑1 (MQ=255)
taataaTAGC.GG.CCCAAATTCTGAGCTGAACATTTa < 1:5058543/36‑1 (MQ=255)
taataaTAGC.GG.CCCAAATTCTGAGCTGAACATTTa < 1:4353484/36‑1 (MQ=255)
ataataGC.GG.CCCCAATTCTGAGCTGAACATTTAcc > 1:7837103/1‑36 (MQ=255)
aataGC.GG.CCCAAATTCTGAGCTGAACATTTACCAg > 1:3717395/1‑36 (MQ=255)
taGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTa > 1:3266733/1‑36 (MQ=255)
aGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTAc < 1:7418774/36‑1 (MQ=255)
gC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTACg > 1:6069622/1‑36 (MQ=255)
|
CGGCGGTAACCACGGTGGTGGCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTACG > REL606/1119099‑1119166
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 6 ≤
ATCG
< 13 ≤
ATCG
< 27 ≤
ATCG
< 37 ≤
ATCG
|