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 I228 R1
|
214 |
26.5 |
1448192 |
97.1% |
1406194 |
86.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
2,102,646 |
T→C |
G140G (GGA→GGG) |
wzx ← |
O16 family O‑antigen flippase |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,102,646 | 0 | T | C | 96.8%
| 95.4
/ ‑3.4
| 31 | G140G (GGA→GGG) | wzx | O16 family O‑antigen flippase |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (14/16); total (15/16) |
| Fisher's exact test for biased strand distribution p-value = 4.84e-01 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01 |
TTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTA > NZ_CP009273/2102561‑2102727
|
ttGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCt > 1:521684/1‑90 (MQ=255)
gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa > 1:163874/1‑90 (MQ=255)
tgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGaa < 1:190059/90‑1 (MQ=255)
ttCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAAc < 1:712829/90‑1 (MQ=255)
ttCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAAc < 2:163874/90‑1 (MQ=255)
acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc < 2:207663/90‑1 (MQ=255)
aatGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAggg > 1:333607/1‑90 (MQ=255)
atGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGt < 2:33782/90‑1 (MQ=255)
agGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTaa < 1:722686/90‑1 (MQ=255)
gATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAAtt > 1:173025/1‑90 (MQ=255)
ccAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTa > 2:124492/1‑90 (MQ=255)
cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa > 1:446291/1‑90 (MQ=255)
cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa > 1:449490/1‑90 (MQ=255)
cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa > 2:292914/1‑90 (MQ=255)
cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa > 2:549280/1‑90 (MQ=255)
aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAAt > 1:63812/1‑90 (MQ=255)
gCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGc < 1:331312/90‑1 (MQ=255)
aTACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAAc < 2:278702/90‑1 (MQ=255)
tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg < 1:1351/90‑1 (MQ=255)
tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg < 2:417156/90‑1 (MQ=255)
tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg < 1:124492/90‑1 (MQ=255)
gctTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCa > 1:172799/2‑90 (MQ=255)
aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg > 1:457287/1‑90 (MQ=255)
aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg > 2:365809/1‑90 (MQ=255)
gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg > 1:540426/1‑90 (MQ=255)
gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg < 2:449490/90‑1 (MQ=255)
gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg < 2:691473/90‑1 (MQ=255)
ttcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGc < 2:716271/90‑1 (MQ=255)
tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt < 1:317182/90‑1 (MQ=255)
tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt < 1:176428/90‑1 (MQ=255)
atttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTg > 2:267976/1‑90 (MQ=255)
ttttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTta < 2:157064/90‑1 (MQ=255)
ttttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTta < 1:479469/90‑1 (MQ=255)
|
TTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTA > NZ_CP009273/2102561‑2102727
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
ATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAA > NZ_CP009273/2102551‑2102737
|
ATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCT > SRR3722112.529355/1‑100 (MQ=60)
ATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGA > SRR3722112.165893/1‑100 (MQ=60)
TGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGAT < SRR3722112.192389/100‑1 (MQ=60)
TTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAG < SRR3722112.723970/100‑1 (MQ=60)
CACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGG > SRR3722112.337655/1‑100 (MQ=60)
ATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATT > SRR3722112.175149/1‑100 (MQ=60)
AGGGATAATAGCCAAGCGGGATAAAATACTCCTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACT < SRR3722112.733966/100‑1 (MQ=60)
GGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAA > SRR3722112.452502/1‑100 (MQ=60)
GGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAA > SRR3722112.455762/1‑100 (MQ=60)
GATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAAT > SRR3722112.64584/1‑100 (MQ=60)
GCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAAC < SRR3722112.335310/100‑1 (MQ=60)
ATAAAATACTGCTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCA > SRR3722112.174923/1‑100 (MQ=60)
TACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAAT < SRR3722112.1367/100‑1 (MQ=60)
TACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAAT < SRR3722112.126032/100‑1 (MQ=60)
ACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACG > SRR3722112.463723/1‑100 (MQ=60)
CTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGG > SRR3722112.548507/1‑100 (MQ=60)
TCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC < SRR3722112.320863/100‑1 (MQ=60)
TCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC < SRR3722112.178592/100‑1 (MQ=60)
CATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCA > SRR3722112.401735/1‑100 (MQ=60)
TTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAA < SRR3722112.486338/100‑1 (MQ=60)
|
ATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAA > NZ_CP009273/2102551‑2102737
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |