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 |
2,526,050 |
A→G |
T95A (ACC→GCC) |
cysK → |
cysteine synthase A |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,526,050 | 0 | A | G | 94.4%
| 103.0
/ ‑1.3
| 36 | T95A (ACC→GCC) | cysK | cysteine synthase A |
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (22/12); total (24/12) |
Fisher's exact test for biased strand distribution p-value = 5.43e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01 |
ACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA > NZ_CP009273/2525972‑2526130
|
aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTAtt > 2:335278/1‑90 (MQ=255)
aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACACTGACCATGCCAGAAGCCATGAGTAtt > 2:70613/1‑90 (MQ=255)
cGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct > 1:795726/1‑90 (MQ=255)
gATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctgct > 2:456013/1‑90 (MQ=255)
ttGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCTAGCTGCTGa < 1:362055/90‑1 (MQ=255)
tGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGaa > 1:671497/1‑90 (MQ=255)
gCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGaaa > 1:802451/1‑90 (MQ=255)
ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt > 2:320948/1‑90 (MQ=255)
ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt > 2:167688/1‑90 (MQ=255)
aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg > 2:324976/1‑90 (MQ=255)
aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg > 2:231849/1‑90 (MQ=255)
aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg > 1:831646/1‑90 (MQ=255)
aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg > 2:473230/1‑90 (MQ=255)
aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg > 2:759608/1‑90 (MQ=255)
gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc > 1:826533/1‑90 (MQ=255)
gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc > 1:554697/1‑90 (MQ=255)
gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc > 2:678039/1‑90 (MQ=255)
gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc > 2:319876/1‑90 (MQ=255)
tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCt > 2:605305/1‑90 (MQ=255)
gCCGCTCGCGGTTACAAACTCACCATGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTg > 2:585574/1‑90 (MQ=255)
cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt > 1:693970/1‑90 (MQ=255)
gTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACgg < 1:335869/90‑1 (MQ=255)
tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa < 2:54514/90‑1 (MQ=255)
tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa < 1:571678/90‑1 (MQ=255)
tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa < 1:456013/90‑1 (MQ=255)
tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa < 1:324976/90‑1 (MQ=255)
aCAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAg < 1:444169/90‑1 (MQ=255)
cAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGCACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAgg < 1:813746/90‑1 (MQ=255)
cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa < 1:759608/90‑1 (MQ=255)
cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa > 2:702301/1‑90 (MQ=255)
cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa < 2:795726/90‑1 (MQ=255)
gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGc < 2:353336/57‑1 (MQ=255)
gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGc > 1:353336/1‑57 (MQ=255)
gACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt > 1:63226/1‑90 (MQ=255)
aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa > 2:732679/1‑90 (MQ=255)
gAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa < 2:731777/42‑1 (MQ=255)
|
ACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA > NZ_CP009273/2525972‑2526130
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
GCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAGC > NZ_CP009273/2525976‑2526147
|
GCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCT > SRR3722079.804262/1‑100 (MQ=60)
ATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAA > SRR3722079.678518/1‑100 (MQ=60)
TACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAA > SRR3722079.811060/1‑100 (MQ=60)
TTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCTAGCTGCTGAAAGCGTTAGG < SRR3722079.364885/100‑1 (MQ=60)
GCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTG > SRR3722079.840549/1‑100 (MQ=60)
GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAC > SRR3722079.835381/1‑100 (MQ=60)
GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAC > SRR3722079.560037/1‑100 (MQ=60)
ATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGT > SRR3722079.701275/1‑100 (MQ=60)
GTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAA < SRR3722079.338527/100‑1 (MQ=60)
agatgtgtataagagacagTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAG < SRR3722079.323318/81‑1 (MQ=60)
TACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAG < SRR3722079.577294/100‑1 (MQ=60)
TACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAG < SRR3722079.460004/100‑1 (MQ=60)
TACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAG < SRR3722079.327584/100‑1 (MQ=60)
ACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGG < SRR3722079.448051/100‑1 (MQ=60)
CAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGCACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGC < SRR3722079.822469/100‑1 (MQ=60)
AACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGtctcttatacacatctccga > SRR3722079.356111/1‑80 (MQ=60)
AACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAT > SRR3722079.63743/1‑100 (MQ=60)
CACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA < SRR3722079.767732/100‑1 (MQ=60)
GACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCctgtctcttatacacatctccgagcccacgagaccaga > SRR3722079.739608/1‑62 (MQ=60)
GACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCA > SRR3722079.754450/1‑100 (MQ=60)
TGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCctgtct > SRR3722079.787811/1‑94 (MQ=60)
TGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCA > SRR3722079.755027/1‑100 (MQ=60)
TGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCA > SRR3722079.204350/1‑100 (MQ=60)
GAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAG > SRR3722079.489556/1‑100 (MQ=60)
AAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAGC < SRR3722079.575473/100‑1 (MQ=60)
|
GCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAGC > NZ_CP009273/2525976‑2526147
|
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: ‑ |