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,259,173 |
A→G |
T2A (ACC→GCC) |
hemA → |
glutamyl‑tRNA reductase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | NZ_CP009273 | 1,259,173 | 0 | A | G | 95.8%
| 69.7
/ ‑3.5
| 24 | T2A (ACC→GCC) | hemA | glutamyl‑tRNA reductase |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (12/11); total (13/11) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01 |
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATA > NZ_CP009273/1259086‑1259248
|
gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGcc < 2:189031/90‑1 (MQ=255)
ttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTc < 2:330310/90‑1 (MQ=255)
tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa < 1:206748/90‑1 (MQ=255)
tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa > 2:352897/1‑90 (MQ=255)
tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCa < 1:544838/90‑1 (MQ=255)
aTCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATaaa > 2:105696/1‑90 (MQ=255)
gTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGc > 2:677304/1‑90 (MQ=255)
gatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCt > 1:227678/1‑90 (MQ=255)
gatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCt > 2:274502/1‑90 (MQ=255)
atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg > 1:836162/1‑90 (MQ=255)
aaGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATc < 1:105696/90‑1 (MQ=255)
aGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCg > 1:448480/1‑90 (MQ=255)
aGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCg > 2:34300/1‑90 (MQ=255)
gCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGc < 2:777216/90‑1 (MQ=255)
aCTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCg > 1:695300/1‑90 (MQ=255)
aaCCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga < 1:689528/90‑1 (MQ=255)
aaCCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga < 1:146698/90‑1 (MQ=255)
aCCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGaa < 2:160804/90‑1 (MQ=255)
aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt > 1:481997/1‑90 (MQ=255)
aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt > 2:434671/1‑90 (MQ=255)
aTCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt > 2:661788/1‑90 (MQ=255)
tGGTATTATTTCCCGCAGACATGGCCC‑TTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTc > 1:283398/1‑90 (MQ=255)
tattatTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGc < 1:630494/90‑1 (MQ=255)
ttCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATa < 1:571582/90‑1 (MQ=255)
|
GGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATA > NZ_CP009273/1259086‑1259248
|
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
AGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTG > NZ_CP009273/1259083‑1259266
|
AGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAG < SRR3722079.216826/100‑1 (MQ=60)
TCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAAC < SRR3722079.208525/100‑1 (MQ=60)
TCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAAC < SRR3722079.550026/100‑1 (MQ=60)
CATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCT > SRR3722079.229609/1‑100 (MQ=60)
ATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTG > SRR3722079.845111/1‑100 (MQ=60)
GTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCG > SRR3722079.452400/1‑100 (MQ=60)
ATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCG > SRR3722079.702619/1‑100 (MQ=60)
AAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAA < SRR3722079.106554/100‑1 (MQ=60)
GACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGT > SRR3722079.486210/1‑100 (MQ=60)
AACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCG < SRR3722079.147966/100‑1 (MQ=60)
AACCCTATCAACGTTGGTATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATAG < SRR3722079.696777/100‑1 (MQ=60)
CTATCAACGTTGGTATTATTTCCCGCAGACATGGCCC‑TTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTC > SRR3722079.285749/1‑100 (MQ=60)
TATTATTTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCT < SRR3722079.636939/100‑1 (MQ=60)
TTCCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCA < SRR3722079.577196/100‑1 (MQ=60)
CCCGCAGACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGG > SRR3722079.406566/1‑100 (MQ=60)
GACATGGCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTG > SRR3722079.416489/1‑100 (MQ=60)
|
AGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTG > NZ_CP009273/1259083‑1259266
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |