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 |
A72 F1 I1 R1
|
16 |
60.6 |
3810144 |
97.0% |
3695839 |
75.8 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
RA |
REL606 |
4,161,105 |
G→A |
E84K (GAA→AAA) |
rpoB → |
DNA‑directed RNA polymerase subunit beta |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | REL606 | 4,161,105 | 0 | G | A | 97.7%
| 163.9
/ ‑3.9
| 45 | E84K (GAA→AAA) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (20/23); total (21/24) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01 |
TCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGC > REL606/4161032‑4161176
|
tCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGaaa > 1:1052425/1‑76 (MQ=255)
tCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGaaa < 1:40210/76‑1 (MQ=255)
agagCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATg > 1:222732/1‑76 (MQ=255)
agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc > 2:368864/1‑76 (MQ=255)
agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc < 1:713823/76‑1 (MQ=255)
agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc > 1:1807076/1‑76 (MQ=255)
tACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGCAATGTCaaa > 1:1282670/1‑76 (MQ=255)
tACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCaaa < 1:613074/76‑1 (MQ=255)
aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt < 2:925780/76‑1 (MQ=255)
aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt > 2:1690714/1‑76 (MQ=255)
aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt < 2:825590/76‑1 (MQ=255)
tAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGc > 2:524633/1‑76 (MQ=255)
aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCg > 2:75262/1‑76 (MQ=255)
aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAACCCGTGGCg > 2:75245/1‑76 (MQ=255)
ttCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTg < 2:1468207/76‑1 (MQ=255)
ttCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTg < 2:1856855/76‑1 (MQ=255)
tCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGa < 2:583249/76‑1 (MQ=255)
ccGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGAc > 1:149582/1‑76 (MQ=255)
gAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCt > 2:573768/1‑76 (MQ=255)
gAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCt > 2:956096/1‑76 (MQ=255)
aGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTa > 2:164592/1‑76 (MQ=255)
cTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTAtt > 2:1508739/1‑76 (MQ=255)
tCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCt < 2:509168/76‑1 (MQ=255)
aGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgc > 1:333893/1‑76 (MQ=255)
aGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgc > 2:1469339/1‑76 (MQ=255)
ccGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGtt < 2:1117745/76‑1 (MQ=255)
gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa > 1:1201179/1‑76 (MQ=255)
gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa < 2:563303/76‑1 (MQ=255)
gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa > 2:376307/1‑76 (MQ=255)
ttGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACt < 1:1000082/76‑1 (MQ=255)
gCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCg > 1:124138/1‑76 (MQ=255)
aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt < 1:791103/76‑1 (MQ=255)
aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt > 1:1497435/1‑76 (MQ=255)
aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCATGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt > 1:1345740/1‑76 (MQ=255)
ccGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTgg < 2:1151066/76‑1 (MQ=255)
cGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGt < 2:1509742/76‑1 (MQ=255)
ggTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg < 1:402648/76‑1 (MQ=255)
ggTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg < 2:376228/76‑1 (MQ=255)
tgtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGAt < 1:1615596/76‑1 (MQ=255)
gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc < 1:376307/76‑1 (MQ=255)
gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc < 1:1071087/76‑1 (MQ=255)
gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc > 1:1131320/1‑76 (MQ=255)
tGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTAt < 2:814515/76‑1 (MQ=255)
gACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATg < 2:963818/76‑1 (MQ=255)
cGTGCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg < 1:387146/76‑1 (MQ=255)
cGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg < 1:1500552/76‑1 (MQ=255)
tCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcg < 1:111121/76‑1 (MQ=255)
ccAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcgc < 1:166771/76‑1 (MQ=255)
|
TCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGC > REL606/4161032‑4161176
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A