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 I231 R1
|
218 |
12.8 |
703680 |
97.2% |
683976 |
86.9 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
1,807,028 |
T→C |
E121G (GAA→GGA) |
BW25113_RS09050 ← |
hypothetical protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,807,028 | 0 | T | C | 100.0%
| 39.1
/ NA
| 13 | E121G (GAA→GGA) | BW25113_RS09050 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (7/6); total (7/6) |
ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA > NZ_CP009273/1806947‑1807108
|
atataATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGa < 1:223594/90‑1 (MQ=255)
taATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAAt < 2:207476/90‑1 (MQ=255)
ttGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGa < 1:87443/90‑1 (MQ=255)
taGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTc < 1:293163/90‑1 (MQ=255)
ccTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAAt > 2:332097/1‑90 (MQ=255)
cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa > 1:27959/1‑90 (MQ=255)
cTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATa > 2:201280/1‑90 (MQ=255)
tAAACAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATaaa > 1:8999/1‑90 (MQ=255)
tGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCAtt < 1:18586/75‑1 (MQ=255)
tGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCAtt > 2:18586/1‑75 (MQ=255)
aaccaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTttatttat > 1:258890/1‑90 (MQ=255)
accaaccTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTtata < 1:258543/90‑1 (MQ=255)
gTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGa > 1:315516/1‑90 (MQ=255)
|
ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA > NZ_CP009273/1806947‑1807108
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA > NZ_CP009273/1806947‑1807108
|
ggcgaccaccgagatctacacagccaacttcgtcggcagcgtcagatgtgtataagagacagGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAG < SRR3722116.352492/38‑1 (MQ=60)
ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAG < SRR3722116.226566/100‑1 (MQ=60)
TTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCC < SRR3722116.88529/100‑1 (MQ=60)
GTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATA > SRR3722116.28305/1‑100 (MQ=60)
ATATAGTCCTTAAACAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAA > SRR3722116.9113/1‑100 (MQ=60)
gacagCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAAT < SRR3722116.18836/95‑1 (MQ=60)
TAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATC < SRR3722116.297451/100‑1 (MQ=60)
AATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTAT > SRR3722116.262580/1‑100 (MQ=60)
ACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCC < SRR3722116.262220/100‑1 (MQ=60)
GATCAAATTGGTCATTTTTCCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA > SRR3722116.320195/1‑100 (MQ=60)
|
ATATAATTTTGTTCATAGTATATAGTCCTTAAATAGTGAATAACCAACCTGTTTTTTATACTGATCAAATTGGTCATTTTTTCTGGGCGAAATACAGAAGTGACATCCATTCCAATAAAATCTTTATTTATACCTCTTATCCGAATTTCACTACCATAAAGA > NZ_CP009273/1806947‑1807108
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |