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 I199 R1
|
220 |
57.5 |
3163680 |
97.0% |
3068769 |
86.1 |
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.9%
| 145.7
/ ‑1.0
| 49 | E237G (GAA→GGA) | glyQ | glycine‑‑tRNA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (26/21); total (26/23) |
| Fisher's exact test for biased strand distribution p-value = 2.15e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01 |
GAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717885‑3718051
|
gagaTGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAgc > 1:915053/1‑90 (MQ=255)
gaTGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAgcgc > 2:614135/1‑90 (MQ=255)
tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGGTTTCCCAGCGCca > 2:841121/1‑90 (MQ=255)
tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 2:1099307/1‑90 (MQ=255)
tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 2:1529341/1‑90 (MQ=255)
gCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCcagc > 2:436668/1‑90 (MQ=255)
ttACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 2:408907/1‑85 (MQ=255)
ttACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca < 1:408907/85‑1 (MQ=255)
cgcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGctgctg > 2:1540013/1‑90 (MQ=255)
gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCtgctgc > 2:1117084/1‑90 (MQ=255)
gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCtgctgc < 1:712349/90‑1 (MQ=255)
gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCtgctgc > 1:239437/1‑90 (MQ=255)
gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCtgctgc > 1:138553/1‑90 (MQ=255)
cATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTgcg > 1:643945/1‑90 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATCCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:88129/90‑1 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGcttc < 1:166052/90‑1 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:1538587/90‑1 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:1441599/90‑1 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:1388077/90‑1 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:1353048/90‑1 (MQ=255)
gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTtctc > 2:400734/1‑90 (MQ=255)
ttGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCg > 1:38882/1‑90 (MQ=255)
aaGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAc < 1:577323/90‑1 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc < 1:868137/90‑1 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:410661/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc < 2:12467/90‑1 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:1131393/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:91788/1‑90 (MQ=255)
ggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCa < 1:168575/90‑1 (MQ=255)
ggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCa > 2:602310/1‑90 (MQ=255)
gAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa < 1:1117084/90‑1 (MQ=255)
gAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa < 1:904049/90‑1 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAgca < 1:1144355/90‑3 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 2:1243690/90‑1 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa > 2:1358191/1‑90 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 1:436668/90‑1 (MQ=255)
aCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGaa > 1:194294/1‑90 (MQ=255)
aCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGaa > 2:354096/1‑90 (MQ=255)
cGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAg > 2:984028/1‑90 (MQ=255)
gTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCcac > 2:1000217/1‑90 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 1:1290579/1‑90 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 2:852358/1‑90 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAATCCACATCcgcg < 1:1131393/90‑1 (MQ=255)
gcggcagcggATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTAtt < 1:498684/90‑1 (MQ=255)
cggcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTc < 1:1131337/90‑1 (MQ=255)
ggcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCg > 2:786503/1‑90 (MQ=255)
gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 2:1056435/90‑1 (MQ=255)
gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 1:1099307/90‑1 (MQ=255)
cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 1:786503/90‑1 (MQ=255)
|
GAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717885‑3718051
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
CTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTCCA > NZ_CP009273/3717875‑3718068
|
CTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGC > SRR3722078.927471/1‑100 (MQ=60)
TGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTG < SRR3722078.413852/100‑1 (MQ=60)
GGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGC > SRR3722078.140164/1‑100 (MQ=60)
GGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGC > SRR3722078.242359/1‑100 (MQ=60)
GCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCG > SRR3722078.651342/1‑100 (MQ=60)
GCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCT < SRR3722078.721073/100‑1 (MQ=60)
ATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCG > SRR3722078.39350/1‑100 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATCCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722078.89169/100‑1 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722078.167984/100‑1 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722078.1374767/100‑1 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722078.1465072/100‑1 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722078.1563670/100‑1 (MQ=60)
CAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGCAC < SRR3722078.1410501/100‑1 (MQ=60)
AAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGC < SRR3722078.584055/100‑1 (MQ=60)
GTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTG < SRR3722078.879758/100‑1 (MQ=60)
GGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGG < SRR3722078.170528/100‑1 (MQ=60)
CTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAA > SRR3722078.196614/1‑100 (MQ=60)
GAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCAC < SRR3722078.916173/100‑1 (MQ=60)
GAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCAC < SRR3722078.1134142/100‑1 (MQ=60)
TACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATC < SRR3722078.441910/100‑1 (MQ=60)
TACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGCAAGTCCACATC < SRR3722078.1161971/100‑1 (MQ=60)
CTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCG > SRR3722078.1311120/1‑100 (MQ=60)
GGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAATCCACATCCGCGTATTCGAAGT < SRR3722078.1148717/100‑1 (MQ=60)
GCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGATGTA < SRR3722078.504587/100‑1 (MQ=60)
CGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAA < SRR3722078.1148659/100‑1 (MQ=60)
GCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGT < SRR3722078.1116042/100‑1 (MQ=60)
CGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGAC < SRR3722078.796766/100‑1 (MQ=60)
CCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTCCA > SRR3722078.1210126/1‑100 (MQ=60)
|
CTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTCCA > NZ_CP009273/3717875‑3718068
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |