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,174,381 |
T→G |
C270G (TGC→GGC) |
yegU → |
ADP‑ribosylglycohydrolase family protein |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,174,381 | 0 | T | G | 93.3%
| 91.3
/ ‑0.4
| 30 | C270G (TGC→GGC) | yegU | ADP‑ribosylglycohydrolase family protein |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base G (17/11); total (18/12) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01 |
GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATG > NZ_CP009273/2174297‑2174462
|
gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc > 2:211956/1‑90 (MQ=255)
gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc > 1:388571/1‑90 (MQ=255)
tCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCATAAACCGACCCGAATCGCGGCGCCg < 1:15033/90‑1 (MQ=255)
tCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCg < 2:90156/90‑1 (MQ=255)
gTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTc > 2:61313/1‑90 (MQ=255)
gcAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTgc > 2:458393/1‑90 (MQ=255)
cAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCt > 2:637631/1‑90 (MQ=255)
gAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGac < 2:388571/90‑1 (MQ=255)
gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac > 1:485809/1‑90 (MQ=255)
gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac > 1:658681/1‑90 (MQ=255)
gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac > 2:60351/1‑90 (MQ=255)
gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac > 2:373222/1‑90 (MQ=255)
cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa > 2:625047/1‑90 (MQ=255)
cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa > 2:289575/1‑90 (MQ=255)
cgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCa < 1:61313/90‑1 (MQ=255)
gtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCAt > 2:497607/1‑90 (MQ=255)
tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc < 1:637631/90‑1 (MQ=255)
tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc < 1:142458/90‑1 (MQ=255)
cgcCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCgg < 2:523756/90‑1 (MQ=255)
gcgcTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGgcggcg > 2:510859/1‑67 (MQ=255)
gcgcTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGgcggcg < 1:510859/67‑1 (MQ=255)
gCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc < 1:16479/66‑1 (MQ=255)
gCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc > 2:16479/1‑66 (MQ=255)
gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc > 2:229001/1‑90 (MQ=255)
gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc > 2:477848/1‑90 (MQ=255)
acaAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgcgg < 2:187037/90‑1 (MQ=255)
aaCCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGgcgc > 2:18726/1‑90 (MQ=255)
ccGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGt < 1:625047/90‑1 (MQ=255)
gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTg > 1:546335/1‑90 (MQ=255)
cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg < 1:152335/90‑1 (MQ=255)
|
GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATG > NZ_CP009273/2174297‑2174462
|
| 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
AGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCGA > NZ_CP009273/2174286‑2174478
|
ccgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacaGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGC < SRR3722112.660922/47‑1 (MQ=60)
GCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCC > SRR3722112.393635/1‑100 (MQ=60)
TCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCATAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGC < SRR3722112.15225/100‑1 (MQ=60)
GCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACAC > SRR3722112.492802/1‑100 (MQ=60)
GCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACAC > SRR3722112.668961/1‑100 (MQ=60)
gtataagagacagGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACAC < SRR3722112.518322/87‑1 (MQ=60)
CGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTAT < SRR3722112.62052/100‑1 (MQ=60)
tgtataagagacagGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG < SRR3722112.16686/86‑1 (MQ=60)
tcgtcggcagcgtcagatgtgtataagagacagGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG < SRR3722112.27060/67‑1 (MQ=60)
TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG < SRR3722112.144242/100‑1 (MQ=60)
TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG < SRR3722112.647538/100‑1 (MQ=60)
GGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTG > SRR3722112.554553/1‑100 (MQ=60)
CCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGT < SRR3722112.634697/100‑1 (MQ=60)
CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC < SRR3722112.154214/100‑1 (MQ=60)
GCGGCGCCGGCCTGTGCGCTAACCTTGGCGGGGACACCGGCACCATCGGTGCTATGGCGACGTCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCGA < SRR3722112.292928/100‑1 (MQ=60)
|
AGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCGA > NZ_CP009273/2174286‑2174478
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |