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 |
|
A4 F1 I1 R1
|
73 |
24.7 |
1166784 |
93.9% |
1095610 |
107.1 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NC_000913 |
635,896 |
A→G |
L225P (CTC→CCC) |
ybdN ← |
PAPS reductase‑like domain protein |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NC_000913 | 635,896 | 0 | A | G | 95.5%
| 68.2
/ ‑3.9
| 22 | L225P (CTC→CCC) | ybdN | PAPS reductase‑like domain protein |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (10/11); total (10/12) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01 |
CGATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAA > NC_000913/635770‑636029
|
cGATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTgatgat > 2:57867/1‑139 (MQ=255)
gATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGTCTTtgaggttt < 1:337321/139‑6 (MQ=255)
cGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCa > 2:320829/1‑139 (MQ=255)
tCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTc < 2:386699/139‑1 (MQ=255)
tGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTaaa < 1:466435/139‑1 (MQ=255)
gCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAAt > 1:131235/1‑139 (MQ=255)
tCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATgg < 1:57867/139‑1 (MQ=255)
aaGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCTCTTTCCAGTCGTAAATGGGGTAAATGTa < 1:320829/139‑1 (MQ=255)
aaGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTCGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa < 1:413053/139‑1 (MQ=255)
agagGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGTACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCg < 2:452992/139‑1 (MQ=255)
ctgTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATTGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACgttgt < 1:386497/137‑1 (MQ=255)
tGGGTTGCAGGGACTTTGCTGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCATGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGttta > 1:15852/1‑139 (MQ=255)
tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTc < 1:372201/55‑1 (MQ=255)
tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTc > 2:372201/1‑55 (MQ=255)
tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTa < 1:47479/72‑1 (MQ=255)
tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTa > 2:47479/1‑72 (MQ=255)
gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa > 2:351424/1‑77 (MQ=255)
gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa < 1:491059/77‑1 (MQ=255)
gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa < 1:351424/77‑1 (MQ=255)
gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa > 2:491059/1‑77 (MQ=255)
tGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTAt > 2:123267/1‑81 (MQ=255)
tGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTAt < 1:123267/81‑1 (MQ=255)
tGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTaa > 1:126747/1‑139 (MQ=255)
|
CGATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAA > NC_000913/635770‑636029
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A