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 I228 R1
|
214 |
26.5 |
1448192 |
97.1% |
1406194 |
86.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
1,782,276 |
T→G |
V192G (GTG→GGG) |
ppsR → |
posphoenolpyruvate synthetase regulatory kinase/phosphorylase PpsR |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 1,782,276 | 0 | T | G | 90.9%
| 28.1
/ ‑3.3
| 11 | V192G (GTG→GGG) | ppsR | posphoenolpyruvate synthetase regulatory kinase/phosphorylase PpsR |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (8/2); total (8/3) |
| Fisher's exact test for biased strand distribution p-value = 2.73e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01 |
GTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGA > NZ_CP009273/1782205‑1782360
|
gTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCaa > 1:676222/1‑90 (MQ=255)
gTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCaa > 2:658133/1‑90 (MQ=255)
ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt > 1:102658/1‑90 (MQ=255)
ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt > 2:87076/1‑90 (MQ=255)
cgcgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGc > 1:648648/1‑90 (MQ=255)
ccTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTg > 1:657958/1‑77 (MQ=255)
ccTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTg < 2:657958/77‑1 (MQ=255)
tttATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGaa < 1:441404/90‑1 (MQ=255)
ttATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAAc < 1:114872/90‑1 (MQ=255)
tGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCt > 1:16033/1‑90 (MQ=255)
tCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa > 2:536520/1‑90 (MQ=255)
|
GTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGA > NZ_CP009273/1782205‑1782360
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ 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
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA > NZ_CP009273/1782179‑1782369
|
ccaccgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacaGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC < SRR3722112.671362/44‑1 (MQ=60)
CCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAA > SRR3722112.686790/1‑100 (MQ=60)
GTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCT > SRR3722112.103934/1‑100 (MQ=60)
GTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGC > SRR3722112.658755/1‑100 (MQ=60)
GCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCTG > SRR3722112.668227/1‑100 (MQ=60)
ACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCT > SRR3722112.16233/1‑100 (MQ=60)
TTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGG < SRR3722112.447514/100‑1 (MQ=60)
TTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGC < SRR3722112.116298/100‑1 (MQ=60)
ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA > SRR3722112.245928/1‑100 (MQ=60)
|
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA > NZ_CP009273/1782179‑1782369
|
| 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: ‑ |
