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 F8 I0 R2
|
422 |
46.8 |
3726540 |
94.7% |
3529033 |
61.7 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
| RA |
W3110S.gb |
2,501,495 |
G→A |
100% |
S135S (TCG→TCA) |
ddg → |
palmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | W3110S.gb | 2,501,495 | 0 | G | A | 87.5%
| 79.3
/ 5.0
| 32 | S135S (TCG→TCA) | ddg | palmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (1/3); new base A (9/19); total (10/22) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01 |
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACATAAT > W3110S.gb/2501430‑2501561
|
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTgg < 1:1703996/70‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTgg < 1:2754895/70‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa < 1:902784/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa < 1:1719010/71‑1 (MQ=255)
aCCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACt < 1:3570427/70‑1 (MQ=255)
ttAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGggc > 1:3423192/1‑71 (MQ=255)
aaaCGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGggcg < 1:83082/70‑1 (MQ=255)
cgcACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACt < 1:1038624/60‑1 (MQ=255)
gcACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGt < 1:942785/71‑1 (MQ=255)
aaaTGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATg > 1:1000736/1‑71 (MQ=255)
gTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat < 1:3277750/71‑1 (MQ=255)
tAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat < 1:3303343/70‑1 (MQ=255)
aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatgat > 1:2628096/1‑71 (MQ=255)
aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatgat > 1:749727/1‑71 (MQ=255)
aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatgat > 1:1715156/1‑71 (MQ=255)
aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatga > 1:1643787/1‑70 (MQ=255)
aTGGATGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGGCAACCAatgct > 1:708437/1‑69 (MQ=255)
gTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTa < 1:1668088/70‑1 (MQ=255)
cGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATgg < 1:2670228/41‑1 (MQ=255)
gTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat < 1:936895/56‑1 (MQ=255)
gTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat < 1:2483506/56‑1 (MQ=255)
gTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat < 1:1721976/56‑1 (MQ=255)
tCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGccaacc < 1:2034521/52‑1 (MQ=255)
tCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGccaacc < 1:258681/52‑1 (MQ=255)
cATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTc > 1:2976872/1‑70 (MQ=255)
cATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTc > 1:3704029/1‑70 (MQ=255)
ttttATGTCACTGGAACTGGGCGGCCGCGTGATgg < 1:2037899/35‑1 (MQ=255)
ttttATGTCACTGGAACTGGGCGGCCGCGTGATgg < 1:1783406/35‑1 (MQ=255)
ttttATGTCACTGGAACTGGGCGGCCGCGTGATgg < 1:672379/35‑1 (MQ=255)
ttATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCcaca > 1:2980910/1‑70 (MQ=255)
tGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACataat < 1:379510/71‑1 (MQ=255)
tGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACataat < 1:2908696/71‑1 (MQ=255)
|
GATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACATAAT > W3110S.gb/2501430‑2501561
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A