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 I227 R1
|
223 |
21.3 |
1167702 |
97.0% |
1132670 |
86.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
2,189,312 |
A→G |
T512A (ACC→GCC) |
metG → |
methionine‑‑tRNA ligase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,189,312 | 0 | A | G | 94.7%
| 52.1
/ ‑3.9
| 19 | T512A (ACC→GCC) | metG | methionine‑‑tRNA ligase |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (10/8); total (11/8) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01 |
ACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAG > NZ_CP009273/2189226‑2189395
|
aCCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCt < 1:440017/90‑1 (MQ=255)
gcgTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTa > 1:262809/1‑90 (MQ=255)
gcgTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTa > 1:449721/1‑90 (MQ=255)
gTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATc > 1:204983/1‑90 (MQ=255)
tACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCgctg > 1:582320/1‑90 (MQ=255)
ggTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACaa > 1:335601/1‑90 (MQ=255)
gCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGaa > 1:126917/1‑90 (MQ=255)
cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc < 1:337883/90‑1 (MQ=255)
cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc < 2:204983/90‑1 (MQ=255)
aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc < 2:557649/90‑1 (MQ=255)
aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc < 2:262809/90‑1 (MQ=255)
gTGCAGGAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGt > 2:367428/1‑90 (MQ=255)
aCGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATc > 1:522283/1‑55 (MQ=255)
aCGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATc < 2:522283/55‑1 (MQ=255)
gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc < 1:529739/62‑1 (MQ=255)
gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc > 2:529739/1‑62 (MQ=255)
gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGc < 2:201773/88‑1 (MQ=255)
gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGc > 1:201773/1‑88 (MQ=255)
gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAg > 1:392138/1‑90 (MQ=255)
|
ACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAG > NZ_CP009273/2189226‑2189395
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
AACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGAAGCAC > NZ_CP009273/2189225‑2189405
|
AACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTA > SRR3722111.265923/1‑100 (MQ=60)
AACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTA > SRR3722111.456537/1‑100 (MQ=60)
ACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTAT < SRR3722111.446620/100‑1 (MQ=60)
CCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATC > SRR3722111.207512/1‑100 (MQ=60)
GCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTG > SRR3722111.591451/1‑100 (MQ=60)
ACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAA > SRR3722111.340088/1‑100 (MQ=60)
AGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAA > SRR3722111.128509/1‑100 (MQ=60)
gatgtgtataagagacagCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAAC < SRR3722111.538054/82‑1 (MQ=60)
CGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAAC < SRR3722111.342407/100‑1 (MQ=60)
ATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCctgtctcttatacacatctgacgct > SRR3722111.530450/1‑75 (MQ=60)
CCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAG > SRR3722111.204275/1‑100 (MQ=60)
CCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAG > SRR3722111.397821/1‑100 (MQ=60)
ATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCctgtctc > SRR3722111.248563/1‑93 (MQ=60)
ATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCctgtctc > SRR3722111.531553/1‑93 (MQ=60)
ATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGA > SRR3722111.78081/1‑100 (MQ=60)
ATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGA > SRR3722111.165836/1‑100 (MQ=60)
GAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCctgtctcttata > SRR3722111.145089/1‑88 (MQ=60)
|
AACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGAAGCAC > NZ_CP009273/2189225‑2189405
|
| 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: ‑ |