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 I209 R1
|
216 |
14.4 |
815930 |
96.1% |
784108 |
85.2 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
1,838,813 |
T→C |
I106V (ATT→GTT) |
ynjI ← |
DUF1266 domain‑containing protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,838,813 | 0 | T | C | 90.9%
| 26.5
/ ‑4.7
| 11 | I106V (ATT→GTT) | ynjI | DUF1266 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (6/4); total (7/4) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAATACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAATTAATAGAT > NZ_CP009273/1838736‑1838896
|
aGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTg < 2:276143/90‑1 (MQ=255)
aTTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGataat < 1:172890/90‑1 (MQ=255)
aTTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGataat < 2:218042/90‑1 (MQ=255)
ttGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAATACCTTTAGCCTGataata > 1:16515/1‑90 (MQ=255)
ttGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGataata > 1:48726/1‑90 (MQ=255)
gcCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTa < 1:55039/39‑1 (MQ=255)
gcCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTa > 2:55039/1‑39 (MQ=39)
gcCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACa > 1:173995/1‑90 (MQ=255)
cTGAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCAtaatt > 2:184058/1‑90 (MQ=255)
gAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCAtaattaa > 2:256178/1‑90 (MQ=255)
gCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAATTAATAGAt > 1:89144/1‑90 (MQ=255)
|
AGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAATACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAATTAATAGAT > NZ_CP009273/1838736‑1838896
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ 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
GGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAATACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGA > NZ_CP009273/1838716‑1838901
|
GGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACAC < SRR3722091.54486/100‑1 (MQ=60)
taagatggtcgtcggcagcgtcagatgtgtataagagacagCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGATAAT < SRR3722091.55863/59‑1 (MQ=60)
TTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAATACCTTTAGCCTGATAATA > SRR3722091.16751/1‑100 (MQ=60)
TTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGATAATA > SRR3722091.49430/1‑100 (MQ=60)
ATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTG < SRR3722091.175369/100‑1 (MQ=60)
CAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACA > SRR3722091.176488/1‑100 (MQ=60)
ATACTGAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAATTAATAGAT > SRR3722091.90526/1‑100 (MQ=60)
GAGCGGCTCAACACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAActgtctcttataca > SRR3722091.204971/1‑86 (MQ=60)
|
GGCCAGAACTCCAGCGTTTCAGACCATTGCTTATGATAGAATCGATAAGGTGCGATCAATTGTAGCGCCTGTAACTTCTCGATACTGAGCGGCTCAATACCTTTAGCCTGATAATAATGCAGTTGTTCTTTTTTTGCTTTAAAACCGGCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGA > NZ_CP009273/1838716‑1838901
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |