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 I207 R1
|
218 |
45.0 |
2482936 |
97.1% |
2410930 |
86.1 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
4,180,586 |
(A)7→6 |
intergenic (+333/+77) |
yjaZ → / ← thiH |
heat‑shock protein C/2‑iminoacetate synthase ThiH |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 4,180,580 | 0 | A | . | 96.5%
| 114.1
/ ‑1.1
| 29 | intergenic (+327/+83) | yjaZ/thiH | heat‑shock protein C/2‑iminoacetate synthase ThiH |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (9/19); total (9/20) |
| 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.98e-01 |
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG > NZ_CP009273/4180492‑4180660
|
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGC‑TTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaa > 1:1188275/1‑90 (MQ=255)
tGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaa > 1:635534/1‑90 (MQ=255)
ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACaa < 1:1129842/90‑1 (MQ=255)
aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGc < 2:276185/90‑1 (MQ=255)
aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGc < 2:579509/90‑1 (MQ=255)
gAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa < 1:317979/72‑1 (MQ=255)
gAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa > 2:317979/1‑72 (MQ=255)
gaagCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGt > 2:1211839/1‑90 (MQ=255)
gCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc > 2:80846/1‑90 (MQ=255)
tttCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGc > 1:634240/1‑90 (MQ=255)
tttCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGc > 1:1230623/1‑90 (MQ=255)
tCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc < 1:60835/90‑1 (MQ=255)
ccTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGaa < 2:941115/90‑1 (MQ=255)
cTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAAt < 2:734901/90‑1 (MQ=255)
aacaTTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAAttt < 2:195341/90‑1 (MQ=255)
acaTTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTa > 2:33928/1‑90 (MQ=255)
caTTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc < 1:540062/90‑1 (MQ=255)
caTTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc < 2:348119/90‑1 (MQ=255)
caTTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc < 1:80846/90‑1 (MQ=255)
ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt < 2:1177268/90‑1 (MQ=255)
ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt < 1:570762/90‑1 (MQ=255)
ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt < 1:546502/90‑1 (MQ=255)
ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt < 1:373789/90‑1 (MQ=255)
ccTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc > 2:567679/1‑90 (MQ=255)
ctgctgACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCt < 2:330468/90‑1 (MQ=255)
tgACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa > 1:458287/1‑90 (MQ=255)
gACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt > 2:1021463/1‑90 (MQ=255)
gACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt < 1:702169/90‑1 (MQ=255)
gACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt > 2:988008/1‑90 (MQ=255)
cGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc < 2:271223/90‑1 (MQ=255)
cGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc < 1:328156/90‑1 (MQ=255)
ggCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc < 2:792262/90‑1 (MQ=255)
tcACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATTTTCCGTAGCGGCTCTTAATAAATTTACGTTCTTACATTCCg < 2:316822/89‑1 (MQ=255)
|
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG > NZ_CP009273/4180492‑4180660
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAG > NZ_CP009273/4180494‑4180678
|
agagacaGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGC < SRR3722089.322007/93‑1 (MQ=60)
GGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAAC < SRR3722089.1147012/100‑1 (MQ=60)
nnnnnnnnnnnnnnnnnnnnnAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAAC < SRR3722089.110540/79‑1 (MQ=60)
GGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGC > SRR3722089.1249330/1‑100 (MQ=60)
GGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGC > SRR3722089.642316/1‑100 (MQ=60)
TCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT < SRR3722089.61577/100‑1 (MQ=60)
ACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACA > SRR3722089.463767/1‑100 (MQ=60)
CATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAT < SRR3722089.81831/100‑1 (MQ=60)
CATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAT < SRR3722089.546352/100‑1 (MQ=60)
TTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC < SRR3722089.552856/100‑1 (MQ=60)
TTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC < SRR3722089.577512/100‑1 (MQ=60)
TTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC < SRR3722089.378372/100‑1 (MQ=60)
GACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCAT < SRR3722089.711321/100‑1 (MQ=60)
CGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAG < SRR3722089.332275/100‑1 (MQ=60)
GCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC > SRR3722089.638664/1‑100 (MQ=60)
ACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATctgtctcttttacacatctgacgctgccgacgaccat > SRR3722089.1241903/1‑63 (MQ=60)
ACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGC > SRR3722089.948447/1‑100 (MQ=60)
atAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAG > SRR3722089.764776/3‑100 (MQ=60)
|
GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAG > NZ_CP009273/4180494‑4180678
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
