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 F8 I0 R1
|
82 |
77.9 |
2158605 |
86.8% |
1873669 |
62.5 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
| RA |
minE |
617,556 |
T→A |
100% |
D198E (GAT→GAA) |
ycaL → |
predicted peptidase with chaperone function |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | minE | 617,556 | 0 | T | A | 86.4%
| 57.5
/ 3.0
| 22 | D198E (GAT→GAA) | ycaL | predicted peptidase with chaperone function |
| Reads supporting (aligned to +/- strand): ref base T (0/3); new base A (19/0); total (19/3) |
| Fisher's exact test for biased strand distribution p-value = 6.49e-04 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 |
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTGAAGAAACGTGGC > minE/617510‑617586
|
cTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTgaa < 1:993110/67‑1 (MQ=255)
cTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTgaa < 1:821171/67‑1 (MQ=255)
cTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTgaa < 1:1368766/67‑1 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1938560/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:823712/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:285726/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:2091175/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:2074938/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1990742/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1967839/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1102114/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1916056/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:177870/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1689714/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1609770/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1484677/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1422592/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1214750/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1170407/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1108946/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCGGTTGAAGAAACGTGGc > 1:796116/1‑66 (MQ=255)
aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCACCTTTGATCTGTTGAAGAAACGTGGc > 1:652938/1‑66 (MQ=255)
|
CTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTGAAGAAACGTGGC > minE/617510‑617586
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A