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 F5 I98 R1
|
10 |
14.2 |
839128 |
94.4% |
792136 |
86.3 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
859,587 |
G→A |
A164T (GCA→ACA) |
fsa → |
fructose‑6‑phosphate aldolase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 859,587 | 0 | G | A | 100.0%
| 43.7
/ NA
| 13 | A164T (GCA→ACA) | fsa | fructose‑6‑phosphate aldolase |
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (4/9); total (4/9) |
AGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGT > NZ_CP009273/859510‑859666
|
aGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCaa > 2:270396/1‑91 (MQ=255)
gggCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCaaa > 1:123887/1‑91 (MQ=255)
ggCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTc > 1:311454/1‑91 (MQ=255)
ggCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTc > 2:57869/1‑91 (MQ=255)
gCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCa < 1:264423/91‑1 (MQ=255)
gCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCa < 2:226426/91‑1 (MQ=255)
gACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTAc < 1:14091/91‑1 (MQ=255)
gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtgt > 2:50453/1‑91 (MQ=255)
aTTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGaa < 1:143793/91‑1 (MQ=255)
catgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTActc < 1:266235/91‑1 (MQ=255)
catgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTActc < 2:261548/91‑1 (MQ=255)
gcAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTg < 2:182444/91‑1 (MQ=255)
ggCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAt < 1:290475/91‑1 (MQ=255)
cGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAtgt < 2:145138/91‑1 (MQ=255)
|
AGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGT > NZ_CP009273/859510‑859666
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
TTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCT > NZ_CP009273/859501‑859686
|
TTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAA > SRR3722121.126807/1‑101 (MQ=60)
GGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTC > SRR3722121.318998/1‑101 (MQ=60)
GCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGAC < SRR3722121.270723/101‑1 (MQ=60)
GACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATG < SRR3722121.14399/101‑1 (MQ=60)
ATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTC < SRR3722121.147180/101‑1 (MQ=60)
CATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGA < SRR3722121.272576/101‑1 (MQ=60)
GGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAAC < SRR3722121.297481/101‑1 (MQ=60)
GACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGctgtctcttatacacatctccgagcccacgagactcctgagcatctc > SRR3722121.338453/1‑54 (MQ=60)
|
TTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACAACAGATGATTAGCT > NZ_CP009273/859501‑859686
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |