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 I228 R1
|
214 |
26.5 |
1448192 |
97.1% |
1406194 |
86.4 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
1,174,665 |
T→G |
L206W (TTG→TGG) |
nagK → |
N‑acetylglucosamine kinase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,174,665 | 0 | T | G | 92.3%
| 35.6
/ ‑3.1
| 13 | L206W (TTG→TGG) | nagK | N‑acetylglucosamine kinase |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (7/5); total (8/5) |
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.34e-01 |
GCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTA > NZ_CP009273/1174591‑1174747
|
gCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAat > 1:166862/1‑90 (MQ=255)
gCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAat > 2:572551/1‑90 (MQ=255)
cTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgc < 1:287449/90‑1 (MQ=255)
tATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCaa < 2:398873/90‑1 (MQ=255)
gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGa < 1:338341/90‑1 (MQ=255)
gTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGc < 2:39300/90‑1 (MQ=255)
cGTGGCTGAATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGgcgc < 1:281865/90‑1 (MQ=255)
aaCACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGtt > 1:234994/1‑90 (MQ=255)
cacTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTAt > 2:679372/1‑90 (MQ=255)
tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg > 1:290060/1‑90 (MQ=255)
caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta > 2:104766/1‑90 (MQ=255)
caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta > 2:282591/1‑90 (MQ=255)
caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta > 2:591853/1‑90 (MQ=255)
|
GCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTA > NZ_CP009273/1174591‑1174747
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
GCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAA > NZ_CP009273/1174581‑1174764
|
GCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAAT > SRR3722112.168912/1‑100 (MQ=60)
CTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGAT < SRR3722112.290784/100‑1 (MQ=60)
GTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCA < SRR3722112.342492/100‑1 (MQ=60)
CGTGGCTGAATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCG < SRR3722112.285139/100‑1 (MQ=60)
TGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTT > SRR3722112.237804/1‑100 (MQ=60)
GTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTG > SRR3722112.293418/1‑100 (MQ=60)
CATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTT > SRR3722112.79322/1‑100 (MQ=60)
GGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAA < SRR3722112.257416/100‑1 (MQ=60)
|
GCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAA > NZ_CP009273/1174581‑1174764
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |