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,661,508 |
T→G |
T250P (ACC→CCC) |
hcaR ← |
DNA‑binding transcriptional regulator HcaR |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,661,508 | 0 | T | G | 97.1%
| 103.4
/ ‑3.7
| 34 | T250P (ACC→CCC) | hcaR | DNA‑binding transcriptional regulator HcaR |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (15/18); total (16/18) |
Fisher's exact test for biased strand distribution p-value = 4.71e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
GGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT > NZ_CP009273/2661423‑2661590
|
ggCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTa > 2:822016/1‑90 (MQ=255)
gCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTaa < 1:455619/90‑1 (MQ=255)
ttCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTaaaa < 1:925026/90‑1 (MQ=255)
cATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatt > 1:948775/1‑90 (MQ=255)
aTCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatta < 2:987006/90‑1 (MQ=255)
tctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAattat < 1:1097578/90‑1 (MQ=255)
tttCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCAtata < 2:483946/90‑1 (MQ=255)
ttcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATaa > 1:542806/1‑90 (MQ=255)
tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc < 2:884117/90‑1 (MQ=255)
tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc < 1:640412/90‑1 (MQ=255)
ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg < 1:623295/90‑1 (MQ=255)
ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg < 2:846036/90‑1 (MQ=255)
ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg < 1:110462/90‑1 (MQ=255)
ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg < 2:345244/90‑1 (MQ=255)
ccACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCggg > 1:1093101/1‑90 (MQ=255)
ccATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCa > 2:521741/1‑90 (MQ=255)
aTCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCaaa > 1:1578168/1‑90 (MQ=255)
gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa < 2:674011/90‑1 (MQ=255)
gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa < 2:1130019/90‑1 (MQ=255)
aaaGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACg > 1:1503148/1‑90 (MQ=255)
aaaGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACg > 2:17161/1‑90 (MQ=255)
gCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACACCCTGTCCGGGATTAAAATTATTCATATAACTGGGTATCAAAGTGACGccc < 1:1479447/90‑1 (MQ=255)
cAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCa > 1:1460233/1‑90 (MQ=255)
aaTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAg > 2:1044315/1‑90 (MQ=255)
aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc > 2:1552280/1‑90 (MQ=255)
aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc > 2:1485974/1‑90 (MQ=255)
aGGCACATTTCCGGCAATGGGACGAAAAACAACCTTTCCGGGATTAAAATTATTCATAGAACCGGGTATCAAAGTGACGCCCACCCCCAt > 2:1479447/1‑90 (MQ=255)
aaTGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAt < 2:962133/90‑1 (MQ=255)
cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc < 2:1460233/90‑1 (MQ=255)
cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc < 2:839980/90‑1 (MQ=255)
gAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCa < 2:246119/90‑1 (MQ=255)
aaaaCAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGa > 1:935478/1‑90 (MQ=255)
tGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGtt > 1:1202669/1‑90 (MQ=255)
tGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGtt > 1:1502261/1‑90 (MQ=255)
|
GGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT > NZ_CP009273/2661423‑2661590
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
ATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC > NZ_CP009273/2661410‑2661601
|
ATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGA > SRR3722078.1537928/1‑100 (MQ=60)
TCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGAT < SRR3722078.299176/100‑1 (MQ=60)
TAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAA < SRR3722078.1425325/100‑1 (MQ=60)
ATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATT > SRR3722078.962072/1‑100 (MQ=60)
GCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCAT < SRR3722078.461063/100‑1 (MQ=60)
TTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATAT < SRR3722078.937715/100‑1 (MQ=60)
CATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAA > SRR3722078.549218/1‑100 (MQ=60)
TCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACC < SRR3722078.1114276/100‑1 (MQ=60)
CCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGG > SRR3722078.1109708/1‑100 (MQ=60)
TCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAA < SRR3722078.647781/100‑1 (MQ=60)
CTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAA > SRR3722078.1603854/1‑100 (MQ=60)
TTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAG < SRR3722078.111766/100‑1 (MQ=60)
TTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAG < SRR3722078.630477/100‑1 (MQ=60)
CGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACG > SRR3722078.1527666/1‑100 (MQ=60)
ATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCA > SRR3722078.1484050/1‑100 (MQ=60)
nnnnnnnnnnnnnnnnttaaGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGC < SRR3722078.1452875/80‑1 (MQ=60)
GCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACACCCTGTCCGGGATTAAAATTATTCATATAACTGGGTATCAAAGTGACGCCCAGCCCCATGC < SRR3722078.1503573/100‑1 (MQ=60)
tcgtcggcagcgtcagatgtgtataagagacaGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGG < SRR3722078.1522781/68‑1 (MQ=60)
AATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGA > SRR3722078.948422/1‑100 (MQ=60)
AAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT > SRR3722078.1526764/1‑100 (MQ=60)
AAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT > SRR3722078.1221531/1‑100 (MQ=60)
GTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC > SRR3722078.34327/1‑100 (MQ=60)
|
ATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGCCACCTGGAC > NZ_CP009273/2661410‑2661601
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |