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 |
354,025 |
T→C |
E126G (GAA→GGA) |
cynR ← |
transcriptional regulator CynR |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 354,025 | 0 | T | C | 77.3%
| 32.5
/ 8.8
| 22 | E126G (GAA→GGA) | cynR | transcriptional regulator CynR |
Reads supporting (aligned to +/- strand): ref base T (0/5); new base C (5/12); total (5/17) |
Fisher's exact test for biased strand distribution p-value = 2.90e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01 |
TCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAA > NZ_CP009273/353941‑354100
|
tCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGt < 2:202025/90‑1 (MQ=255)
gAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTgg > 1:680999/1‑90 (MQ=255)
cacaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGt < 2:510876/90‑1 (MQ=255)
acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTg < 2:471502/90‑1 (MQ=255)
acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg < 2:636746/90‑1 (MQ=255)
acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg < 1:57534/90‑1 (MQ=255)
acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg < 1:702142/90‑1 (MQ=255)
gcgcGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGc < 1:370805/90‑1 (MQ=255)
cAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATc < 1:809092/90‑1 (MQ=255)
tcgtcgCGGCAAAGCATATCCTGGATTTTCCCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGcc < 2:668259/90‑1 (MQ=255)
tcgtcgCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGcc < 1:383976/90‑1 (MQ=255)
cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg < 2:297407/90‑1 (MQ=255)
cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg < 2:418418/90‑1 (MQ=255)
cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg < 1:173485/90‑1 (MQ=255)
tCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAg > 2:291368/1‑90 (MQ=255)
tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg > 1:487236/1‑71 (MQ=255)
tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg < 2:487236/71‑1 (MQ=255)
tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa > 1:639768/1‑90 (MQ=255)
tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa < 1:469055/90‑1 (MQ=255)
tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa < 1:290942/90‑1 (MQ=255)
ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAg > 2:226559/1‑67 (MQ=255)
ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAg < 1:226559/67‑1 (MQ=255)
|
TCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAA > NZ_CP009273/353941‑354100
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
GGAATTGCCTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGCGAT > NZ_CP009273/353926‑354119
|
GGAATTGCCTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTC < SRR3722079.622498/100‑1 (MQ=60)
CTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGC < SRR3722079.654476/100‑1 (MQ=60)
CAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGG > SRR3722079.688150/1‑100 (MQ=60)
ACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGAT < SRR3722079.709583/100‑1 (MQ=60)
ACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGAT < SRR3722079.57992/100‑1 (MQ=60)
GCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCG < SRR3722079.373781/100‑1 (MQ=60)
CAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAG < SRR3722079.817774/100‑1 (MQ=60)
TCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGC < SRR3722079.387123/100‑1 (MQ=60)
gtataagagacagCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAA < SRR3722079.228485/87‑1 (MQ=60)
CAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGT < SRR3722079.174994/100‑1 (MQ=60)
TCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTctgtctctt > SRR3722079.491507/1‑91 (MQ=60)
TCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAA > SRR3722079.646357/1‑100 (MQ=60)
TCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGT < SRR3722079.473156/100‑1 (MQ=60)
TCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGT < SRR3722079.293340/100‑1 (MQ=60)
TCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGA > SRR3722079.124946/1‑100 (MQ=60)
GCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGC > SRR3722079.412082/1‑100 (MQ=60)
ACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGCGAT > SRR3722079.547379/1‑100 (MQ=60)
|
GGAATTGCCTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGCGAT > NZ_CP009273/353926‑354119
|
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: ‑ |