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 |
661,892 |
A→C |
G594G (GGT→GGG) |
mrdA ← |
peptidoglycan DD‑transpeptidase MrdA |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 661,892 | 0 | A | C | 85.7%
| 39.5
/ 1.8
| 15 | coding (1782/1902 nt) | mrdA | peptidoglycan DD‑transpeptidase MrdA |
| Reads supporting (aligned to +/- strand): ref base A (1/0); major base C (8/4); minor base . (1/1); total (10/5) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01 |
TGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC > NZ_CP009273/661822‑661966
|
tgttgtTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt > 1:244806/1‑75 (MQ=255)
tgttgtTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt < 2:244806/75‑1 (MQ=255)
gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt > 1:97917/1‑83 (MQ=255)
gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt < 2:97917/83‑1 (MQ=255)
gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa > 2:81549/1‑90 (MQ=255)
gAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGTGGATTGTTGTATgg > 1:303338/1‑90 (MQ=255)
aaCCGCCGGACCCGCACC‑CCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg < 1:283109/53‑1 (MQ=255)
aaCCGCCGGACCCGCACC‑CCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg > 2:283109/1‑53 (MQ=255)
cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt > 1:26891/1‑53 (MQ=255)
cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt < 2:26891/53‑1 (MQ=255)
cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAGAGGCGGTCATCAGTTTGTGGTc > 1:10353/1‑90 (MQ=255)
cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc > 1:280003/1‑90 (MQ=255)
cgccgGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc > 2:316618/1‑90 (MQ=255)
ccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt > 1:243913/1‑59 (MQ=255)
ccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt < 2:243913/59‑1 (MQ=255)
|
TGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC > NZ_CP009273/661822‑661966
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ 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
GGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGC > NZ_CP009273/661812‑661981
|
GGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCtgtc > SRR3722090.247745/1‑96 (MQ=60)
gggctgggagatgtgtataagagacagTCAGTGTACCACCCGCCGGACCCGC‑CCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT < SRR3722090.286647/73‑1 (MQ=60)
GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA > SRR3722090.98974/1‑100 (MQ=60)
TAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGTGGATTGTTGTATGG > SRR3722090.307196/1‑100 (MQ=60)
GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCtgtctcttatacacatctccgagccc > SRR3722090.27163/1‑74 (MQ=60)
GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAGAGGCGGTCATCAGTTTGTGGTC > SRR3722090.10466/1‑100 (MQ=60)
GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC > SRR3722090.283485/1‑100 (MQ=60)
GTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCtgtctcttatacacatctcc > SRR3722090.246837/1‑80 (MQ=60)
GACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGC > SRR3722090.47629/1‑100 (MQ=60)
|
GGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGC > NZ_CP009273/661812‑661981
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
