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 |
2,814,267 |
A→G |
A368A (GCT→GCC) |
alaS ← |
alanine‑‑tRNA ligase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,814,267 | 0 | A | G | 84.0%
| 51.1
/ 5.9
| 25 | A368A (GCT→GCC) | alaS | alanine‑‑tRNA ligase |
| Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (10/11); total (12/13) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-02 |
CGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGT > NZ_CP009273/2814181‑2814352
|
cGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCa > 1:401062/1‑90 (MQ=255)
ggAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaa < 1:372211/90‑1 (MQ=255)
aGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctc < 1:14428/90‑1 (MQ=255)
cATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcag < 1:3308/90‑1 (MQ=255)
cATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcag < 1:382642/90‑1 (MQ=255)
aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt > 2:57314/1‑90 (MQ=255)
ccAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtct > 1:453002/1‑90 (MQ=255)
gCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtcttca > 2:58525/1‑90 (MQ=255)
gCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtcttca > 1:700076/1‑90 (MQ=255)
tATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCa < 1:623815/90‑1 (MQ=255)
aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc < 1:58525/90‑1 (MQ=255)
tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAcc > 2:347389/1‑90 (MQ=255)
aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc > 1:103417/1‑90 (MQ=255)
ctTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGc < 2:654739/90‑1 (MQ=255)
cAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAg > 1:321303/1‑90 (MQ=255)
aaCGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCt < 2:103417/90‑1 (MQ=255)
aaCGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCt < 2:423389/90‑1 (MQ=255)
gCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCa > 1:362731/1‑90 (MQ=255)
cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag < 2:378903/90‑1 (MQ=255)
cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag < 2:453002/90‑1 (MQ=255)
cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag < 2:641048/90‑1 (MQ=255)
gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga > 2:581081/1‑90 (MQ=255)
gcgcTCCAGAGTACGAGCAAACTGCTCTtcttcagt > 1:534528/1‑36 (MQ=255)
gcgcTCCAGAGTACGAGCAAACTGCTCTtcttcagt < 2:534528/36‑1 (MQ=255)
tACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGt > 1:152216/1‑90 (MQ=255)
|
CGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGT > NZ_CP009273/2814181‑2814352
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
GTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGAT > NZ_CP009273/2814171‑2814356
|
GTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCA > SRR3722112.406339/1‑100 (MQ=60)
GGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTC < SRR3722112.377028/100‑1 (MQ=60)
AGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTC < SRR3722112.14615/100‑1 (MQ=60)
GTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCT > SRR3722112.459365/1‑100 (MQ=60)
TCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCA > SRR3722112.711016/1‑100 (MQ=60)
CATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAC < SRR3722112.3346/100‑1 (MQ=60)
CATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAC < SRR3722112.387615/100‑1 (MQ=60)
GCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTC > SRR3722112.104706/1‑100 (MQ=60)
TATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAAC < SRR3722112.633450/100‑1 (MQ=60)
ATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACC < SRR3722112.59236/100‑1 (MQ=60)
GCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAG > SRR3722112.325075/1‑100 (MQ=60)
ATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCA > SRR3722112.367378/1‑100 (MQ=60)
ACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCctgtctcttatacacatctgacgctgccgacgaatgtcaatgtg > SRR3722112.542473/1‑56 (MQ=60)
CGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGT > SRR3722112.154095/1‑100 (MQ=60)
CCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGAT > SRR3722112.556861/1‑100 (MQ=60)
|
GTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGAT > NZ_CP009273/2814171‑2814356
|
| 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: ‑ |