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 |
|
A2 F1 I0 R1
|
4577 |
57.4 |
3436498 |
77.7% |
2670158 |
63.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
| RA |
minE |
1,532,228 |
A→T |
100% |
D166E (GAT→GAA) |
aroC ← |
chorismate synthase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | minE | 1,532,228 | 0 | A | T | 95.3%
| 43.4
/ ‑4.1
| 21 | D166E (GAT→GAA) | aroC | chorismate synthase |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (0/20); total (1/20) |
| Fisher's exact test for biased strand distribution p-value = 4.76e-02 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
| Rejected as polymorphism: E-value score below prediction cutoff. |
| Rejected as polymorphism: Frequency below/above cutoff threshold. |
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCG > minE/1532206‑1532288
|
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2528438/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:461837/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:384659/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:3433037/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:3195152/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2996289/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2994602/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2938176/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:277767/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2646980/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:1098693/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2307573/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:1796947/69‑1 (MQ=255)
cGACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:1785016/69‑1 (MQ=255)
gACCTGAGACCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCaa < 1:2470535/68‑1 (MQ=255)
aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt < 1:1585249/68‑1 (MQ=255)
aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt < 1:1511539/68‑1 (MQ=255)
aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt < 1:3226501/68‑1 (MQ=255)
aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt < 1:1375292/68‑1 (MQ=255)
aCCAGTCTTTGATTTCCAGTGGAATGTCGCCCATCTGGGTCAGGCAACCGCGAATCTCAATGCCAAAt < 1:1358441/68‑1 (MQ=255)
tGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCg > 1:1286996/1‑65 (MQ=255)
|
CGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCG > minE/1532206‑1532288
|
| 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