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 |
A23 F58 I0 R1
|
19 |
25.6 |
590480 |
94.5% |
558003 |
140.9 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
RA |
USA300TCH1516_ALE |
2,010,363 |
C→T |
100% |
C146Y (TGT→TAT) |
USA300HOU_RS10125 ← |
Putative multidrug export ATP‑binding/permease protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | USA300TCH1516_ALE | 2,010,363 | 0 | C | T | 96.6%
| 45.7
/ ‑6.0
| 19 | C146Y (TGT→TAT) | USA300HOU_RS10125 | Putative multidrug export ATP‑binding/permease protein |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (11/7); minor base G (1/0); total (12/7) |
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 |
Rejected as polymorphism: E-value score below prediction cutoff. |
Rejected as polymorphism: Frequency below/above cutoff threshold. |
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCT > USA300TCH1516_ALE/2010227‑2010485
|
ttCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACAGTTAAAATGTAAAATGAGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACAGTATGGATAGTGCAATAATAATTGTTATATAATc > 1:112535/1‑141 (MQ=255)
tCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCaaa > 1:38526/1‑141 (MQ=255)
tCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGCTAACATGTAAAATGGGAAGAAAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACAATATGAATAGTGCAATACTAATTCTTCTAGCATCTACCCaaa > 2:27997/1‑141 (MQ=255)
cAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAat < 2:90833/141‑1 (MQ=255)
ttttCTTAATCGGCACAAGCAAACGTACACCGTTAAAAGGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATAtt < 2:30034/141‑1 (MQ=255)
tCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCAt > 1:232867/1‑141 (MQ=255)
ttAATCTTCCAAAGAAAACGTACACCGTTAATATGTAAAATGGGAAGATAAACAGTGCTGCTAAAATCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAACAAATCTTCATTa > 2:24160/1‑141 (MQ=255)
aaaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg > 1:187563/1‑141 (MQ=255)
aaaatgtaaaaggGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTttgtttgt < 2:130643/141‑1 (MQ=255)
cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa > 2:293593/1‑141 (MQ=255)
tgctgcTAAAGCCAAGTTCACAGCTAAAAAGAACATTATTGATAGTGCAATAATAATTGTTATATAATCTAACCATATATTCAGTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat < 1:130962/141‑1 (MQ=255)
tgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATTTTTTGTTAGTACAACATCATTAATCACTCTAGATATTAc > 1:111259/1‑141 (MQ=255)
aaTTTCATATATAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt < 1:131009/141‑1 (MQ=255)
cacaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGatt > 2:135840/1‑141 (MQ=255)
tAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGc > 1:184842/1‑141 (MQ=255)
aaaaaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCa < 1:77323/141‑1 (MQ=255)
aaaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATa < 1:86545/141‑1 (MQ=255)
tAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACACGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTaa > 2:58043/1‑141 (MQ=255)
gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTTGCACTTAAAGCt > 2:289066/1‑141 (MQ=255)
|
TTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCT > USA300TCH1516_ALE/2010227‑2010485
|
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