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,818,451 |
T→C |
E162G (GAA→GGA) |
mltB ← |
lytic murein transglycosylase B |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,818,451 | 0 | T | C | 90.9%
| 52.9
/ ‑0.9
| 22 | E162G (GAA→GGA) | mltB | lytic murein transglycosylase B |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (2/18); total (3/19) |
| Fisher's exact test for biased strand distribution p-value = 2.60e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 |
CCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAA > NZ_CP009273/2818368‑2818540
|
ccGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGc > 2:260350/1‑90 (MQ=255)
ggTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCcga < 1:690408/90‑1 (MQ=255)
tAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATaa < 1:807536/90‑1 (MQ=255)
cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg < 1:572360/90‑1 (MQ=255)
cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg < 2:124730/90‑1 (MQ=255)
cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg < 1:102597/90‑1 (MQ=255)
cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg < 1:305978/90‑1 (MQ=255)
ttGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTa < 1:458397/90‑1 (MQ=255)
ttGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTa < 1:834426/90‑1 (MQ=255)
tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc < 2:220487/90‑1 (MQ=255)
tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc < 1:231874/90‑1 (MQ=255)
cAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTcc < 2:585335/90‑1 (MQ=255)
gcgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCACCGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAt < 1:498725/90‑1 (MQ=255)
gcgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAt < 2:131755/90‑1 (MQ=255)
gcgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAt < 1:260350/90‑1 (MQ=255)
cgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATa > 1:456/1‑90 (MQ=255)
aTGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATac < 2:788614/81‑1 (MQ=255)
aTGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATac > 1:788614/1‑81 (MQ=255)
gAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCAcgcgc < 1:771357/90‑1 (MQ=255)
gAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCAcgcgc < 1:412535/90‑1 (MQ=255)
cATCCCGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGc < 1:626283/90‑1 (MQ=255)
ccccAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATAc < 1:788915/90‑1 (MQ=255)
tCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGaaa < 2:60290/90‑1 (MQ=255)
|
CCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAA > NZ_CP009273/2818368‑2818540
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ 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
GGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGA > NZ_CP009273/2818380‑2818538
|
GGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTC < SRR3722079.697664/100‑1 (MQ=60)
TAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGG < SRR3722079.816198/100‑1 (MQ=60)
GCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATA > SRR3722079.461/1‑100 (MQ=60)
CGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATAC < SRR3722079.103420/100‑1 (MQ=60)
CGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATAC < SRR3722079.308467/100‑1 (MQ=60)
CGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATAC < SRR3722079.577990/100‑1 (MQ=60)
TTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACAC < SRR3722079.462407/100‑1 (MQ=60)
TTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACAC < SRR3722079.843359/100‑1 (MQ=60)
TGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACC < SRR3722079.233841/100‑1 (MQ=60)
GCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCACCGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCA < SRR3722079.503185/100‑1 (MQ=60)
GCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCA < SRR3722079.262535/100‑1 (MQ=60)
CGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCAC > SRR3722079.797064/1‑100 (MQ=60)
GAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGC < SRR3722079.416066/100‑1 (MQ=60)
GAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGC < SRR3722079.779617/100‑1 (MQ=60)
CATCCCGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATAC < SRR3722079.632672/100‑1 (MQ=60)
CCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGA < SRR3722079.797367/100‑1 (MQ=60)
|
GGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGA > NZ_CP009273/2818380‑2818538
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
