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 I199 R1
|
220 |
57.5 |
3163680 |
97.0% |
3068769 |
86.1 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
2,751,237 |
A→G |
V241V (GTT→GTC) |
BW25113_RS13685 ← |
hypothetical protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,751,237 | 0 | A | G | 92.5%
| 108.4
/ 2.2
| 40 | V241V (GTT→GTC) | BW25113_RS13685 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (17/20); total (18/22) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01 |
AACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCAC > NZ_CP009273/2751154‑2751322
|
aaCTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACagag > 1:348772/1‑90 (MQ=255)
aaCTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACagag > 1:840964/1‑90 (MQ=255)
aTTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCattaatt < 2:876114/90‑1 (MQ=255)
tGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGc < 2:700612/90‑1 (MQ=255)
ccATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACtt > 1:1013126/1‑90 (MQ=255)
cATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTc < 2:621469/90‑1 (MQ=255)
aGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGt < 2:296209/90‑1 (MQ=255)
ggTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGa > 2:320069/1‑90 (MQ=255)
ggTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGa > 2:804685/1‑90 (MQ=255)
gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGaa > 1:960669/1‑90 (MQ=255)
gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGaa > 2:743414/1‑90 (MQ=255)
gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGaa > 2:764374/1‑90 (MQ=255)
ttCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAAt < 1:570471/90‑1 (MQ=255)
ttCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAAt < 1:1353522/90‑1 (MQ=255)
tCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATg > 2:707544/1‑90 (MQ=255)
tAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGatat < 2:333611/90‑1 (MQ=255)
tAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGatat < 1:1552570/90‑1 (MQ=255)
tAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGatat < 2:293138/90‑1 (MQ=255)
tCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTac < 2:1559963/90‑1 (MQ=255)
tCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTac < 2:244530/90‑1 (MQ=255)
tGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgc > 2:698890/1‑90 (MQ=255)
tGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgc > 1:1414339/1‑90 (MQ=255)
tGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgc > 2:48107/1‑90 (MQ=255)
tGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgc > 2:1279610/1‑90 (MQ=255)
cTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctt < 2:960669/90‑1 (MQ=255)
cTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctt < 1:1243839/90‑1 (MQ=255)
cTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctt < 2:840964/90‑1 (MQ=255)
aaCCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctttgct < 1:246862/90‑1 (MQ=255)
aaCCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctttgct < 2:844344/90‑1 (MQ=255)
ccTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctttgcttt > 1:975903/1‑90 (MQ=255)
tATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGc > 2:469201/1‑90 (MQ=255)
tATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGc > 2:550485/1‑90 (MQ=255)
ttttCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCAc < 1:804685/90‑1 (MQ=255)
ttttCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCAc < 1:698890/90‑1 (MQ=255)
ttttCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCAc < 1:1113393/90‑1 (MQ=255)
tCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGatat < 2:1048663/60‑1 (MQ=255)
tCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGatat > 1:1048663/1‑60 (MQ=255)
tCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCct > 1:1432931/1‑90 (MQ=255)
tCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCAc < 2:1414339/90‑1 (MQ=255)
tCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCAc < 1:839598/90‑1 (MQ=255)
|
AACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCAC > NZ_CP009273/2751154‑2751322
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
TTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTC > NZ_CP009273/2751138‑2751335
|
TTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGG < SRR3722078.599737/100‑1 (MQ=60)
TTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGG < SRR3722078.388404/100‑1 (MQ=60)
GATTATTTATCTCAAACTCTGCCATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAGCCTTATTTTTCCACCCATCATTCGGGAC < SRR3722078.1096517/100‑1 (MQ=60)
GATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGAC < SRR3722078.1063022/100‑1 (MQ=60)
GATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGAC < SRR3722078.13050/100‑1 (MQ=60)
ATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAG > SRR3722078.353045/1‑100 (MQ=60)
ATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAG > SRR3722078.852173/1‑100 (MQ=60)
GAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTT > SRR3722078.1028050/1‑100 (MQ=60)
TCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAA > SRR3722078.974304/1‑100 (MQ=60)
TTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTAC < SRR3722078.577151/100‑1 (MQ=60)
TTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTAC < SRR3722078.1375247/100‑1 (MQ=60)
acgaccgccatcgtcggcagcgtcagatgtgtataagagacagATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACAC < SRR3722078.823166/57‑1 (MQ=60)
TAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTC < SRR3722078.1577853/100‑1 (MQ=60)
GTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGC > SRR3722078.1437305/1‑100 (MQ=60)
ACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTT > SRR3722078.989924/1‑100 (MQ=60)
GCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTctgtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatctcggtg > SRR3722078.1020836/1‑44 (MQ=60)
CTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCC < SRR3722078.1263503/100‑1 (MQ=60)
AACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCC < SRR3722078.249876/100‑1 (MQ=60)
ACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCtgtctcttatacacatctg > SRR3722078.1064375/1‑81 (MQ=60)
ACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCT > SRR3722078.1456233/1‑100 (MQ=60)
TTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCC < SRR3722078.707293/100‑1 (MQ=60)
TTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCC < SRR3722078.815297/100‑1 (MQ=60)
TTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCC < SRR3722078.1130391/100‑1 (MQ=60)
ATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTG > SRR3722078.1172770/1‑100 (MQ=60)
TCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGAC < SRR3722078.850776/100‑1 (MQ=60)
GGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACctgtctcttatac > SRR3722078.1491015/1‑87 (MQ=60)
GGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTC < SRR3722078.752829/100‑1 (MQ=60)
GGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTC < SRR3722078.1026869/100‑1 (MQ=60)
|
TTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTC > NZ_CP009273/2751138‑2751335
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |