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 I209 R1
|
216 |
14.4 |
815930 |
96.1% |
784108 |
85.2 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
NZ_CP009273 |
2,131,004 |
A→G |
intergenic (‑280/‑379) |
wza ← / → yegH |
polysaccharide export protein Wza/TerC family protein |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 2,131,004 | 0 | A | G | 100.0%
| 32.3
/ NA
| 11 | intergenic (‑280/‑379) | wza/yegH | polysaccharide export protein Wza/TerC family protein |
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (3/8); total (3/8) |
TGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGC > NZ_CP009273/2130917‑2131082
|
tGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGgaa > 2:242134/1‑90 (MQ=255)
aGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTcc < 1:281738/90‑1 (MQ=255)
aTACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATaa < 2:262434/90‑1 (MQ=255)
cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt < 2:154760/90‑1 (MQ=255)
cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt < 2:52830/90‑1 (MQ=255)
gCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCaat < 1:242134/90‑1 (MQ=255)
aTGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTAt < 1:142451/63‑1 (MQ=255)
aTGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTAt > 2:142451/1‑63 (MQ=255)
atgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAg < 2:19809/90‑1 (MQ=255)
atgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAg < 2:74642/90‑1 (MQ=255)
aaGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTAc > 2:80939/1‑90 (MQ=255)
ggTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGc > 2:170545/1‑90 (MQ=255)
|
TGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGC > NZ_CP009273/2130917‑2131082
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
AGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAG > NZ_CP009273/2130926‑2131069
|
AGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGG < SRR3722091.286422/100‑1 (MQ=60)
atgtgtataagagacagCCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAA < SRR3722091.144557/83‑1 (MQ=60)
CTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATnnnnnnnnnnnnnnnnnnnn > SRR3722091.268467/1‑80 (MQ=60)
GCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAG < SRR3722091.245883/100‑1 (MQ=60)
|
AGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAG > NZ_CP009273/2130926‑2131069
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |