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 I205 R1
|
220 |
17.5 |
961422 |
97.1% |
933540 |
86.5 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
4,306,692 |
A→T |
I115N (ATC→AAC) |
phnM ← |
alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 4,306,692 | 0 | A | T | 92.3%
| 36.6
/ ‑3.3
| 13 | I115N (ATC→AAC) | phnM | alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base T (5/7); total (5/8) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 |
ATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGAT > NZ_CP009273/4306613‑4306771
|
aTGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTtctc > 1:248562/1‑90 (MQ=255)
gCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGGTGTTCATCTTCTCCa > 2:314953/1‑90 (MQ=255)
cGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTCCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAg < 2:335542/90‑1 (MQ=255)
cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATc < 2:248562/90‑1 (MQ=255)
cGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTTTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcc > 2:384915/1‑89 (MQ=255)
cgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa > 1:79581/1‑90 (MQ=255)
cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 1:143396/90‑1 (MQ=255)
cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 2:131212/90‑1 (MQ=255)
cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa > 1:144739/1‑90 (MQ=255)
cgTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa < 2:201074/90‑1 (MQ=255)
ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc < 2:41487/90‑1 (MQ=255)
tGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGt < 1:15500/90‑1 (MQ=255)
gATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGAt < 2:79581/90‑1 (MQ=255)
|
ATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGAT > NZ_CP009273/4306613‑4306771
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
CAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGA > NZ_CP009273/4306601‑4306770
|
gggctcggagatgtgtataagagacagGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTC < SRR3722087.426327/73‑1 (MQ=60)
GCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTC > SRR3722087.250961/1‑100 (MQ=60)
gagattcgccttagtctcgtgggctcgagatgtgtataagagacaGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCG < SRR3722087.122310/55‑1 (MQ=60)
GTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCA > SRR3722087.80320/1‑100 (MQ=60)
GTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCA > SRR3722087.146073/1‑100 (MQ=60)
CGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAC < SRR3722087.144722/100‑1 (MQ=60)
TGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGA < SRR3722087.15623/100‑1 (MQ=60)
|
CAGCGTGGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGA > NZ_CP009273/4306601‑4306770
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
