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 F5 I97 R1
|
275 |
44.0 |
1733530 |
96.5% |
1672856 |
122.0 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
678,913 |
A→G |
intergenic (‑64/+20) |
hscC ← / ← rihA |
molecular chaperone HscC/pyrimidine‑specific ribonucleoside hydrolase RihA |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 678,913 | 0 | A | G | 93.8%
| 95.7
/ ‑0.8
| 32 | intergenic (‑64/+20) | hscC/rihA | molecular chaperone HscC/pyrimidine‑specific ribonucleoside hydrolase RihA |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (19/11); total (20/12) |
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.01e-01 |
AGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAAC > NZ_CP009273/678833‑679043
|
aGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATcaa > 2:73173/1‑138 (MQ=255)
tCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACg < 2:763669/84‑1 (MQ=255)
tCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACg > 1:763669/1‑84 (MQ=255)
tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCTGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 2:224962/1‑139 (MQ=255)
tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 2:256259/1‑139 (MQ=255)
tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 2:18819/1‑139 (MQ=255)
tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 1:829883/1‑139 (MQ=255)
tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 1:633160/1‑139 (MQ=255)
cATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAATCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACa < 2:847354/139‑1 (MQ=255)
cATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACa < 2:98116/139‑1 (MQ=255)
ttATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACa > 2:602755/1‑136 (MQ=255)
tGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTAc > 1:90932/1‑139 (MQ=255)
gaaagaaaAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTa < 1:253780/137‑1 (MQ=255)
gaaagaaaAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTa < 2:633160/137‑1 (MQ=255)
aaGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCAtt < 1:73173/139‑1 (MQ=255)
aaaCAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTa > 2:285310/1‑139 (MQ=255)
aTTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCt > 1:299049/1‑139 (MQ=255)
cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:761864/1‑138 (MQ=255)
cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:346311/1‑138 (MQ=255)
cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 2:357924/1‑138 (MQ=255)
cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:27311/1‑138 (MQ=255)
cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:632998/1‑138 (MQ=255)
ccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataata < 2:619723/139‑1 (MQ=255)
cAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataataa < 1:285310/139‑1 (MQ=255)
cgcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGAc > 2:817291/1‑92 (MQ=255)
cgcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGAc < 1:817291/92‑1 (MQ=255)
gcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac < 1:41423/139‑1 (MQ=255)
gcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac < 2:346311/139‑1 (MQ=255)
cgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac > 2:529070/1‑138 (MQ=255)
cgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac > 2:163516/1‑138 (MQ=255)
gtgGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACTATCAAc < 1:436962/80‑1 (MQ=255)
gtgGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAAc > 2:436962/1‑80 (MQ=255)
|
AGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAAC > NZ_CP009273/678833‑679043
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
GCGGCACCGTCTTTCCAGACGGCAATTAAACTATTGGTAGTACCGAGATCGATACCAATGGCGAGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAACAACCGTCATACCCTGGGTAT > NZ_CP009273/678770‑679063
|
GCGGCACCGTCTTTCCAGACGGCAATTAAACTATTGGTAGTACCGAGATCGATACCAATGGCGAGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAAT < SRR3722120.54342/151‑1 (MQ=60)
CGGCACCGGCTTTCCAGACGGCAATTAAACTATTGGTAGTACCGAGATCGATACCAATGGCGAGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATC < SRR3722120.353093/151‑1 (MQ=60)
GTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGctgtctcttatacacatctccgagcccacgagactaaggcgaatc > SRR3722120.773959/1‑106 (MQ=60)
GCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAAC > SRR3722120.841039/1‑151 (MQ=60)
GCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAAC > SRR3722120.641906/1‑151 (MQ=60)
acgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACTATCAACATCGACCATTA < SRR3722120.443349/102‑1 (MQ=60)
CATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACG > SRR3722120.92102/1‑151 (MQ=60)
tagtctcgtgggctcggagatgtgtataagagacagCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGG < SRR3722120.828296/114‑1 (MQ=60)
TGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTC < SRR3722120.257456/149‑1 (MQ=60)
AAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGC < SRR3722120.74137/151‑1 (MQ=60)
GAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTG > SRR3722120.303337/1‑151 (MQ=60)
GATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAAT > SRR3722120.351321/1‑150 (MQ=60)
GATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAAT > SRR3722120.641743/1‑150 (MQ=60)
GATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAAT > SRR3722120.772135/1‑150 (MQ=60)
GATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAAT > SRR3722120.27619/1‑150 (MQ=60)
GATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGCGTGTTAAGCGCAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCACTACGGTGGCATTCGGTTTAGTGCCTGTCAGATAAT > SRR3722120.240633/1‑150 (MQ=60)
gtgtataagagacagAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAAT < SRR3722120.38136/135‑1 (MQ=60)
CCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAACAAC < SRR3722120.289415/151‑1 (MQ=60)
CGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAACAACCGTCATAC < SRR3722120.41924/151‑1 (MQ=60)
GTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAACAACCGTCATACCCTGGGTAT > SRR3722120.331495/1‑151 (MQ=60)
|
GCGGCACCGTCTTTCCAGACGGCAATTAAACTATTGGTAGTACCGAGATCGATACCAATGGCGAGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAACAACCGTCATACCCTGGGTAT > NZ_CP009273/678770‑679063
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |