Sample Resequencing Stats
Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate |
Predicted Mutations |
Mean Coverage |
Total Reads |
Percent Mapped |
Mapped Reads |
Average Read Length |
A1 F3 I0 R1
|
560 |
56.4 |
4650532 |
94.9% |
4413354 |
59.3 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
RA |
W3110S.gb |
4,581,599 |
T→C |
10.9% |
K348E (AAA→GAA) |
mcrC ← |
5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | W3110S.gb | 4,581,599 | 0 | T | C | 10.9%
| 143.1
/ 8.6
| 55 | K348E (AAA→GAA) | mcrC | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
Reads supporting (aligned to +/- strand): ref base T (23/26); new base C (2/4); total (25/30) |
Fisher's exact test for biased strand distribution p-value = 6.78e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTA > W3110S.gb/4581529‑4581669
|
gCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt > 1:3756323/1‑71 (MQ=255)
tcctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt > 1:694125/1‑68 (MQ=255)
ctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgaga > 1:936136/1‑70 (MQ=255)
ctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgaga > 1:2275287/1‑70 (MQ=255)
atgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGatat > 1:2426245/1‑70 (MQ=255)
atgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATAtt > 1:2108631/1‑71 (MQ=255)
gatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa > 1:4333961/1‑71 (MQ=255)
ccGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa < 1:776889/71‑1 (MQ=255)
ccGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa < 1:4388709/71‑1 (MQ=255)
gCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAAt < 1:4116816/71‑1 (MQ=255)
cccGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTga > 1:3438904/1‑52 (MQ=255)
gATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATAtt > 1:945778/1‑55 (MQ=255)
aTGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:3752367/71‑1 (MQ=255)
aTGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:1776614/71‑1 (MQ=255)
aTGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:622682/71‑1 (MQ=255)
tgtgGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTc < 1:4494239/54‑1 (MQ=255)
tggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTaa < 1:2250166/70‑1 (MQ=255)
ggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa > 1:2271852/1‑52 (MQ=255)
ggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa > 1:778583/1‑52 (MQ=255)
gtCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCtt > 1:3540215/1‑71 (MQ=255)
tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGGAAATGTCGAGTAATTCtt < 1:3281450/70‑1 (MQ=255)
tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:2204341/70‑1 (MQ=255)
tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:649048/70‑1 (MQ=255)
tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:1213369/70‑1 (MQ=255)
gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTg > 1:494803/1‑69 (MQ=255)
gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt < 1:4550045/71‑1 (MQ=255)
gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt > 1:471925/1‑71 (MQ=255)
ggAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg < 1:1651094/71‑1 (MQ=255)
gAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg > 1:2595198/1‑70 (MQ=255)
gAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGt < 1:726730/71‑1 (MQ=255)
ggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata < 1:1098585/70‑1 (MQ=255)
gTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt < 1:4284900/64‑1 (MQ=255)
gTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt < 1:3374154/64‑1 (MQ=255)
aTTTTCTATAAAACTTACCGCTTATTTGAGATAtt > 1:2860220/1‑35 (MQ=255)
aTTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa > 1:1622165/1‑61 (MQ=255)
aTTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACg > 1:397517/1‑71 (MQ=255)
ttttCTATAAAACTTACCGCTTATTCGAGATATTCATCGaaaa > 1:3097301/1‑43 (MQ=255)
tttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATg < 1:146590/44‑1 (MQ=255)
ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:2362039/56‑1 (MQ=255)
ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa > 1:4354395/1‑58 (MQ=255)
ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata > 1:1183585/1‑64 (MQ=255)
ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:1096210/1‑70 (MQ=255)
ttCTATAAAACTTACCGCTTATTTGAGATATACATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:4419928/1‑70 (MQ=255)
ttCTATAAAACTTACCGCTTATTCGAGATATTCATCGaa < 1:1804939/39‑1 (MQ=255)
tCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCg > 1:3023943/1‑45 (MQ=255)
cTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAt > 1:2252905/1‑71 (MQ=255)
tataAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAtt > 1:285438/1‑71 (MQ=255)
ataAAACTTACCGCTTATTTGAGATATTCATCGaa < 1:3321150/35‑1 (MQ=255)
ataAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATaca < 1:556463/62‑1 (MQ=255)
taAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTc < 1:381531/70‑1 (MQ=255)
taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTaa < 1:163904/46‑1 (MQ=255)
aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa > 1:128944/1‑51 (MQ=255)
aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTg < 1:2349331/71‑1 (MQ=255)
aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTg < 1:4231436/71‑1 (MQ=255)
aaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata < 1:2749528/56‑1 (MQ=255)
aCCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:628796/70‑1 (MQ=255)
cTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAAtt < 1:4313888/70‑1 (MQ=255)
tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc > 1:419654/1‑55 (MQ=255)
tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACgg > 1:160713/1‑71 (MQ=255)
tGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTa > 1:4574057/1‑71 (MQ=255)
|
GCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTA > W3110S.gb/4581529‑4581669
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A