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 |
A15 F54 I0 R1
|
17 |
21.8 |
494604 |
94.9% |
469379 |
140.9 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
RA |
USA300TCH1516_ALE |
1,744,428 |
T→C |
37.0% |
K82E (AAA→GAA) |
apt ← |
Adenine phosphoribosyltransferase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | USA300TCH1516_ALE | 1,744,428 | 0 | T | C | 37.0%
| 6.3
/ 15.3
| 19 | K82E (AAA→GAA) | apt | Adenine phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base T (6/6); new base C (4/3); total (10/9) |
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.49e-01 |
CAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGT > USA300TCH1516_ALE/1744302‑1744563
|
caGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGACGTCATTTCCCTTCTTTACtca > 2:33470/1‑139 (MQ=255)
taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGc < 1:9134/141‑1 (MQ=255)
taataaATCATCTGTAATTAACACACGTTGACCGGGTGTAAGTGCATGTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGc < 1:232294/141‑1 (MQ=255)
aaTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTTCAAAGCCAATCCCCAt > 1:184095/1‑141 (MQ=255)
aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATACCCATTGaa > 1:66407/1‑141 (MQ=255)
ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTCCAGGa > 1:9753/1‑141 (MQ=255)
cATCTTTGTGCATTGTTAAAACATTTGTATCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCACAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa > 2:182492/1‑141 (MQ=255)
aCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATACAGCCACGCGCTTCAGGTCCt > 1:87641/1‑141 (MQ=255)
aCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt < 2:186782/141‑1 (MQ=255)
cTAGGTCATACTCATAACGAATGACTTCACGAGGTACTTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACACGACAGCCAACGATAAAGCCACGCCCTTCACATCCTACAACAATACCAAcct > 2:58590/1‑139 (MQ=255)
cTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt < 1:147417/141‑1 (MQ=255)
cATAACGAATGACTTCACGAGGTAATTTCCCTACTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGTCAGCCAATGATAAAGCCACGCGCTGCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg < 1:30918/141‑1 (MQ=255)
tAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCg > 1:242204/1‑141 (MQ=255)
taCTTCATGAGGTAATTCTCTTTCTTTTCGAACATGTTCAAAGCCAATCCCAATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc < 2:154534/140‑1 (MQ=255)
ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACAATATCAACATCTCTGTCTTTTGCGTCTTCTACAAt > 2:144535/1‑141 (MQ=255)
cGAGGTAATTTACCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACCGGACAGACAATGATAAAGCCACGCGCTTCAGGTCCTACAACTATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTa < 1:95647/141‑1 (MQ=255)
tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt < 1:52293/141‑1 (MQ=255)
tAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt < 1:52286/141‑1 (MQ=255)
tAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGt > 1:192930/1‑141 (MQ=255)
|
CAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGT > USA300TCH1516_ALE/1744302‑1744563
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A