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 I208 R1
|
222 |
13.1 |
731276 |
96.6% |
706412 |
85.9 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
4,613,222 |
C→G |
R238P (CGC→CCC) |
lplA ← |
lipoate‑‑protein ligase LplA |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 4,613,222 | 0 | C | G | 92.3%
| 39.7
/ ‑2.6
| 13 | R238P (CGC→CCC) | lplA | lipoate‑‑protein ligase LplA |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base G (11/1); total (11/2) |
Fisher's exact test for biased strand distribution p-value = 1.54e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01 |
CCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATG > NZ_CP009273/4613139‑4613304
|
ccACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCaa < 1:132271/90‑1 (MQ=255)
aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt > 1:286379/1‑90 (MQ=255)
aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt > 2:90161/1‑90 (MQ=255)
gTTCATCCAGCAGATGCGAGAATGCCGGAACCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa < 1:122617/90‑1 (MQ=255)
ttCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCaa > 2:143977/1‑90 (MQ=255)
tCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAg > 2:132709/1‑90 (MQ=255)
gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg > 2:316460/1‑90 (MQ=255)
gACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCg > 1:326502/1‑90 (MQ=255)
gAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCtta > 1:210220/1‑89 (MQ=255)
aTTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCt > 2:21011/1‑90 (MQ=255)
aTTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCt > 2:57078/1‑90 (MQ=255)
cAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCa > 2:321482/1‑90 (MQ=255)
aCTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATg > 2:158830/1‑90 (MQ=255)
|
CCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATG > NZ_CP009273/4613139‑4613304
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGA > NZ_CP009273/4613124‑4613309
|
ggctcggagatgtgtataagagacagCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGG < SRR3722090.318956/74‑1 (MQ=60)
CCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGC < SRR3722090.133671/100‑1 (MQ=60)
GCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTT > SRR3722090.289957/1‑100 (MQ=60)
GTTCATCCAGCAGATGCGAGAATGCCGGAACCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGT < SRR3722090.123921/100‑1 (MQ=60)
GCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCG > SRR3722090.330740/1‑100 (MQ=60)
GAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTA > SRR3722090.212558/1‑100 (MQ=60)
CAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATctgtctctt > SRR3722090.155987/1‑91 (MQ=60)
|
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGA > NZ_CP009273/4613124‑4613309
|
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: ‑ |