Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 3,328,657 | G→A | 11.6% | P108S (CCG→TCG) | mreB ← | cell wall structural complex MreBCD, actin‑like component MreB |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,328,657 | 0 | G→A | 11.6% | 6.7 ‑1.0 ‑0.1 | 43 | 322 | mreB | cell wall structural complex MreBCD, actin‑like component MreB |
Reads supporting (aligned to +/- strand): new base (4/1): ref base (21/17): total (25/18) | ||||||||||
Fisher's exact test strand distribution p-value = 3.80e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.37e-01 |
GGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTT > REL606/3328622‑3328692
|
ggTCGCGCCAACCGGCACACAAACCAGAACGCGCgg > 1:1795307/1‑36 (MQ=255)
cgcgCCAACCGGCACACAAACCAGAACGCGCGGGCt > 1:5470941/1‑36 (MQ=255)
gcgcCAACCGGCACACAAACCAGAACGCGCGGGCtt < 1:5725481/36‑1 (MQ=255)
cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTg > 1:2119324/1‑36 (MQ=255)
gcCAAACGGCACACAAACCAGAACGCGCGGGCTTgg > 1:6082736/1‑36 (MQ=255)
ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGa < 1:4990215/36‑1 (MQ=255)
ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGa > 1:403876/1‑36 (MQ=255)
ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGa > 1:3542606/1‑36 (MQ=255)
aaCCGGCACACAAACCAGAACGCGCGGGCTTGGACg > 1:2489126/1‑36 (MQ=255)
aCCGGCACACAAACCAGAACGCGCGGGCTTGGACGc < 1:5130182/36‑1 (MQ=255)
aCCGGCACACAAACCAGAACGCGCGAGCTTGGACGc > 1:211998/1‑36 (MQ=255)
ccGGCACACAAACCAGAACGCGCGGGCTTGGACGCa < 1:3772938/36‑1 (MQ=255)
ccGGCACACAAACCAGAACGCGCGGGCTTGGACGCa < 1:1022897/36‑1 (MQ=255)
acGGCACACAAACCAGAACGCGCGGGCTTGGACGCa > 1:7439000/2‑36 (MQ=255)
cGGCACACAAACCAGAACGCGCGAGCTTGGACGCAt > 1:1458314/1‑36 (MQ=255)
ggCACACAAACCAGAACGCGCGGGCTTGGACGCATa > 1:6951080/1‑36 (MQ=255)
ggCACACAAACCAGAACGCGCGAGCTTGGACGCATa > 1:3355141/1‑36 (MQ=255)
ggCACAAAAACCAGAACGCGCGGGCTTGGACGCATa > 1:6166301/1‑36 (MQ=255)
cacacaAACCAGAACGCGCGGGCTTGGACGCATaaa > 1:3611465/1‑36 (MQ=255)
aacacaAACCAGAACGCGCGGGCTTGGACGCATaaa < 1:5544331/35‑1 (MQ=255)
acacaAACCAGAACGCGCGGGCTTGGACGCATAAAg > 1:4210173/1‑36 (MQ=255)
acacaAACCAGAACGCGCGGGCTTGGACGCATAAAg > 1:7445058/1‑36 (MQ=255)
cacaAACCAGAACGCGCGGGCTTGGACGCATAAAGc < 1:5690676/36‑1 (MQ=255)
acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCt < 1:1198604/36‑1 (MQ=255)
acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCt < 1:5255795/36‑1 (MQ=255)
acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCt > 1:5809425/1‑36 (MQ=255)
caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTg > 1:678080/1‑36 (MQ=255)
aaaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGt < 1:3160869/36‑1 (MQ=255)
aaaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGt > 1:1266324/1‑36 (MQ=255)
aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGtt > 1:1266661/1‑36 (MQ=255)
aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGtt > 1:3955541/1‑36 (MQ=255)
aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGtt < 1:1918724/36‑1 (MQ=255)
aaCCAGAACGCGCGAGCTTGGACGCATAAAGCTGtt < 1:6675898/36‑1 (MQ=255)
cAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCt < 1:1192542/36‑1 (MQ=255)
gAACGCGCGGGCTTGGACGCATAAAGCTGTTGCtgt > 1:2340392/1‑36 (MQ=255)
aaCGCGCGGGCTTGGACGCATAAAGCTGTTGCtgtg < 1:4987362/36‑1 (MQ=255)
aCGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGc < 1:441703/36‑1 (MQ=255)
aCGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGc < 1:6892817/36‑1 (MQ=255)
cgcgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCa > 1:528843/1‑36 (MQ=255)
gcgcgAGCTTGGACGCATAAAGCTGTTGCTGTGCAc > 1:2864785/1‑36 (MQ=255)
cgcgGGCTTGGACGCATAAAGCTGTTTCTGTGCCCt < 1:5496209/36‑1 (MQ=37)
cgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACt > 1:3855699/1‑36 (MQ=255)
cgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACt < 1:2136337/36‑1 (MQ=255)
gcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACtt > 1:5524428/1‑36 (MQ=255)
gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt > 1:5107378/1‑36 (MQ=255)
gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt > 1:5028427/1‑36 (MQ=255)
gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt < 1:6134219/36‑1 (MQ=255)
gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt < 1:1544599/36‑1 (MQ=255)
gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt > 1:293831/1‑36 (MQ=255)
gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt < 1:2063131/36‑1 (MQ=255)
ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt < 1:631009/36‑1 (MQ=255)
ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt < 1:3728441/36‑1 (MQ=255)
ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt < 1:2795548/36‑1 (MQ=255)
ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt < 1:2640500/36‑1 (MQ=255)
|
GGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTT > REL606/3328622‑3328692
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 17 ≤
ATCG
< 26 ≤
ATCG
< 32 ≤
ATCG
< 37 ≤
ATCG
|