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,671,899 |
T→G |
T88P (ACG→CCG) |
pntA ← |
Re/Si‑specific NAD(P)(+) transhydrogenase subunit alpha |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,671,899 | 0 | T | G | 90.0%
| 86.7
/ 2.7
| 30 | T88P (ACG→CCG) | pntA | Re/Si‑specific NAD(P)(+) transhydrogenase subunit alpha |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base G (13/14); total (14/16) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01 |
CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC > NZ_CP009273/1671816‑1671981
|
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc > 2:672183/1‑90 (MQ=255)
aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAGTTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCACAGGATTc > 2:198096/1‑90 (MQ=255)
aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGTGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc < 1:758312/90‑1 (MQ=255)
aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc > 1:28665/1‑90 (MQ=255)
cgttacgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCaa < 2:66373/90‑1 (MQ=255)
ttacgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAAtt > 1:314509/1‑90 (MQ=255)
gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTtcat < 1:437078/90‑1 (MQ=255)
gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat < 2:611882/90‑1 (MQ=255)
gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat > 2:14790/1‑90 (MQ=255)
gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat < 2:13108/90‑1 (MQ=255)
ttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcatc < 1:14790/90‑1 (MQ=255)
cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCatcatc > 2:488264/1‑90 (MQ=255)
tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg < 1:729047/90‑1 (MQ=255)
tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg < 1:443148/90‑1 (MQ=255)
cATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCAtt < 1:138315/90‑1 (MQ=255)
aTTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTg < 1:89974/90‑1 (MQ=255)
ttAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGa > 2:570175/1‑90 (MQ=255)
tAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAc > 1:573631/1‑90 (MQ=255)
aaTTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAcc > 2:814096/1‑90 (MQ=255)
ttCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCtt < 2:60170/90‑1 (MQ=255)
cTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAt > 1:755313/1‑90 (MQ=255)
cTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAt > 2:780980/1‑90 (MQ=255)
ggCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc < 1:563366/90‑1 (MQ=255)
ggCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc < 2:755313/90‑1 (MQ=255)
gagaaaaCTCACCAGCGGTGTCCCAGGATTCAGTACCGCAATTTCATCATCTACCGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa < 2:522458/87‑1 (MQ=255)
gATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa < 2:682072/90‑1 (MQ=255)
aTAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAg < 1:420090/90‑1 (MQ=255)
aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc > 2:556363/1‑90 (MQ=255)
aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc > 1:731588/1‑90 (MQ=255)
aCCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTc > 1:297183/1‑90 (MQ=255)
|
CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC > NZ_CP009273/1671816‑1671981
|
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
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTTCTACAATTTCAGC > NZ_CP009273/1671801‑1671998
|
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGT > SRR3722079.254653/1‑100 (MQ=60)
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGT > SRR3722079.354428/1‑100 (MQ=60)
AGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGA < SRR3722079.656550/100‑1 (MQ=60)
GTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTC > SRR3722079.28878/1‑100 (MQ=60)
ATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGTGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAA < SRR3722079.766425/100‑1 (MQ=60)
CACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATT > SRR3722079.317050/1‑100 (MQ=60)
GTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGG < SRR3722079.440871/100‑1 (MQ=60)
TTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGC < SRR3722079.14896/100‑1 (MQ=60)
GTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAC > SRR3722079.579282/1‑100 (MQ=60)
TTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACC < SRR3722079.447014/100‑1 (MQ=60)
TTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACC < SRR3722079.736839/100‑1 (MQ=60)
CATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGA < SRR3722079.139486/100‑1 (MQ=60)
ATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAA < SRR3722079.90686/100‑1 (MQ=60)
ATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAT > SRR3722079.763388/1‑100 (MQ=60)
GGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGC < SRR3722079.568852/100‑1 (MQ=60)
GGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGC > SRR3722079.739417/1‑100 (MQ=60)
GATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC > SRR3722079.299619/1‑100 (MQ=60)
ATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCC < SRR3722079.423708/100‑1 (MQ=60)
GTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTTCTACAATTTCAGC > SRR3722079.504313/1‑100 (MQ=60)
|
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTTCTACAATTTCAGC > NZ_CP009273/1671801‑1671998
|
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: ‑ |