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 I200 R1
|
216 |
30.8 |
1683622 |
97.2% |
1636480 |
87.0 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
4,367,563 |
T→C |
intergenic (‑24/+65) |
blc ← / ← ampC |
lipocalin Blc/BlaEC family class C beta‑lactamase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 4,367,563 | 0 | T | C | 84.0%
| 51.6
/ 4.8
| 25 | intergenic (‑24/+65) | blc/ampC | lipocalin Blc/BlaEC family class C beta‑lactamase |
| Reads supporting (aligned to +/- strand): ref base T (1/3); new base C (9/12); total (10/15) |
| Fisher's exact test for biased strand distribution p-value = 6.26e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 |
AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA > NZ_CP009273/4367477‑4367647
|
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTa < 2:280909/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa < 2:505515/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa < 2:491389/90‑1 (MQ=255)
cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg > 2:562465/1‑89 (MQ=255)
cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg < 1:562465/89‑1 (MQ=255)
cGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAg < 2:179832/90‑1 (MQ=255)
cGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAg < 2:231461/90‑1 (MQ=255)
ccAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCaa < 1:650495/90‑1 (MQ=255)
cAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCg > 2:562817/1‑90 (MQ=255)
tGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGgaga > 2:18261/1‑90 (MQ=255)
ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc < 1:38538/90‑1 (MQ=255)
ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc > 2:104153/1‑90 (MQ=255)
cAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa > 1:397582/1‑90 (MQ=255)
aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt < 1:780602/90‑1 (MQ=255)
cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGCAAATTCGGACCCGATGGAAt < 1:526418/90‑1 (MQ=255)
cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt < 1:104153/90‑1 (MQ=255)
aTAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt > 2:248303/1‑90 (MQ=255)
aTAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt > 2:580817/1‑90 (MQ=255)
tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt < 1:102909/90‑1 (MQ=255)
ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa < 1:248303/90‑1 (MQ=255)
tCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag > 2:625568/1‑90 (MQ=255)
cTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag < 1:832123/90‑1 (MQ=255)
cTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag > 2:470564/1‑90 (MQ=255)
cTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag < 2:761959/90‑1 (MQ=255)
tACTGGTTTTTTCCTAAGTGTAGATGAAAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCg > 2:702809/1‑90 (MQ=255)
ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGaa < 1:331407/90‑1 (MQ=255)
ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGaa < 2:710857/90‑1 (MQ=255)
|
AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA > NZ_CP009273/4367477‑4367647
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
GGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGT > NZ_CP009273/4367475‑4367662
|
GGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGA < SRR3722079.159667/100‑1 (MQ=60)
GCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGA < SRR3722079.567938/100‑1 (MQ=60)
CCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGG < SRR3722079.657211/100‑1 (MQ=60)
CGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGA > SRR3722079.400915/1‑100 (MQ=60)
GGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAT < SRR3722079.38824/100‑1 (MQ=60)
AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC < SRR3722079.788968/100‑1 (MQ=60)
CATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGCAAATTCGGACCCGATGGAATTTTACTGTAG < SRR3722079.531312/100‑1 (MQ=60)
CATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAG < SRR3722079.104986/100‑1 (MQ=60)
TAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAG < SRR3722079.103734/100‑1 (MQ=60)
TTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAG < SRR3722079.250384/100‑1 (MQ=60)
CTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAT < SRR3722079.841030/100‑1 (MQ=60)
TTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGC < SRR3722079.334048/100‑1 (MQ=60)
TTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG > SRR3722079.100745/1‑100 (MQ=60)
CCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGT > SRR3722079.749226/1‑100 (MQ=60)
|
GGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGT > NZ_CP009273/4367475‑4367662
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |