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 |
1,897,230 |
A→G |
D309G (GAT→GGT) |
manX → |
PTS mannose transporter subunit IIAB |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 1,897,230 | 0 | A | G | 89.7%
| 95.8
/ 5.6
| 39 | D309G (GAT→GGT) | manX | PTS mannose transporter subunit IIAB |
| Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (15/20); total (17/22) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01 |
TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC > NZ_CP009273/1897146‑1897316
|
tttCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGa < 1:66137/90‑1 (MQ=255)
aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat < 1:338433/90‑1 (MQ=255)
aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat < 2:551420/90‑1 (MQ=255)
aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat > 1:542363/1‑90 (MQ=255)
atatCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATc < 2:5852/90‑1 (MQ=255)
atatCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATc > 1:532970/1‑90 (MQ=255)
gCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGc > 1:428004/1‑90 (MQ=255)
gTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaa > 2:729790/1‑90 (MQ=255)
gAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGa < 1:265451/90‑1 (MQ=255)
aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa < 2:430834/90‑1 (MQ=255)
aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa < 1:97609/90‑1 (MQ=255)
gAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTa < 1:380413/90‑1 (MQ=255)
aaTGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTaa < 2:163741/90‑1 (MQ=255)
aaTGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTaa < 1:839343/90‑1 (MQ=255)
gcgcgGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTAt < 2:330189/90‑1 (MQ=255)
gTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTgtgt < 2:122379/90‑1 (MQ=255)
gTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTgtgt < 2:139102/90‑1 (MQ=255)
cTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATc > 2:267128/1‑90 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat > 1:202908/1‑42 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat < 2:202908/42‑1 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAg > 2:442043/1‑67 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAg < 1:442043/67‑1 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc > 2:464051/1‑90 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc > 2:773180/1‑90 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc > 1:89079/1‑90 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc > 1:836064/1‑90 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc > 1:648614/1‑90 (MQ=255)
ggAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc > 1:679042/1‑90 (MQ=255)
aGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCa < 1:184755/90‑1 (MQ=255)
tCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGt > 1:652144/1‑90 (MQ=255)
tCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGt > 1:457248/1‑90 (MQ=255)
aaGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTca < 1:436535/90‑1 (MQ=255)
aGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTcac < 2:648614/90‑1 (MQ=255)
aGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTcac < 1:779055/90‑1 (MQ=255)
ttCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACtt > 2:296295/1‑90 (MQ=255)
ttCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACtt > 1:392319/1‑90 (MQ=255)
tCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTa < 2:552627/90‑1 (MQ=255)
ccGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTc < 2:392319/90‑1 (MQ=255)
ccGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTc < 2:554271/90‑1 (MQ=255)
|
TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC > NZ_CP009273/1897146‑1897316
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTCTTA > NZ_CP009273/1897146‑1897319
|
TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAAT < SRR3722079.66676/100‑1 (MQ=60)
CGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGAT > SRR3722079.547505/1‑100 (MQ=60)
GATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATC > SRR3722079.537951/1‑100 (MQ=60)
AAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATC < SRR3722079.341108/100‑1 (MQ=60)
AGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGC > SRR3722079.431708/1‑100 (MQ=60)
GAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGT < SRR3722079.267676/100‑1 (MQ=60)
AACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTAT < SRR3722079.98383/100‑1 (MQ=60)
gtataagagacaGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTG < SRR3722079.445899/88‑1 (MQ=60)
GAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTG < SRR3722079.383510/100‑1 (MQ=60)
AATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGT < SRR3722079.848324/100‑1 (MQ=60)
GTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCctgtctcttatacacatctccgagcccacgagaccaga > SRR3722079.204663/1‑62 (MQ=60)
GTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAC > SRR3722079.655314/1‑100 (MQ=60)
GTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAC > SRR3722079.845012/1‑100 (MQ=60)
GTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAC > SRR3722079.89785/1‑100 (MQ=60)
GTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAC > SRR3722079.686168/1‑100 (MQ=60)
GAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGT > SRR3722079.461251/1‑100 (MQ=60)
GAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGT > SRR3722079.658895/1‑100 (MQ=60)
AGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACAC < SRR3722079.186371/100‑1 (MQ=60)
TCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTT > SRR3722079.395584/1‑100 (MQ=60)
AAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC < SRR3722079.440321/100‑1 (MQ=60)
AGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTCT < SRR3722079.787404/100‑1 (MQ=60)
GTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTCTTA > SRR3722079.571487/1‑100 (MQ=60)
|
TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTCTTA > NZ_CP009273/1897146‑1897319
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |