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 I200 R1
|
216 |
30.8 |
1683622 |
97.2% |
1636480 |
87.0 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
1,325,886 |
T→C |
G194G (GGT→GGC) |
topA → |
type I DNA topoisomerase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,325,886 | 0 | T | C | 92.9%
| 33.2
/ ‑3.2
| 14 | G194G (GGT→GGC) | topA | type I DNA topoisomerase |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (9/4); total (9/5) |
Fisher's exact test for biased strand distribution p-value = 3.57e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01 |
CGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGT > NZ_CP009273/1325804‑1325959
|
cgcgTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGc > 1:192686/1‑90 (MQ=255)
cgTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAg < 2:125250/90‑1 (MQ=255)
tttATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTg < 1:75552/90‑1 (MQ=255)
cgcgTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTg < 1:11482/81‑1 (MQ=255)
cgcgTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTg > 2:11482/1‑81 (MQ=255)
cgcgTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGc > 2:698720/1‑90 (MQ=255)
gTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCgagcgtgag > 1:699676/1‑90 (MQ=255)
tttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCgagcgtgagc > 1:765893/1‑90 (MQ=255)
ttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCgagcgtgagcg > 2:190476/1‑90 (MQ=255)
ttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCgagcgtgagcg > 2:504375/1‑90 (MQ=255)
ctgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGGGAAAtt > 2:166586/1‑90 (MQ=255)
gATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGcc < 1:393086/90‑1 (MQ=255)
gATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGcc < 2:479834/90‑1 (MQ=255)
gTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGt > 2:425208/1‑90 (MQ=255)
|
CGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGT > NZ_CP009273/1325804‑1325959
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
GTTAATGCCCAGCAGGCGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTC > NZ_CP009273/1325788‑1325961
|
GTTAATGCCCAGCAGGCGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCC > SRR3722079.504629/1‑100 (MQ=60)
GCCCAGCAGGCGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGC > SRR3722079.194355/1‑100 (MQ=60)
TTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCC < SRR3722079.11568/100‑1 (MQ=60)
TTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCC < SRR3722079.76171/100‑1 (MQ=60)
GGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAG > SRR3722079.707079/1‑100 (MQ=60)
GGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGC > SRR3722079.774089/1‑100 (MQ=60)
GATCGCTCGTGGCCTGTCTGCCGGCCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTC < SRR3722079.396362/100‑1 (MQ=60)
|
GTTAATGCCCAGCAGGCGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTC > NZ_CP009273/1325788‑1325961
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |