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 I1 R1
|
748 |
32.0 |
1685426 |
92.6% |
1560704 |
103.9 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NC_000913 |
4,424,386 |
G→A |
*229* (TAG→TAA) |
ulaF → |
L‑ribulose 5‑phosphate 4‑epimerase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NC_000913 | 4,424,386 | 0 | G | A | 95.5%
| 60.0
/ ‑3.6
| 22 | *229* (TAG→TAA) | ulaF | L‑ribulose 5‑phosphate 4‑epimerase |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (9/12); total (10/12) |
| Fisher's exact test for biased strand distribution p-value = 4.55e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01 |
AGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTAT > NC_000913/4424260‑4424507
|
agaagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGc > 1:324987/1‑139 (MQ=255)
ggATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga < 1:305460/139‑1 (MQ=255)
tGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATa < 2:324987/139‑1 (MQ=255)
cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 1:348817/1‑139 (MQ=255)
cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 1:646086/1‑139 (MQ=255)
caACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc < 1:25328/122‑1 (MQ=255)
caACTCAATCACATCGACAGCTTCCTGATGAAAAAACACTTCATGCGTAAACACGGTCCGAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 2:25328/1‑122 (MQ=255)
cTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGccc < 2:646086/139‑1 (MQ=255)
cTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCCGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag < 1:123655/139‑1 (MQ=255)
ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaa < 2:490906/111‑1 (MQ=255)
ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaa > 1:490906/1‑111 (MQ=255)
tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca < 1:287073/139‑1 (MQ=255)
tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca > 2:475742/1‑139 (MQ=255)
gCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACt > 1:136915/1‑139 (MQ=255)
cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg < 2:348817/139‑1 (MQ=255)
cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg < 1:475742/139‑1 (MQ=255)
aaCGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 2:344423/1‑62 (MQ=255)
aaCGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc < 1:344423/62‑1 (MQ=255)
gCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt < 2:710079/139‑1 (MQ=255)
tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCg < 2:352417/47‑1 (MQ=255)
tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCg > 1:352417/1‑47 (MQ=255)
tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat > 1:483983/1‑139 (MQ=255)
|
AGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTAT > NC_000913/4424260‑4424507
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A