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 I228 R1
|
214 |
26.5 |
1448192 |
97.1% |
1406194 |
86.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
3,901,121 |
A→G |
D202D (GAT→GAC) |
pstB ← |
phosphate ABC transporter ATP‑binding protein PstB |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 3,901,121 | 0 | A | G | 96.8%
| 96.4
/ ‑3.7
| 31 | D202D (GAT→GAC) | pstB | phosphate ABC transporter ATP‑binding protein PstB |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (17/13); total (17/14) |
| Fisher's exact test for biased strand distribution p-value = 4.52e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01 |
TCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGC > NZ_CP009273/3901041‑3901205
|
tCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTc < 1:424823/90‑1 (MQ=255)
gCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTcgg > 2:239860/1‑88 (MQ=255)
aaCGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGc < 2:367808/90‑1 (MQ=255)
cGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctc > 1:378373/1‑90 (MQ=255)
cGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctc > 2:116131/1‑90 (MQ=255)
cGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctc > 1:71588/1‑90 (MQ=255)
gCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGctct > 1:188663/1‑90 (MQ=255)
gtgtGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCa > 1:108313/1‑90 (MQ=255)
gtgGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAAt > 2:531315/1‑90 (MQ=255)
gtgGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAAt > 2:397387/1‑90 (MQ=255)
ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCGGTTCGGTGATCAGCTCTTCAATACGccc < 1:239860/90‑1 (MQ=255)
ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc < 2:36494/90‑1 (MQ=255)
ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc < 1:480938/90‑1 (MQ=255)
ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc < 1:720108/90‑1 (MQ=255)
ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc < 1:73290/90‑1 (MQ=255)
acGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCt < 2:85556/50‑1 (MQ=255)
acGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCt > 1:85556/1‑50 (MQ=255)
cAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAg > 2:498765/1‑90 (MQ=255)
ctgcATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAg < 2:309184/90‑1 (MQ=255)
aTGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGcc < 1:127511/90‑1 (MQ=255)
ttgtGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGaa > 1:695783/1‑90 (MQ=255)
ggTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACgg > 2:539680/1‑90 (MQ=255)
tGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGtt > 1:348776/1‑90 (MQ=255)
cGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtcgt < 1:271662/90‑1 (MQ=255)
tCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGa < 1:33324/90‑1 (MQ=255)
caccacGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc < 1:679146/51‑1 (MQ=255)
caccacGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc > 2:679146/1‑51 (MQ=255)
accacGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagc > 1:281024/1‑90 (MQ=255)
cGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagcagca > 2:499619/1‑90 (MQ=255)
gtgtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTTGAGCGCCGAACACGGTTCGTCGagcagcagc < 2:549037/90‑1 (MQ=255)
gtgtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagcagcagc > 2:123186/1‑90 (MQ=255)
|
TCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGC > NZ_CP009273/3901041‑3901205
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC > NZ_CP009273/3901023‑3901215
|
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGT < SRR3722112.713774/100‑1 (MQ=60)
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGT > SRR3722112.366874/1‑100 (MQ=60)
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC < SRR3722112.436254/100‑1 (MQ=60)
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC < SRR3722112.365105/100‑1 (MQ=60)
TCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGA < SRR3722112.430500/100‑1 (MQ=60)
GGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTC > SRR3722112.383284/1‑100 (MQ=60)
GGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTC > SRR3722112.72461/1‑100 (MQ=60)
GTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCT > SRR3722112.190976/1‑100 (MQ=60)
CATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCA > SRR3722112.109662/1‑100 (MQ=60)
agatctacacactgacattcgtcggcagcgtcagatgtgtataagagacaGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAG < SRR3722112.605854/50‑1 (MQ=60)
GGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCctgtctcttatacacatctgacgctgccga > SRR3722112.86617/1‑70 (MQ=60)
cggcagcgtcagatgtgtataagagacagGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA < SRR3722112.689764/71‑1 (MQ=60)
GGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCGGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA < SRR3722112.242709/100‑1 (MQ=60)
GGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA < SRR3722112.731354/100‑1 (MQ=60)
GGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA < SRR3722112.74182/100‑1 (MQ=60)
GGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA < SRR3722112.487826/100‑1 (MQ=60)
CTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAA > SRR3722112.706646/1‑100 (MQ=60)
ATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTT < SRR3722112.129092/100‑1 (MQ=60)
ATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTT > SRR3722112.353156/1‑100 (MQ=60)
GGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGC > SRR3722112.284287/1‑100 (MQ=60)
CGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAG < SRR3722112.274839/100‑1 (MQ=60)
TCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCAC < SRR3722112.33744/100‑1 (MQ=60)
GGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGG > SRR3722112.388483/1‑100 (MQ=60)
GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC > SRR3722112.4149/1‑100 (MQ=60)
GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC > SRR3722112.413643/1‑100 (MQ=60)
|
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC > NZ_CP009273/3901023‑3901215
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
