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 I225 R1
|
227 |
21.4 |
1179702 |
97.1% |
1145490 |
86.5 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
2,726,562 |
(T)5→4 |
coding (971/2574 nt) |
clpB ← |
ATP‑dependent chaperone ClpB |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,726,558 | 0 | T | . | 88.9%
| 50.3
/ 4.8
| 18 | coding (975/2574 nt) | clpB | ATP‑dependent chaperone ClpB |
Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (9/7); total (10/8) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA > NZ_CP009273/2726470‑2726640
|
cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt < 2:198266/90‑1 (MQ=255)
tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttc < 1:308475/90‑2 (MQ=255)
gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATg < 2:225916/90‑1 (MQ=255)
ccACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGt < 1:81688/90‑1 (MQ=255)
aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTAc > 2:526581/1‑90 (MQ=255)
gCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTg < 2:192247/90‑1 (MQ=255)
ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc > 2:541406/1‑90 (MQ=255)
ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCa > 2:50238/1‑90 (MQ=255)
ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCa > 2:256582/1‑90 (MQ=255)
aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg < 1:1227/90‑1 (MQ=255)
gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg > 2:551556/1‑86 (MQ=255)
gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg < 1:551556/86‑1 (MQ=255)
aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc < 2:29325/90‑1 (MQ=255)
aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt > 1:272977/1‑90 (MQ=255)
ggCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa < 1:489812/90‑1 (MQ=255)
cTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACg > 2:464277/1‑90 (MQ=255)
aacaaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGc < 2:211609/90‑1 (MQ=255)
acgacgTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgc > 2:210963/1‑90 (MQ=255)
ttCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGt > 2:193552/1‑90 (MQ=255)
gcagcaTC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa > 2:336422/1‑90 (MQ=255)
|
CTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA > NZ_CP009273/2726470‑2726640
|
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
TCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGC > NZ_CP009273/2726473‑2726661
|
TCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGG < SRR3722109.312513/100‑1 (MQ=60)
CCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATA < SRR3722109.82693/100‑1 (MQ=60)
CTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCT > SRR3722109.276388/1‑100 (MQ=60)
AACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGC < SRR3722109.1242/100‑1 (MQ=60)
AACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGC < SRR3722109.559940/100‑1 (MQ=60)
GGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATC‑TTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAA < SRR3722109.497196/100‑1 (MQ=60)
TCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGC > SRR3722109.471307/1‑100 (MQ=60)
|
TCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGC > NZ_CP009273/2726473‑2726661
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |