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 I200 R1
|
216 |
30.8 |
1683622 |
97.2% |
1636480 |
87.0 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
2,554,121 |
A→C |
N169T (AAC→ACC) |
BW25113_RS12740 → |
hypothetical protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,554,121 | 0 | A | C | 97.3%
| 118.7
/ ‑3.9
| 37 | N169T (AAC→ACC) | BW25113_RS12740 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (24/12); total (25/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 = 6.15e-01 |
CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAA > NZ_CP009273/2554035‑2554207
|
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac < 1:178431/90‑1 (MQ=255)
tCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACgg > 2:514823/1‑90 (MQ=255)
ggCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGc > 2:220188/1‑90 (MQ=255)
cacaGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTc > 2:566330/1‑90 (MQ=255)
tGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCa > 1:676318/1‑90 (MQ=255)
ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt > 1:324616/1‑90 (MQ=255)
ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt > 1:36408/1‑90 (MQ=255)
ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt > 1:403969/1‑90 (MQ=255)
tGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTc < 1:456658/90‑1 (MQ=255)
cAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCaa > 2:152499/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:12611/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:345015/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:395279/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:788820/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:751488/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:719475/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:50906/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:574587/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:280153/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:604821/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:690430/1‑90 (MQ=255)
ggAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGt > 2:444745/1‑90 (MQ=255)
gAACAATCGACAAACGAGCCAGGCTATCAAAAAccac > 1:471200/1‑37 (MQ=255)
gAACAATCGACAAACGAGCCAGGCTATCAAAAAccac < 2:471200/37‑1 (MQ=255)
gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTaaaaa > 2:561369/1‑46 (MQ=255)
gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTaaaaa < 1:561369/46‑1 (MQ=255)
gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTa > 1:761057/1‑90 (MQ=255)
gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTa > 1:185693/1‑90 (MQ=255)
aTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGaa < 1:604821/90‑1 (MQ=255)
tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg < 2:751488/90‑1 (MQ=255)
gACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGc < 1:597757/90‑1 (MQ=255)
cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta < 1:566330/90‑1 (MQ=255)
gAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATAc < 1:397052/90‑1 (MQ=255)
caaTCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg < 2:551332/88‑1 (MQ=255)
cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg < 2:324616/90‑1 (MQ=255)
cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg < 2:12611/90‑1 (MQ=255)
cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg < 1:12365/90‑1 (MQ=255)
aTCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGc < 2:302566/90‑1 (MQ=255)
aTCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGc > 1:191874/1‑90 (MQ=255)
tCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCa > 2:704447/1‑90 (MQ=255)
aaacCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGaaa < 1:364604/90‑1 (MQ=255)
|
CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAA > NZ_CP009273/2554035‑2554207
|
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
GTGCTCTTTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAACGA > NZ_CP009273/2554022‑2554220
|
GTGCTCTTTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAC < SRR3722079.716960/100‑1 (MQ=60)
CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAA < SRR3722079.179989/100‑1 (MQ=60)
gtctcgtgggctcggagatgtgtataagagacaGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGAT < SRR3722079.566822/67‑1 (MQ=60)
AACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCA > SRR3722079.683415/1‑100 (MQ=60)
CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT > SRR3722079.327224/1‑100 (MQ=60)
CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT > SRR3722079.36674/1‑100 (MQ=60)
CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT > SRR3722079.407389/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722079.12701/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722079.797271/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722079.759521/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722079.727159/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722079.282485/1‑100 (MQ=60)
TGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTG < SRR3722079.460655/100‑1 (MQ=60)
GTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAActgtctcttatacacatctccgagcccacgagaccagagagga > SRR3722079.475320/1‑57 (MQ=60)
GTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTA > SRR3722079.187314/1‑100 (MQ=60)
GTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTA > SRR3722079.769201/1‑100 (MQ=60)
ATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATAC < SRR3722079.610920/100‑1 (MQ=60)
GACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAG < SRR3722079.603754/100‑1 (MQ=60)
CAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGA < SRR3722079.571864/100‑1 (MQ=60)
GAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGC < SRR3722079.400376/100‑1 (MQ=60)
GAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGC > SRR3722079.193538/1‑100 (MQ=60)
AGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAA > SRR3722079.639408/1‑100 (MQ=60)
CTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGG < SRR3722079.12454/100‑1 (MQ=60)
ATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAG > SRR3722079.615882/1‑100 (MQ=60)
AAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAA < SRR3722079.367471/100‑1 (MQ=60)
CCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAACGA < SRR3722079.398580/100‑1 (MQ=60)
|
GTGCTCTTTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAACGA > NZ_CP009273/2554022‑2554220
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: ‑ |