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 I200 R1
|
216 |
30.8 |
1683622 |
97.2% |
1636480 |
87.0 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
2,859,998 |
T→C |
I305V (ATT→GTT) |
rpoS ← |
RNA polymerase sigma factor RpoS |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,859,998 | 0 | T | C | 83.9%
| 62.9
/ 9.1
| 31 | I305V (ATT→GTT) | rpoS | RNA polymerase sigma factor RpoS |
| Reads supporting (aligned to +/- strand): ref base T (1/4); new base C (11/15); total (12/19) |
| Fisher's exact test for biased strand distribution p-value = 6.24e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01 |
ACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC > NZ_CP009273/2859916‑2860084
|
acttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa < 2:595550/90‑1 (MQ=255)
acttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGa < 2:36845/90‑1 (MQ=255)
acttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGa < 2:670189/90‑1 (MQ=255)
tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca < 1:334007/90‑1 (MQ=255)
ctCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACg > 1:426847/1‑90 (MQ=255)
ggAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCAc < 1:431328/90‑1 (MQ=255)
aGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTg < 1:351129/90‑1 (MQ=255)
tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCa > 2:814279/1‑90 (MQ=255)
tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:678269/90‑1 (MQ=255)
tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:458469/90‑1 (MQ=255)
cAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg > 2:489886/1‑90 (MQ=255)
gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc > 2:774462/1‑90 (MQ=255)
gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc > 2:470530/1‑90 (MQ=255)
ccccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc < 2:531873/90‑1 (MQ=255)
tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:812106/82‑1 (MQ=255)
tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc > 2:812106/1‑82 (MQ=255)
tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 1:449951/90‑1 (MQ=255)
tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 1:398953/90‑1 (MQ=255)
ttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc > 1:244200/1‑90 (MQ=255)
tCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc > 1:702725/1‑90 (MQ=255)
tCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc < 2:470679/90‑1 (MQ=255)
gcgcAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc < 2:89787/90‑1 (MQ=255)
aaCGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGt < 1:558268/90‑1 (MQ=255)
gcgcAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAccc > 2:162262/1‑90 (MQ=255)
gcgcAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAccc > 2:87098/1‑90 (MQ=255)
cAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAg < 1:368052/90‑1 (MQ=255)
tCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCg > 1:828945/1‑90 (MQ=255)
tCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCg > 2:837317/1‑90 (MQ=255)
gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc < 2:684022/61‑1 (MQ=255)
gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc < 2:682271/90‑1 (MQ=255)
gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc < 1:246636/90‑1 (MQ=255)
|
ACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC > NZ_CP009273/2859916‑2860084
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
ATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCAC > NZ_CP009273/2859910‑2860094
|
ATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACAC < SRR3722079.76544/100‑1 (MQ=60)
TNCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACG > SRR3722079.572182/1‑100 (MQ=60)
TGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACG > SRR3722079.430532/1‑100 (MQ=60)
TACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGG < SRR3722079.336658/100‑1 (MQ=60)
GGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCC < SRR3722079.435063/100‑1 (MQ=60)
AGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTT < SRR3722079.353893/100‑1 (MQ=60)
TCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATC < SRR3722079.820813/100‑1 (MQ=60)
TCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATC < SRR3722079.685389/100‑1 (MQ=60)
TCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATC < SRR3722079.462479/100‑1 (MQ=60)
GCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTC > SRR3722079.246257/1‑100 (MQ=60)
TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTG < SRR3722079.402305/100‑1 (MQ=60)
TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTG < SRR3722079.453889/100‑1 (MQ=60)
TTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACctgtctcttatacacatctccgagcccacgagaccagagaggatct > SRR3722079.424605/1‑54 (MQ=60)
TTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCC > SRR3722079.710169/1‑100 (MQ=60)
AACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAA < SRR3722079.563670/100‑1 (MQ=60)
GCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCG > SRR3722079.837818/1‑100 (MQ=60)
CAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAAT < SRR3722079.370973/100‑1 (MQ=60)
CCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCtgtctcttatacacatct > SRR3722079.691213/1‑82 (MQ=60)
ACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGctgtc > SRR3722079.522582/1‑95 (MQ=60)
ACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCA > SRR3722079.383/1‑100 (MQ=60)
GAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCAC < SRR3722079.248707/100‑1 (MQ=60)
|
ATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCAC > NZ_CP009273/2859910‑2860094
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
