| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,687,242 | C→T | 75.1% | E126E (GAG→GAA) | sseB ← | rhodanase‑like enzyme, sulfur transfer from thiosulfate |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,687,242 | 0 | C | T | 75.1% | 38.7 / 9.5 | 36 | E126E (GAG→GAA) | sseB | rhodanase‑like enzyme, sulfur transfer from thiosulfate |
| Reads supporting (aligned to +/- strand): ref base C (5/4); new base T (16/11); total (21/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTA > minE/1687199‑1687308 | aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTGCCTCTTCACCAATCAACAAACTGAttt < 1:1244807/69‑1 (MQ=255)aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt < 1:56674/69‑1 (MQ=255)aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt < 1:3115155/69‑1 (MQ=255)aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt < 1:3263548/69‑1 (MQ=255) ccAGTACTTCCTGGCTACTCAGGGAATTCCCTTCTTCACCAATc > 1:223189/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2923367/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2926496/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2648455/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:2946256/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:3225919/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:3261280/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:1531145/1‑44 (MQ=255) ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc > 1:81742/1‑44 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:3212823/60‑1 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:3170187/60‑1 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:3367441/60‑1 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:367206/60‑1 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:860812/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:1397517/60‑1 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:1418632/60‑1 (MQ=255) aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg < 1:2907634/60‑1 (MQ=255) tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:1577977/1‑68 (MQ=255) tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:1237730/1‑68 (MQ=255) tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:455036/1‑68 (MQ=255) tttCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:2261847/1‑68 (MQ=255) tttCCCTCTTCAACAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc > 1:2348793/1‑68 (MQ=255) ccTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa < 1:1443239/68‑1 (MQ=255) ccTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa < 1:1637306/68‑1 (MQ=255) ccTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa < 1:648057/68‑1 (MQ=255) ccTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa < 1:794224/68‑1 (MQ=255) ccTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa < 1:2726768/68‑1 (MQ=255) ctctTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa < 1:1486231/67‑1 (MQ=255) | ATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTA > minE/1687199‑1687308 |
| 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: ‑ |