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 |
3,717,969 |
T→C |
E237G (GAA→GGA) |
glyQ ← |
glycine‑‑tRNA ligase subunit alpha |
|
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 3,717,969 | 0 | T | C | 95.8%
| 69.2
/ ‑3.3
| 24 | E237G (GAA→GGA) | glyQ | glycine‑‑tRNA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (13/10); total (14/10) |
| 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.68e-01 |
GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717887‑3718051
|
gaTGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAgcgc > 2:692580/1‑90 (MQ=255)
cATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTgcg > 2:716898/1‑90 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:483699/90‑1 (MQ=255)
gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTtctc > 2:620007/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:125362/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:350241/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 1:40810/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc < 1:365953/90‑1 (MQ=255)
gggCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGa > 2:434652/1‑90 (MQ=255)
ggcggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGaa < 1:643017/90‑1 (MQ=255)
cTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAgg < 1:704673/90‑1 (MQ=255)
cAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGaa > 2:398757/1‑90 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 1:711992/90‑1 (MQ=255)
cGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAg > 1:448137/1‑90 (MQ=255)
gTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCcac > 2:641329/1‑90 (MQ=255)
tGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgc < 1:263588/90‑1 (MQ=255)
ggcagcggcagcGGATTTCCCAGCTCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 2:285274/1‑90 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 2:22952/1‑90 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 2:645268/1‑90 (MQ=255)
ggcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCg > 2:166842/1‑90 (MQ=255)
gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 2:229342/90‑1 (MQ=255)
cagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGaa < 1:125362/90‑1 (MQ=255)
cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 2:448137/90‑1 (MQ=255)
cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 1:184514/90‑1 (MQ=255)
|
GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717887‑3718051
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
CGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGAC > NZ_CP009273/3717877‑3718061
|
CGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGC < SRR3722112.597577/100‑1 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722112.490650/100‑1 (MQ=60)
GGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTC > SRR3722112.41302/1‑100 (MQ=60)
GTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTG < SRR3722112.370653/100‑1 (MQ=60)
GGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAAC < SRR3722112.653022/100‑1 (MQ=60)
CTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAA < SRR3722112.715675/100‑1 (MQ=60)
TTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAG > SRR3722112.454374/1‑100 (MQ=60)
TACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATC < SRR3722112.723121/100‑1 (MQ=60)
TGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAG < SRR3722112.266689/100‑1 (MQ=60)
CAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTG < SRR3722112.126917/100‑1 (MQ=60)
CGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGAC < SRR3722112.186775/100‑1 (MQ=60)
|
CGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGAC > NZ_CP009273/3717877‑3718061
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: ‑ |
