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 |
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.8%
| 154.2
/ ‑3.7
| 46 | N169T (AAC→ACC) | BW25113_RS12740 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base C (24/21); total (24/22) |
Fisher's exact test for biased strand distribution p-value = 4.78e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGG > NZ_CP009273/2554035‑2554201
|
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac < 2:296301/90‑1 (MQ=255)
cTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAccac < 2:539119/90‑1 (MQ=255)
tCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATGAAAAACCACgg > 2:113155/1‑90 (MQ=255)
tCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACgg < 2:533382/90‑1 (MQ=255)
gCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGtt > 1:579716/1‑90 (MQ=255)
ggCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGc > 2:155760/1‑90 (MQ=255)
aaaCTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGa > 1:429802/1‑90 (MQ=255)
acaGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCt < 1:146468/90‑1 (MQ=255)
ggATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAc < 2:373710/90‑1 (MQ=255)
ggATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAc < 1:551984/90‑1 (MQ=255)
ggATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAc < 1:313786/90‑1 (MQ=255)
gATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAcc < 1:334541/90‑1 (MQ=255)
gATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAcc < 2:561472/90‑1 (MQ=255)
ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt > 1:701341/1‑90 (MQ=255)
ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt > 1:677828/1‑90 (MQ=255)
ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt > 1:226232/1‑90 (MQ=255)
gggTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAc > 2:8199/1‑90 (MQ=255)
gggTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAc > 1:171851/1‑90 (MQ=255)
gggTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAc > 1:184966/1‑90 (MQ=255)
tGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgt > 2:197323/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:544979/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:288029/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:209198/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:509028/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:116889/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 2:105187/1‑90 (MQ=255)
ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt > 1:313503/1‑90 (MQ=255)
gTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgttg < 1:111856/90‑1 (MQ=255)
gTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgttg < 1:509028/90‑1 (MQ=255)
aGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAg < 1:8199/90‑1 (MQ=255)
gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAcc > 2:588709/1‑71 (MQ=255)
gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAAcc < 1:588709/71‑1 (MQ=255)
gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCg > 1:33831/1‑90 (MQ=255)
gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCg > 2:396704/1‑90 (MQ=255)
aaTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGa > 1:101096/1‑90 (MQ=255)
aaTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGa > 1:312460/1‑90 (MQ=255)
tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg < 2:378707/90‑1 (MQ=255)
tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg < 2:116889/90‑1 (MQ=255)
tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg < 2:33831/90‑1 (MQ=255)
tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg < 2:184966/90‑1 (MQ=255)
cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta < 2:444759/90‑1 (MQ=255)
cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta < 2:544979/90‑1 (MQ=255)
aaaCGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCtat < 1:288029/90‑1 (MQ=255)
aaaCGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCtat < 1:197323/90‑1 (MQ=255)
aaaCGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCtat < 1:49962/90‑1 (MQ=255)
cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg < 2:312460/90‑1 (MQ=255)
cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg < 1:87056/90‑1 (MQ=255)
tATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCgg > 1:39994/1‑90 (MQ=255)
|
CTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGG > NZ_CP009273/2554035‑2554201
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
TTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAACGA > NZ_CP009273/2554029‑2554220
|
TTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTT > SRR3722112.588560/1‑100 (MQ=60)
GTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGA > SRR3722112.435619/1‑100 (MQ=60)
ACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCAC < SRR3722112.148284/100‑1 (MQ=60)
CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT > SRR3722112.688421/1‑100 (MQ=60)
CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT > SRR3722112.712294/1‑100 (MQ=60)
CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT > SRR3722112.228958/1‑100 (MQ=60)
GGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAC < SRR3722112.317385/100‑1 (MQ=60)
GGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAC < SRR3722112.560329/100‑1 (MQ=60)
GGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAC > SRR3722112.173964/1‑100 (MQ=60)
GGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAC > SRR3722112.187230/1‑100 (MQ=60)
aagagacagGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACT < SRR3722112.597699/91‑1 (MQ=60)
GATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACT < SRR3722112.338605/100‑1 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722112.118347/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722112.317096/1‑100 (MQ=60)
GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT > SRR3722112.553164/1‑100 (MQ=60)
GGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCG > SRR3722112.34252/1‑100 (MQ=60)
GTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGA > SRR3722112.102355/1‑100 (MQ=60)
GTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGA < SRR3722112.113249/100‑1 (MQ=60)
GTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGA < SRR3722112.516453/100‑1 (MQ=60)
GTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGA > SRR3722112.316035/1‑100 (MQ=60)
AGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAG < SRR3722112.8302/100‑1 (MQ=60)
AAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGAT < SRR3722112.199729/100‑1 (MQ=60)
AAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGAT < SRR3722112.291369/100‑1 (MQ=60)
AAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGAT < SRR3722112.50574/100‑1 (MQ=60)
CGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGG > SRR3722112.40478/1‑100 (MQ=60)
CTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGG < SRR3722112.88137/100‑1 (MQ=60)
ATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAG > SRR3722112.213188/1‑100 (MQ=60)
ATCAAAAACCACGGTTAAAAAATGCCTCGACCGTCTGTTCCACCATTTCAAACTTGTTGTGCGGACTGAAGGCTATATCCCAGGATggtaaggaacgccg > SRR3722112.336516/1‑86 (MQ=50)
CCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAACGA < SRR3722112.211730/100‑1 (MQ=60)
|
TTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAGTGCAACGA > NZ_CP009273/2554029‑2554220
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |