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,536,106 |
T→C |
E261G (GAA→GGA) |
cysP ← |
thiosulfate/sulfate ABC transporter substrate‑binding protein CysP |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,536,106 | 0 | T | C | 90.0%
| 48.1
/ ‑0.0
| 20 | E261G (GAA→GGA) | cysP | thiosulfate/sulfate ABC transporter substrate‑binding protein CysP |
| Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (14/4); total (14/6) |
| Fisher's exact test for biased strand distribution p-value = 7.89e-02 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01 |
GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA > NZ_CP009273/2536021‑2536190
|
gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGt < 1:396620/90‑1 (MQ=255)
tATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAc > 2:499550/1‑90 (MQ=255)
aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc > 1:424204/1‑90 (MQ=255)
aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc > 1:72916/1‑90 (MQ=255)
cgcgGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGGGGCTTTTCCCGTACCCTTgg > 1:50472/1‑90 (MQ=255)
agtaATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCa < 1:486762/90‑1 (MQ=255)
tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:303279/1‑90 (MQ=255)
tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:93549/1‑90 (MQ=255)
tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:611182/1‑90 (MQ=255)
tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:279481/1‑90 (MQ=255)
tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:215404/1‑90 (MQ=255)
tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:57025/1‑90 (MQ=255)
cGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTAt > 2:307690/1‑90 (MQ=255)
ggTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGc < 2:276059/90‑1 (MQ=255)
tGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAt < 2:700722/90‑1 (MQ=255)
aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg > 2:320656/1‑90 (MQ=255)
aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg > 2:60812/1‑90 (MQ=255)
aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg > 2:684206/1‑90 (MQ=255)
aGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGaa < 1:279481/90‑1 (MQ=255)
cTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCaa < 2:132323/90‑1 (MQ=255)
|
GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA > NZ_CP009273/2536021‑2536190
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ 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
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTC > NZ_CP009273/2536008‑2536203
|
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCC < SRR3722112.68621/100‑1 (MQ=60)
ACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAC < SRR3722112.222249/100‑1 (MQ=60)
ACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAC < SRR3722112.340349/100‑1 (MQ=60)
CCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACC > SRR3722112.429864/1‑100 (MQ=60)
CCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACC > SRR3722112.73803/1‑100 (MQ=60)
GGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGGGGCTTTTCCCGTACCCTTGG > SRR3722112.51088/1‑100 (MQ=60)
GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCC < SRR3722112.401828/100‑1 (MQ=60)
ATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATGGGCTTTGGCGGCTTTTCCCGTACCtgtctcttat > SRR3722112.180296/1‑90 (MQ=60)
GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA > SRR3722112.306764/1‑100 (MQ=60)
GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA > SRR3722112.57724/1‑100 (MQ=60)
AGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATC < SRR3722112.493776/100‑1 (MQ=60)
AGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTC < SRR3722112.282729/100‑1 (MQ=60)
TTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCtgtctctta > SRR3722112.117035/1‑91 (MQ=60)
|
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCCTGCGCTTC > NZ_CP009273/2536008‑2536203
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
