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 F14 I0 R1
|
248 |
61.6 |
2601574 |
98.7% |
2567753 |
146.9 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
RA |
NC_002947_CJ‑RC |
4,740,912 |
Δ1 bp |
100% |
intergenic (+42/+711) |
gltA → / ← yeiW |
citrate synthase/putative oxidoreductase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NC_002947_CJ‑RC | 4,740,912 | 0 | T | . | 93.3%
| 53.8
/ ‑1.1
| 15 | intergenic (+42/+711) | gltA/yeiW | citrate synthase/putative oxidoreductase |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (2/12); total (3/12) |
Fisher's exact test for biased strand distribution p-value = 2.00e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
Rejected as polymorphism: E-value score below prediction cutoff. |
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
GAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCA‑GGGTATTGAACGGCGCAT > NC_002947_CJ‑RC/4740771‑4741037
|
gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCCttttttat > 7:327504/1‑147 (MQ=255)
aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTGCGGGCGGCGTTTTTATttttac > 1:144493/1‑147 (MQ=255)
cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg < 1:86851/149‑1 (MQ=255)
aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGc < 5:251132/149‑1 (MQ=255)
gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc < 8:327504/149‑1 (MQ=255)
ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc < 3:127569/149‑1 (MQ=255)
cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg < 8:239801/149‑1 (MQ=255)
gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc < 7:289168/149‑1 (MQ=255)
cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa < 5:257935/148‑1 (MQ=255)
cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga < 1:274192/149‑1 (MQ=255)
gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag < 3:148827/149‑1 (MQ=38)
gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag < 8:114179/149‑1 (MQ=38)
gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag < 8:140720/149‑1 (MQ=255)
ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAgg < 2:187000/149‑1 (MQ=255)
cTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACATAGCCGCCTTTCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCAAAGTTTTTATACGCCCGCTAACCCACCCCATTGA‑gg > 8:273121/1‑149 (MQ=255)
aagaagGCTGCCCT‑CGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACATAGCGGCCGTTCCCGCGAAGAGGCCCGCCCCGACCCCACCAAAAACCACTTTTTGTGCGCCCCCGCACCCACCCCGATGACGGGTATTAAACgcgccaa > 8:313676/1‑143 (MQ=255)
|
GAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCT‑CGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCA‑GGGTATTGAACGGCGCAT > NC_002947_CJ‑RC/4740771‑4741037
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A