| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,335,701 | A→G | 60.0% | E29E (GAA→GAG) | fdoH → | formate dehydrogenase‑O, Fe‑S subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,335,701 | 0 | A | G | 60.0% | 17.0 / 30.7 | 30 | E29E (GAA→GAG) | fdoH | formate dehydrogenase‑O, Fe‑S subunit |
| Reads supporting (aligned to +/- strand): ref base A (4/8); new base G (13/5); total (17/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.10e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 | |||||||||||
CGTTCCGCGACTAACGGTCTGACCCCCGCGCCTCAGGCGCGGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCAGAGTGG > minE/2335645‑2335770 | cgtTCCGCGACTAACGGTCTGACCCCCGCGCCTCAGGCGCGGGACTTCCAGGaagaag > 1:2415139/1‑58 (MQ=255) ccccGCGCCTCAGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCt > 1:1684149/1‑71 (MQ=255) cccGCGCCTCAGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTg > 1:328051/1‑71 (MQ=255) cgcCTCAGGCGCGGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTa > 1:2312981/1‑69 (MQ=255) gcCTCAGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAg > 1:2831995/1‑71 (MQ=255) ccTCAGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGt < 1:3980880/47‑1 (MQ=255) cAGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTg > 1:2414332/1‑71 (MQ=255) cAGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTg > 1:2111511/1‑71 (MQ=255) cggCGGGGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGcc < 1:492119/67‑1 (MQ=255) aGGCGCGGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGt > 1:2320874/1‑71 (MQ=255) cgcgGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTaa > 1:904246/1‑61 (MQ=255) cgcgGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAg > 1:3589703/1‑71 (MQ=255) cgGGACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTacca < 1:1511438/43‑1 (MQ=255) cgGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGt > 1:2382119/1‑71 (MQ=255) gggACTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTg > 1:3378981/1‑71 (MQ=255) ggACTTCCAGGAAGAGGTGGCAAAACTCATCGACGTTACCACCTGt < 1:882368/46‑1 (MQ=38) ggACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGcc < 1:1789427/61‑1 (MQ=255) ggACCTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTgg < 1:3504181/71‑1 (MQ=255) aCTTCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGc > 1:3387237/1‑70 (MQ=255) aCTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGTCg < 1:159888/71‑1 (MQ=255) aCTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCg < 1:2108834/71‑1 (MQ=255) aCTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCg < 1:1868700/71‑1 (MQ=255) aCTTCCAGGAAGAAGTGGCGAAACTCACCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGGGGCg < 1:3738096/71‑1 (MQ=255) tCCAGTAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTAt > 1:3079029/1‑43 (MQ=38) tCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTAt > 1:1747736/1‑43 (MQ=255) tCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGc > 1:3703731/1‑67 (MQ=255) tCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGc > 1:126597/1‑67 (MQ=255) tCCAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCgtgt > 1:1752253/1‑71 (MQ=255) cAGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGtt < 1:683062/70‑1 (MQ=255) aGGAAGAGGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCa < 1:2447720/71‑1 (MQ=255) aGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCa < 1:2644088/71‑1 (MQ=255) gaagaaGTGGCGAAACTCATCGACGTTACCACCTGTATCg < 1:2665925/40‑1 (MQ=255) gaagaaGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCa > 1:1016890/1‑69 (MQ=255) aagaagTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCagag > 1:3308069/1‑71 (MQ=255) agaagTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCa < 1:2962195/54‑1 (MQ=255) agaagTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGt > 1:1319419/1‑57 (MQ=255) agTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCAGAGTgg > 1:2536706/1‑70 (MQ=255) agTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCAGAGTgg > 1:696885/1‑70 (MQ=255) | CGTTCCGCGACTAACGGTCTGACCCCCGCGCCTCAGGCGCGGGACTTCCAGGAAGAAGTGGCGAAACTCATCGACGTTACCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCAGAGTGG > minE/2335645‑2335770 |
| 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: ‑ |