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 I197 R1
|
189 |
12.7 |
697584 |
97.1% |
677354 |
86.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 | 83.3%
| 21.9
/ ‑1.6
| 12 | intergenic (‑280/‑379) | wza/yegH | polysaccharide export protein Wza/TerC family protein |
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (4/6); total (4/8) |
Fisher's exact test for biased strand distribution p-value = 5.15e-01 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCT > NZ_CP009273/2130917‑2131083
|
tGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGgaa > 2:269667/1‑90 (MQ=255)
cAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAatat < 2:179497/90‑1 (MQ=255)
tGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGAc < 1:162769/90‑1 (MQ=255)
gATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACa > 1:335532/1‑90 (MQ=255)
gTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATgg > 2:325386/1‑90 (MQ=255)
tCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAata > 1:257722/1‑90 (MQ=255)
gCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCaat < 2:7846/90‑1 (MQ=255)
cAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGc > 2:71796/1‑90 (MQ=255)
aGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCa < 2:305235/90‑1 (MQ=255)
cGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCagga < 1:41473/90‑1 (MQ=255)
tgatgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAggagga < 2:173663/90‑1 (MQ=255)
atgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAg < 1:269667/90‑1 (MQ=255)
gTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCt < 1:313008/90‑1 (MQ=255)
|
TGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCT > NZ_CP009273/2130917‑2131083
|
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
CCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGA > NZ_CP009273/2130913‑2131102
|
CCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACT < SRR3722076.285522/100‑1 (MQ=60)
CCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACA > SRR3722076.340628/1‑100 (MQ=60)
TGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAA < SRR3722076.164748/100‑1 (MQ=60)
GATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACctgtctcttatacacatctgacgctgccg > SRR3722076.250045/1‑71 (MQ=60)
GTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATA > SRR3722076.261437/1‑100 (MQ=60)
CGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTA < SRR3722076.41963/100‑1 (MQ=60)
ATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCT < SRR3722076.273614/100‑1 (MQ=60)
GTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAG < SRR3722076.317751/100‑1 (MQ=60)
GGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGA > SRR3722076.113539/1‑100 (MQ=60)
|
CCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGA > NZ_CP009273/2130913‑2131102
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |