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 F2 I208 R1
|
222 |
13.1 |
731276 |
96.6% |
706412 |
85.9 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
476,525 |
A→G |
intergenic (‑361/+185) |
tomB ← / ← acrB |
Hha toxicity modulator TomB/efflux RND transporter permease AcrB |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 476,525 | 0 | A | G | 100.0%
| 20.2
/ NA
| 7 | intergenic (‑361/+185) | tomB/acrB | Hha toxicity modulator TomB/efflux RND transporter permease AcrB |
| Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (3/4); total (3/4) |
AGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTT > NZ_CP009273/476439‑476610
|
aGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGgtta > 1:11952/1‑90 (MQ=255)
aGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGgtta > 1:48657/1‑90 (MQ=255)
aGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGgtta > 2:328150/1‑90 (MQ=255)
gATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGgttatt > 1:111000/1‑90 (MQ=255)
tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc < 2:260901/90‑1 (MQ=255)
tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc > 2:3876/1‑90 (MQ=255)
tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc < 2:72461/90‑1 (MQ=255)
gAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATaaa > 2:206135/1‑90 (MQ=255)
gaatAAAGAATTATCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATtaacta < 2:318494/90‑1 (MQ=255)
tAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCg < 2:248481/90‑1 (MQ=255)
aTGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGttt > 2:80835/1‑90 (MQ=255)
|
AGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTT > NZ_CP009273/476439‑476610
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
CACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATT > NZ_CP009273/476429‑476530
|
CACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTA > SRR3722090.12081/1‑100 (MQ=60)
CACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTA > SRR3722090.49177/1‑100 (MQ=60)
CTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATT > SRR3722090.112187/1‑100 (MQ=60)
|
CACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATT > NZ_CP009273/476429‑476530
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |