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 I199 R1
|
220 |
57.5 |
3163680 |
97.0% |
3068769 |
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 | 94.6%
| 106.2
/ ‑0.4
| 37 | F16S (TTC→TCC) | pspD | phage shock protein PspD |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (16/19); total (17/20) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01 |
TCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCT > NZ_CP009273/1363610‑1363781
|
tCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCaa > 2:579441/1‑90 (MQ=255)
tACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGc > 1:479067/1‑90 (MQ=255)
tACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGc > 1:98559/1‑90 (MQ=255)
aCGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCa < 1:1056816/90‑1 (MQ=255)
gCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCaa < 1:187372/90‑1 (MQ=255)
ccGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg < 2:479067/90‑1 (MQ=255)
gTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCt > 2:624565/1‑90 (MQ=255)
gTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCt < 2:1279194/90‑1 (MQ=255)
ttCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTgg < 2:1409687/90‑1 (MQ=255)
cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 2:754372/90‑1 (MQ=255)
cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 1:1346780/90‑1 (MQ=255)
cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 1:457053/90‑1 (MQ=255)
aCTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActtctt > 1:121748/1‑90 (MQ=255)
aCTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActtctt > 1:542423/1‑90 (MQ=255)
cTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTa > 2:1196623/1‑90 (MQ=255)
tgtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTAc < 1:714729/90‑1 (MQ=255)
tgtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTAc < 1:579441/90‑1 (MQ=255)
gAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGc < 1:173297/90‑1 (MQ=255)
ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 1:511519/90‑1 (MQ=255)
ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 1:1284839/90‑1 (MQ=255)
ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 1:1431264/90‑1 (MQ=255)
aaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTgc < 1:1407143/90‑1 (MQ=255)
ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCTAAATTAGCAGGCAAg > 1:524180/1‑65 (MQ=255)
ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCTAAATTAGCAGGCAAg < 2:524180/65‑1 (MQ=255)
ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg < 2:821860/65‑1 (MQ=255)
ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg > 1:821860/1‑65 (MQ=255)
gAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATgg > 1:360459/1‑90 (MQ=255)
cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc > 1:444662/1‑63 (MQ=255)
cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc < 2:444662/63‑1 (MQ=255)
cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCcgg > 1:640300/1‑90 (MQ=255)
tGGCAACAGGGCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCgggg > 2:1187425/1‑89 (MQ=255)
gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTa > 2:738616/1‑65 (MQ=255)
gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTa < 1:738616/65‑1 (MQ=255)
aaaGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt > 2:1326228/1‑37 (MQ=39)
aaaGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt < 1:1326228/37‑1 (MQ=39)
ggTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAACAGTTGCt > 2:472303/1‑90 (MQ=255)
ggTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGGGTGGCGGGGGGGGGGGTAAAAAAAGtttgt > 2:1018529/1‑87 (MQ=255)
|
TCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCT > NZ_CP009273/1363610‑1363781
|
| 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
CGATACTTTCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCG > NZ_CP009273/1363602‑1363748
|
CGATACTTTCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAAT < SRR3722078.1169020/100‑1 (MQ=60)
ACTTTCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGC > SRR3722078.99716/1‑100 (MQ=60)
ACTTTCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGC > SRR3722078.484746/1‑100 (MQ=60)
ACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGG < SRR3722078.1072696/100‑1 (MQ=60)
GCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTT < SRR3722078.189598/100‑1 (MQ=60)
GTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTT > SRR3722078.123167/1‑100 (MQ=60)
GTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTT > SRR3722078.548830/1‑100 (MQ=60)
CGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGC < SRR3722078.462509/100‑1 (MQ=60)
CGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGC < SRR3722078.1368401/100‑1 (MQ=60)
acaaccgccatcgtcggcagcgtcagatgtgtataagagacaGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAC < SRR3722078.1347431/58‑1 (MQ=60)
TGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGC < SRR3722078.586202/100‑1 (MQ=60)
TGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGC < SRR3722078.723511/100‑1 (MQ=60)
GTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCTAAATTAGCAGGCAAGCTGGTACTTCctgtctcttatacac > SRR3722078.530380/1‑85 (MQ=60)
GTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCctgtctcttatacac > SRR3722078.832733/1‑85 (MQ=60)
GAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTAT < SRR3722078.175330/100‑1 (MQ=60)
GGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGG < SRR3722078.1454538/100‑1 (MQ=60)
GGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGG > SRR3722078.364847/1‑100 (MQ=60)
GGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGG < SRR3722078.1305277/100‑1 (MQ=60)
GGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGG < SRR3722078.517570/100‑1 (MQ=60)
AAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCC < SRR3722078.1429962/100‑1 (MQ=60)
GTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCctgtctcttatacacat > SRR3722078.449985/1‑83 (MQ=60)
GTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGG > SRR3722078.647667/1‑100 (MQ=60)
gtgtataagagacaGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCG < SRR3722078.747940/86‑1 (MQ=60)
|
CGATACTTTCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCG > NZ_CP009273/1363602‑1363748
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |