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 I227 R1
|
223 |
21.3 |
1167702 |
97.0% |
1132670 |
86.4 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
2,670,359 |
T→C |
G454G (GGA→GGG) |
yphE ← |
sugar ABC transporter ATP‑binding protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,670,359 | 0 | T | C | 87.5%
| 33.0
/ 0.3
| 16 | G454G (GGA→GGG) | yphE | sugar ABC transporter ATP‑binding protein |
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (9/5); total (10/6) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 |
AATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAA > NZ_CP009273/2670291‑2670445
|
aataaCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACCCCACGCTTTTCCCTTCGGCAGCCAGCTCACgg < 2:36850/90‑1 (MQ=255)
cAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt < 1:179572/90‑1 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACg > 2:325705/1‑90 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGAc < 2:526740/84‑1 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGAc > 1:526740/1‑84 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg < 1:406834/90‑1 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg > 2:479199/1‑90 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg > 2:231013/1‑90 (MQ=255)
gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg > 2:365718/1‑90 (MQ=255)
acacacCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATctg > 2:403039/1‑90 (MQ=255)
aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg < 1:325705/90‑1 (MQ=255)
aaCTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGGTTGGCTTCg > 1:324310/1‑90 (MQ=255)
cctccACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGatat > 2:569228/1‑90 (MQ=255)
ctccACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATc > 2:119857/1‑90 (MQ=255)
cTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGAcgc < 1:307137/90‑1 (MQ=255)
acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc > 2:519583/1‑90 (MQ=255)
ttttCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGcaacaa > 1:265316/1‑90 (MQ=255)
|
AATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAA > NZ_CP009273/2670291‑2670445
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAA > NZ_CP009273/2670289‑2670445
|
GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG > SRR3722111.535005/1‑100 (MQ=60)
CAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATC < SRR3722111.181830/100‑1 (MQ=60)
GATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGC < SRR3722111.412772/100‑1 (MQ=60)
CACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGGTTGGCTTCG > SRR3722111.328629/1‑100 (MQ=60)
AGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATAT < SRR3722111.330049/100‑1 (MQ=60)
CTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGG < SRR3722111.311138/100‑1 (MQ=60)
ACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAA > SRR3722111.268498/1‑100 (MQ=60)
ACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAA > SRR3722111.30307/1‑100 (MQ=60)
|
GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAA > NZ_CP009273/2670289‑2670445
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |