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 I200 R1
|
216 |
30.8 |
1683622 |
97.2% |
1636480 |
87.0 |
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 | 86.8%
| 86.3
/ 8.2
| 38 | E237G (GAA→GGA) | glyQ | glycine‑‑tRNA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base T (3/2); new base C (17/16); total (20/18) |
| 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.22e-01 |
GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717887‑3718051
|
gaTGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAgcgc > 2:301853/1‑90 (MQ=255)
tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 1:199058/1‑90 (MQ=255)
tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 1:653015/1‑90 (MQ=255)
ttACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCagcagc > 1:78719/1‑90 (MQ=255)
gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCtgctgc > 2:781733/1‑90 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATGCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:718390/90‑1 (MQ=255)
cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:810707/90‑1 (MQ=255)
gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc > 1:551131/1‑84 (MQ=255)
gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:551131/84‑1 (MQ=255)
gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTtctc > 1:801587/1‑90 (MQ=255)
aaGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAc < 2:653015/90‑1 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:631718/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:67112/1‑90 (MQ=255)
gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 1:689769/1‑90 (MQ=255)
gggCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGa > 1:387456/1‑90 (MQ=255)
cggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc > 1:523339/1‑90 (MQ=255)
cggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc < 2:517303/90‑1 (MQ=255)
gcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAg > 1:163107/1‑90 (MQ=255)
ttCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTg > 2:327349/1‑90 (MQ=255)
aaTACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAg < 2:501346/90‑1 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 2:288569/90‑1 (MQ=255)
tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 2:10334/90‑1 (MQ=255)
gTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc > 1:631754/1‑71 (MQ=255)
gTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc < 2:631754/71‑1 (MQ=255)
gTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCcac > 2:712828/1‑90 (MQ=255)
aGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACAt < 2:689769/90‑1 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc > 2:340245/1‑44 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:445254/44‑1 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:340245/44‑1 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc > 1:445254/1‑44 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa > 2:825222/1‑73 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa < 1:825222/73‑1 (MQ=255)
ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 1:121735/1‑90 (MQ=255)
gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 1:298582/90‑1 (MQ=255)
gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 2:265751/90‑1 (MQ=255)
cggATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 2:523339/90‑1 (MQ=255)
cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 1:508479/90‑1 (MQ=255)
cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 2:45414/90‑1 (MQ=255)
|
GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717887‑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
TGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTCC > NZ_CP009273/3717870‑3718067
|
TGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTC < SRR3722079.592352/100‑1 (MQ=60)
GTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCA > SRR3722079.200776/1‑100 (MQ=60)
GTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCA > SRR3722079.659775/1‑100 (MQ=60)
GGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGC > SRR3722079.79360/1‑100 (MQ=60)
CATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTC > SRR3722079.810182/1‑100 (MQ=60)
CATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTC > SRR3722079.556414/1‑100 (MQ=60)
tggtctcgtgggctcggagatgtgtataagagacaGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAC < SRR3722079.342929/65‑1 (MQ=60)
GGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTC > SRR3722079.697020/1‑100 (MQ=60)
TTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGA > SRR3722079.390651/1‑100 (MQ=60)
AAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGC > SRR3722079.528186/1‑100 (MQ=60)
GCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAG > SRR3722079.164538/1‑100 (MQ=60)
gagacaGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCAC < SRR3722079.834057/94‑1 (MQ=60)
GAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGctgtctctt > SRR3722079.638215/1‑91 (MQ=60)
CTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACctgtctcttatacacatctccgagcccacgagacca > SRR3722079.449143/1‑64 (MQ=60)
CTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCG > SRR3722079.122734/1‑100 (MQ=60)
GCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGT < SRR3722079.301029/100‑1 (MQ=60)
CGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGAC < SRR3722079.513101/100‑1 (MQ=60)
TTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTC < SRR3722079.322608/100‑1 (MQ=60)
TCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTCC > SRR3722079.333739/1‑100 (MQ=60)
|
TGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTTGTAAGTGGACTGCTCC > NZ_CP009273/3717870‑3718067
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |