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 |
1,903,061 |
A→T |
intergenic (‑38/‑121) |
mgrB ← / → yobH |
PhoP/PhoQ regulator MgrB/YobH family protein |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 1,903,061 | 0 | A | T | 88.9%
| 84.7
/ 0.5
| 27 | intergenic (‑38/‑121) | mgrB/yobH | PhoP/PhoQ regulator MgrB/YobH family protein |
| Reads supporting (aligned to +/- strand): ref base A (0/3); new base T (11/13); total (11/16) |
| Fisher's exact test for biased strand distribution p-value = 2.48e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01 |
AGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTA > NZ_CP009273/1902975‑1903147
|
agcaagcCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCa < 2:840105/90‑1 (MQ=255)
gcaagcCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCAt < 2:681844/90‑1 (MQ=255)
ccAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCgtt < 2:205243/90‑1 (MQ=255)
ccacGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTaa < 2:18815/90‑1 (MQ=255)
acGACACACAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTCCACCTACCTTATGTCACTAAACTAGCATGCGTTAGTTAAAc < 2:536781/90‑1 (MQ=255)
acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc < 1:229040/90‑1 (MQ=255)
acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc > 1:454872/1‑90 (MQ=255)
cGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACa < 2:400589/90‑1 (MQ=255)
gACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAg < 2:781070/90‑1 (MQ=255)
aCCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGcc > 1:159755/1‑90 (MQ=255)
aCGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTa < 2:454872/90‑1 (MQ=255)
aCGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTa < 2:586/90‑1 (MQ=255)
ccATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCg > 2:435980/1‑90 (MQ=255)
cTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc > 2:561475/1‑64 (MQ=255)
cTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc < 1:561475/64‑1 (MQ=255)
tttCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCtt > 2:157475/1‑90 (MQ=255)
tcCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGt < 2:303445/90‑1 (MQ=255)
ccGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTc > 1:22023/1‑90 (MQ=255)
cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg > 2:16731/1‑90 (MQ=255)
ttCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgctt < 2:173236/90‑1 (MQ=255)
tacctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCa > 1:678742/1‑90 (MQ=255)
aTGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCat > 1:147528/1‑90 (MQ=255)
tCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATAtt < 1:763045/90‑1 (MQ=255)
aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc < 2:452589/90‑1 (MQ=255)
aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc > 2:77152/1‑90 (MQ=255)
ttAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTTTGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGcc > 1:141159/1‑90 (MQ=255)
aCTTGCATGTGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTa < 2:195223/90‑1 (MQ=255)
|
AGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTA > NZ_CP009273/1902975‑1903147
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ 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
AAGCAGCAAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGA > NZ_CP009273/1902967‑1903150
|
AAGCAGCAAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATG < SRR3722079.73182/100‑1 (MQ=60)
AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAC > SRR3722079.458853/1‑100 (MQ=60)
CACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCC > SRR3722079.161151/1‑100 (MQ=60)
tgtgtataagagacagCCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAAC < SRR3722079.566928/84‑1 (MQ=60)
ACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAAC < SRR3722079.230983/100‑1 (MQ=60)
ACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTC > SRR3722079.22175/1‑100 (MQ=60)
CCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCA > SRR3722079.685865/1‑100 (MQ=60)
CACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCAT > SRR3722079.148801/1‑100 (MQ=60)
CCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTTTGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCC > SRR3722079.142360/1‑100 (MQ=60)
TCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGA < SRR3722079.771207/100‑1 (MQ=60)
|
AAGCAGCAAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGA > NZ_CP009273/1902967‑1903150
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
