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 I207 R1
|
218 |
45.0 |
2482936 |
97.1% |
2410930 |
86.1 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
1,363,696 |
T→C |
F16S (TTC→TCC) |
pspD → |
phage shock protein PspD |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 1,363,696 | 0 | T | C | 96.2%
| 81.8
/ ‑4.6
| 26 | F16S (TTC→TCC) | pspD | phage shock protein PspD |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (13/12); total (13/13) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTT > NZ_CP009273/1363610‑1363778
|
tCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCaa < 2:1202456/90‑1 (MQ=255)
aCGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCa > 2:349155/1‑90 (MQ=255)
tcgCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGc < 1:568704/88‑1 (MQ=255)
aGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCa > 1:938429/1‑90 (MQ=255)
ttCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTgg > 2:439405/1‑90 (MQ=255)
cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAACGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 1:654740/90‑1 (MQ=255)
cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 2:41544/90‑1 (MQ=255)
cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 2:235390/90‑1 (MQ=255)
tCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt < 1:349155/90‑1 (MQ=255)
tCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt < 1:402681/90‑1 (MQ=255)
gtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTAcc > 1:1123626/1‑90 (MQ=255)
gAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGc < 1:1211681/70‑1 (MQ=255)
gAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGc > 2:1211681/1‑70 (MQ=255)
ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 1:404878/90‑1 (MQ=255)
ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 1:1225445/90‑1 (MQ=255)
aaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTgc > 2:324899/1‑90 (MQ=255)
ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCt > 2:560084/1‑90 (MQ=255)
ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCt > 1:981304/1‑90 (MQ=255)
tATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTa < 1:2368/90‑1 (MQ=255)
gAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATgg > 2:461980/1‑90 (MQ=255)
cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCcgg > 1:745290/1‑90 (MQ=255)
cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCcgg > 1:483217/1‑90 (MQ=255)
gCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCcggcgg > 1:1115501/1‑90 (MQ=255)
cTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCggg < 1:174328/90‑1 (MQ=255)
tGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCgggg > 1:369714/1‑89 (MQ=255)
ccTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGtt < 1:1053806/90‑1 (MQ=255)
|
TCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTT > NZ_CP009273/1363610‑1363778
|
| 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
CACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGAC > NZ_CP009273/1363611‑1363788
|
CACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCA > SRR3722089.952160/1‑100 (MQ=60)
taagagacagCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAC < SRR3722089.1230109/90‑1 (MQ=60)
TCGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAC < SRR3722089.575421/100‑1 (MQ=60)
TCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACC > SRR3722089.1140675/1‑100 (MQ=60)
CGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAACGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGC < SRR3722089.663158/100‑1 (MQ=60)
TCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAC < SRR3722089.407569/100‑1 (MQ=60)
TCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAC < SRR3722089.353491/100‑1 (MQ=60)
GTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCT > SRR3722089.995833/1‑100 (MQ=60)
GGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGG < SRR3722089.409790/100‑1 (MQ=60)
GGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGG < SRR3722089.1244071/100‑1 (MQ=60)
GTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGG > SRR3722089.488957/1‑100 (MQ=60)
GTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGG > SRR3722089.755113/1‑100 (MQ=60)
TATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCG < SRR3722089.2390/100‑1 (MQ=60)
ATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGG > SRR3722089.1132396/1‑100 (MQ=60)
GAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGG > SRR3722089.374257/1‑100 (MQ=60)
CTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGC < SRR3722089.176509/100‑1 (MQ=60)
GGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGGGTGGCGGGgggggggggaaaaaaaatttgt > SRR3722089.501223/1‑78 (MQ=60)
GTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGgggggcggggttggggggaaaaacaaatttgc > SRR3722089.703467/1‑68 (MQ=60)
CCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGAC < SRR3722089.1069616/100‑1 (MQ=60)
|
CACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGAC > NZ_CP009273/1363611‑1363788
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |