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 |
3,579,168 |
T→C |
E339G (GAA→GGA) |
ggt ← |
gamma‑glutamyltransferase |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 3,579,168 | 0 | T | C | 100.0%
| 31.2
/ NA
| 10 | E339G (GAA→GGA) | ggt | gamma‑glutamyltransferase |
| Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (4/6); total (4/6) |
ATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATC > NZ_CP009273/3579089‑3579253
|
aTCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc < 1:282028/90‑1 (MQ=255)
gATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggcgtaggcg > 1:200150/1‑90 (MQ=255)
tAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTcc < 1:45143/90‑1 (MQ=255)
ggCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCAGc > 2:336430/1‑90 (MQ=255)
tttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc < 2:238334/90‑1 (MQ=255)
tttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc < 2:64009/90‑1 (MQ=255)
ttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCt > 2:192538/1‑90 (MQ=255)
tGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCa < 2:141011/90‑1 (MQ=255)
tGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatc > 2:286418/1‑90 (MQ=255)
atatCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATc < 1:286418/90‑1 (MQ=255)
|
ATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATC > NZ_CP009273/3579089‑3579253
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
ATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTTTCC > NZ_CP009273/3579089‑3579263
|
ATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCG > SRR3722090.202359/1‑100 (MQ=60)
ATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCG < SRR3722090.285546/100‑1 (MQ=60)
TAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCA < SRR3722090.45609/100‑1 (MQ=60)
AAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTT > SRR3722090.72778/1‑100 (MQ=60)
ATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTTTCC < SRR3722090.289997/100‑1 (MQ=60)
|
ATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTCATATCGAAGTTTTCC > NZ_CP009273/3579089‑3579263
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
