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 F3 I95 R1
|
196 |
14.2 |
808436 |
96.0% |
776098 |
85.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
3,898,828 |
(T)5→4 |
coding (130/1878 nt) |
bglF ← |
PTS beta‑glucoside transporter subunit IIABC |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 3,898,824 | 0 | T | . | 96.2%
| 97.9
/ ‑1.2
| 26 | coding (134/1878 nt) | bglF | PTS beta‑glucoside transporter subunit IIABC |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (14/11); total (15/11) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 |
TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC > NZ_CP009273/3898735‑3898912
|
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt > 1:346453/1‑90 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt < 1:156178/90‑1 (MQ=255)
ccAGGAAGACATCGGCCACATGGTTACCTATGACCGCCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag > 2:15339/1‑87 (MQ=255)
ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag > 2:349298/1‑87 (MQ=255)
ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag > 1:220387/1‑87 (MQ=255)
acaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCgctt < 1:232682/90‑1 (MQ=255)
tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGTGCTTTGCTTTCATCCtt > 2:223019/1‑90 (MQ=255)
accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTAcc < 1:105257/59‑1 (MQ=255)
accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTAcc > 2:105257/1‑59 (MQ=255)
accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt > 1:343584/1‑90 (MQ=255)
accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt > 2:209829/1‑90 (MQ=255)
ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt > 1:291516/1‑90 (MQ=255)
ccTGAAACTGGCCACCGCTTCCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt > 2:224721/1‑90 (MQ=255)
tGAAACTGGCCACCGCTTTCCACCCCCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatttaa < 2:137139/90‑1 (MQ=255)
ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa > 1:102889/1‑90 (MQ=255)
gCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt < 1:188034/90‑1 (MQ=255)
caccaccaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa < 1:349298/90‑1 (MQ=255)
ccaccaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg > 1:103273/1‑90 (MQ=255)
ccaccaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg > 1:40204/1‑90 (MQ=255)
ccaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt > 2:226081/1‑90 (MQ=255)
taataCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACt < 1:109246/90‑1 (MQ=255)
ataCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca > 1:328551/1‑90 (MQ=255)
ataCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca < 1:224721/90‑1 (MQ=255)
ataCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca < 1:209829/90‑1 (MQ=255)
aCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCGATGCATCAGACTcaca < 2:260606/90‑1 (MQ=255)
ccGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt > 1:315442/1‑67 (MQ=255)
ccGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt < 2:315442/67‑1 (MQ=255)
gggggTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa > 1:247366/1‑72 (MQ=255)
gggggTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa < 2:247366/72‑1 (MQ=255)
gggggTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg > 2:81623/1‑90 (MQ=255)
cTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATc > 1:216332/1‑90 (MQ=255)
|
TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC > NZ_CP009273/3898735‑3898912
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT > NZ_CP009273/3898731‑3898925
|
CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAG > SRR3722235.223159/1‑100 (MQ=60)
tctcgtgggctcggagatgtgtataagagacagCCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTAC < SRR3722235.384983/67‑1 (MQ=60)
TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACC < SRR3722235.157966/100‑1 (MQ=60)
ggagatgtgtataagagacagACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCG < SRR3722235.106488/79‑1 (MQ=60)
ACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATC < SRR3722235.235624/100‑1 (MQ=60)
ACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATT > SRR3722235.348518/1‑100 (MQ=60)
CCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTT > SRR3722235.295573/1‑100 (MQ=60)
ACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAA > SRR3722235.104093/1‑100 (MQ=60)
GCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGC < SRR3722235.190240/100‑1 (MQ=60)
CCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATG > SRR3722235.104483/1‑100 (MQ=60)
CCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATG > SRR3722235.40648/1‑100 (MQ=60)
CACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGAC < SRR3722235.354332/100‑1 (MQ=60)
CACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCA > SRR3722235.333234/1‑100 (MQ=60)
CCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCtgtctcttatacacatctccgagcccacgagacaagaggcaatctcgta > SRR3722235.171872/1‑51 (MQ=60)
CATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCctgtctcttatac > SRR3722235.319891/1‑87 (MQ=60)
TAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACctgtctct > SRR3722235.250541/1‑92 (MQ=60)
TAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTA < SRR3722235.110522/100‑1 (MQ=60)
ATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC < SRR3722235.227554/100‑1 (MQ=60)
ATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAACGTTATC < SRR3722235.212414/100‑1 (MQ=60)
TACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC > SRR3722235.219037/1‑100 (MQ=60)
gacagTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCC < SRR3722235.234989/95‑1 (MQ=60)
tcTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCG < SRR3722235.43091/98‑1 (MQ=60)
TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTctgtctcttatacacatctccgagcccacgagacaagaggcaa > SRR3722235.331593/1‑57 (MQ=60)
TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT > SRR3722235.267290/1‑100 (MQ=60)
TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT > SRR3722235.283768/1‑100 (MQ=60)
TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT > SRR3722235.365230/1‑100 (MQ=60)
TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT > SRR3722235.12571/1‑100 (MQ=60)
|
CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT > NZ_CP009273/3898731‑3898925
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
