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 |
A9 F110 I0 R1
|
28 |
36.3 |
769303 |
98.8% |
760071 |
141.0 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
RA |
USA300TCH1516_ALE |
2,028,269 |
G→C |
100% |
P126A (CCA→GCA) |
USA300HOU_RS10230 ← |
hypothetical protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | USA300TCH1516_ALE | 2,028,269 | 0 | G | C | 96.7%
| 95.2
/ ‑2.8
| 30 | P126A (CCA→GCA) | USA300HOU_RS10230 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base C (11/18); total (11/19) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01 |
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. |
CCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAAC > USA300TCH1516_ALE/2028134‑2028402
|
ccTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATATACTGCTGTAc < 1:554687/141‑1 (MQ=255)
cACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTg < 1:389576/141‑1 (MQ=255)
aaaTGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAAcc < 1:655435/141‑1 (MQ=255)
aaTGTGTCTAACAACAGTGGTATTATTTTCCTTAATATCTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTAATATACATAAACTGCTGTACGTTGCTCACCAAACCa > 1:424274/1‑141 (MQ=255)
caacaaTGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAACTAACTATCTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTtaa < 1:750029/141‑1 (MQ=255)
ggTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGACAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGa < 1:64715/141‑1 (MQ=255)
gTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGAt < 1:490574/141‑1 (MQ=255)
aGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGt > 1:297829/1‑141 (MQ=255)
gCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTc > 1:57908/1‑141 (MQ=255)
ttttGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTGTATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAGCCTCAGCcgccgc > 1:254560/1‑141 (MQ=255)
aaGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactact < 1:102675/141‑1 (MQ=255)
gTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACttt > 1:667623/1‑141 (MQ=255)
tCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtttt < 1:764992/141‑1 (MQ=255)
tatGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTg > 1:308281/1‑141 (MQ=255)
aGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGCCGCAACTACTTTTTTTGGTTTa > 1:595595/1‑141 (MQ=255)
cAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCACCTACTTTTTTTGGTTTAAACTtgt > 1:44771/1‑141 (MQ=255)
gTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAAt < 1:687711/141‑1 (MQ=255)
aTTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATTGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGc < 1:479999/141‑1 (MQ=255)
ttGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGcac < 1:767290/141‑1 (MQ=255)
tACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAAt > 1:61038/1‑141 (MQ=255)
tACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAAt < 1:50184/141‑1 (MQ=255)
tCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGt < 1:661357/141‑1 (MQ=255)
cctCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGta < 1:254838/139‑1 (MQ=255)
ttCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaa < 1:360696/141‑1 (MQ=255)
ttCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaa < 1:751062/141‑1 (MQ=255)
cACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTTTTTTTTTTTTGTGTATTACATAACGCACAATACGtaatact > 1:27583/1‑139 (MQ=255)
ataCAGAAACTGCTGGACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACTTCATCCGTCGCAACTTCTTTTTTTGGTTTAAACTTGTGTCTTCGGTCTCTCACAATAAGTAATAATAAAAAGa < 1:182407/141‑1 (MQ=255)
taCATAAACTGCTGTACGTTGCTCACCAAACCACTTATGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATACGAAATAATAAAAAGAc < 1:423664/141‑1 (MQ=255)
cATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGa > 1:746786/1‑141 (MQ=255)
tacaCTGCTGTACGATGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAAc < 1:72338/138‑1 (MQ=255)
|
CCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAAC > USA300TCH1516_ALE/2028134‑2028402
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A