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 |
1,844,118 |
A→G |
E342G (GAA→GGA) |
sppA → |
signal peptide peptidase SppA |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,844,118 | 0 | A | G | 94.7%
| 51.4
/ ‑3.4
| 19 | E342G (GAA→GGA) | sppA | signal peptide peptidase SppA |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (4/14); total (4/15) |
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 |
ATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCT > NZ_CP009273/1844044‑1844202
|
aTTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGt < 1:134468/90‑1 (MQ=255)
aaaCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCg < 2:87375/90‑1 (MQ=255)
aaaCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCg < 2:201858/90‑1 (MQ=255)
aaaCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCg < 1:245475/90‑1 (MQ=255)
gCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGa > 2:110303/1‑90 (MQ=255)
gCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGa > 2:309698/1‑90 (MQ=255)
gCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATAcc < 1:158858/90‑1 (MQ=255)
gCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATAcc < 1:93051/90‑1 (MQ=255)
tACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGc < 2:309814/90‑1 (MQ=255)
gTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGcac < 2:416837/90‑1 (MQ=255)
aTTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGaaa < 1:521776/90‑1 (MQ=255)
tggatggCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGa < 2:414951/90‑1 (MQ=255)
ggatggCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGAt < 1:551262/90‑1 (MQ=255)
gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCtt < 2:105839/65‑1 (MQ=255)
gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCtt > 1:105839/1‑65 (MQ=255)
gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTcc > 2:473922/1‑90 (MQ=255)
gCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTcc < 2:86464/90‑1 (MQ=255)
cGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCt < 1:90953/90‑1 (MQ=255)
cGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCt < 2:73689/90‑1 (MQ=255)
|
ATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCT > NZ_CP009273/1844044‑1844202
|
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
CGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAAT > NZ_CP009273/1844020‑1844212
|
CGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGA > SRR3722111.87524/1‑100 (MQ=60)
AGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGG < SRR3722111.540535/100‑1 (MQ=60)
ATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGAT < SRR3722111.136167/100‑1 (MQ=60)
AAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCAC < SRR3722111.248378/100‑1 (MQ=60)
GCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCAC < SRR3722111.160868/100‑1 (MQ=60)
GCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCAC < SRR3722111.94197/100‑1 (MQ=60)
ATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGctgtctcttatacac > SRR3722111.107158/1‑85 (MQ=60)
ATTATGGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCC < SRR3722111.529930/100‑1 (MQ=60)
GGATGGCGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGT < SRR3722111.559921/100‑1 (MQ=60)
CGAGGGAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAAT < SRR3722111.92079/100‑1 (MQ=60)
|
CGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAAT > NZ_CP009273/1844020‑1844212
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: ‑ |