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 I225 R1
|
227 |
21.4 |
1179702 |
97.1% |
1145490 |
86.5 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
2,012,918 |
A→G |
D351G (GAC→GGC) |
fliK → |
flagellar hook‑length control protein FliK |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,012,918 | 0 | A | G | 92.3%
| 29.7
/ ‑3.0
| 13 | D351G (GAC→GGC) | fliK | flagellar hook‑length control protein FliK |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (4/8); total (4/9) |
| 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.50e-01 |
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAG > NZ_CP009273/2012838‑2013002
|
cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATAc < 1:166371/90‑1 (MQ=255)
cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGtt > 1:555458/1‑90 (MQ=255)
cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc < 2:555458/90‑1 (MQ=255)
ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 1:212524/90‑1 (MQ=255)
ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 1:398853/90‑1 (MQ=255)
ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 1:514892/90‑1 (MQ=255)
ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc > 2:146167/1‑90 (MQ=255)
ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 2:355368/90‑1 (MQ=255)
aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAg > 2:529173/1‑90 (MQ=255)
aTGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt < 1:155160/90‑1 (MQ=255)
ctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa > 2:33490/1‑90 (MQ=255)
gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 2:36770/90‑1 (MQ=255)
ggAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag < 2:225704/90‑1 (MQ=255)
|
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAG > NZ_CP009273/2012838‑2013002
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ 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
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATC > NZ_CP009273/2012838‑2013016
|
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTT < SRR3722109.168430/100‑1 (MQ=60)
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTT > SRR3722109.563895/1‑100 (MQ=60)
GGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGG < SRR3722109.215080/100‑1 (MQ=60)
GGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGG < SRR3722109.404600/100‑1 (MQ=60)
GGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGG < SRR3722109.522713/100‑1 (MQ=60)
ATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACG < SRR3722109.157077/100‑1 (MQ=60)
GGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATC < SRR3722109.392410/100‑1 (MQ=60)
|
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATC > NZ_CP009273/2012838‑2013016
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
