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 |
|
A2 F1 I0 R1
|
4577 |
57.4 |
3436498 |
77.7% |
2670158 |
63.4 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
| RA |
minE |
1,687,239 |
T→C |
100% |
G127G (GGA→GGG) |
sseB ← |
rhodanase‑like enzyme, sulfur transfer from thiosulfate |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | minE | 1,687,239 | 0 | T | C | 87.3%
| 53.8
/ ‑0.3
| 31 | G127G (GGA→GGG) | sseB | rhodanase‑like enzyme, sulfur transfer from thiosulfate |
| Reads supporting (aligned to +/- strand): ref base T (0/4); new base C (16/11); total (16/15) |
| Fisher's exact test for biased strand distribution p-value = 4.34e-02 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01 |
| Rejected as polymorphism: E-value score below prediction cutoff. |
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
ATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGC > minE/1687199‑1687304
|
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTGCCTCTTCACCAATCAACAAACTGAttt < 1:1244807/69‑1 (MQ=255)
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt < 1:3115155/69‑1 (MQ=255)
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt < 1:3263548/69‑1 (MQ=255)
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt < 1:56674/69‑1 (MQ=255)
ccAGTACTTCCTGGCTACTCAGGGAATTCCCTTCTTCACCAATc > 1:223189/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2923367/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2648455/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2926496/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:3225919/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:3261280/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:81742/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2946256/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2255422/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2159651/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2103404/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:1980290/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:1940686/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:185022/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:1686313/1‑44 (MQ=255)
ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:1531145/1‑44 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:860812/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:3170187/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:3212823/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:3367441/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:367206/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:2907634/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:2636104/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:2633094/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:2450489/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:1418632/60‑1 (MQ=255)
aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:1397517/60‑1 (MQ=255)
tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:2261847/1‑68 (MQ=255)
tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:1577977/1‑68 (MQ=255)
tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:455036/1‑68 (MQ=255)
tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:1237730/1‑68 (MQ=255)
tttCCCTCTTCAACAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:2348793/1‑68 (MQ=255)
|
ATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGC > minE/1687199‑1687304
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A