| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 599,499 | C→T | 46.1% | R614R (CGC→CGT) | dmsA → | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 599,499 | 0 | C | T | 46.1% | 2.6 / 29.4 | 39 | R614R (CGC→CGT) | dmsA | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Reads supporting (aligned to +/- strand): ref base C (9/12); new base T (8/10); total (17/22) | |||||||||||
| 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.63e-01 | |||||||||||
GCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCC > minE/599439‑599560 | gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:2257027/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:378891/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:1204866/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:2238376/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:2638443/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:3219691/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAggg > 1:1460114/1‑69 (MQ=255)gCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGTAAGCAgg > 1:1450067/1‑68 (MQ=255) tGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATc > 1:2698705/1‑68 (MQ=255) tGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATc > 1:2161675/1‑68 (MQ=255) cgcgGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAgcgc > 1:3006312/1‑68 (MQ=255) cgcgGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGcgcg < 1:1333123/69‑1 (MQ=255) gcgGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCa > 1:2551824/1‑48 (MQ=255) gcgGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGa > 1:2751601/1‑69 (MQ=255) gcgGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGa > 1:1547031/1‑69 (MQ=255) gcgGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGAGCTTTAAAAAGcgcg < 1:2604234/68‑1 (MQ=255) tGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTg < 1:3410641/69‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGGAGGGGATCTTTaaa < 1:2044843/37‑1 (MQ=38) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:2915857/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:3388502/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:430774/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:881973/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:2263757/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:1302850/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:1411271/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:2042470/37‑1 (MQ=255) aaCGTTTGAAGAGTTCCGTAAGCAGGGGATCTTTaaa < 1:2271088/37‑1 (MQ=255) tttGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATaaa < 1:2127859/69‑1 (MQ=255) tttGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATaaa < 1:2063997/69‑1 (MQ=255) tttGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATaaa < 1:2219169/69‑1 (MQ=255) gagTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACTTTGCTTATAAAGCCTTc < 1:1034993/69‑1 (MQ=255) gagTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATaaa < 1:834680/63‑1 (MQ=255) gagTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc < 1:2206661/69‑1 (MQ=255) gagTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc < 1:2854030/69‑1 (MQ=255) gagTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc < 1:1167876/69‑1 (MQ=255) agTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc < 1:1113340/68‑1 (MQ=255) gTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTcc > 1:2677660/1‑68 (MQ=255) gTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTcc > 1:337429/1‑68 (MQ=255) gTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTcc > 1:1258369/1‑68 (MQ=255) | GCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCC > minE/599439‑599560 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |