Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,610,889 | G→A | 46.5% | C283C (TGC→TGT) | eutB ← | ethanolamine ammonia‑lyase, large subunit, heavy chain |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,610,889 | 0 | G | A | 46.5% | 2.9 / 33.1 | 43 | C283C (TGC→TGT) | eutB | ethanolamine ammonia‑lyase, large subunit, heavy chain |
Reads supporting (aligned to +/- strand): ref base G (14/9); new base A (7/13); total (21/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01 |
AGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCT > minE/1610844‑1610937 | aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt > 1:2779911/1‑68 (MQ=255) aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt > 1:1353132/1‑68 (MQ=255) aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt > 1:1355978/1‑68 (MQ=255) aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt > 1:3273839/1‑68 (MQ=255) aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt > 1:2245818/1‑68 (MQ=255) aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt > 1:2203472/1‑68 (MQ=255) tAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGcc < 1:2002504/69‑1 (MQ=255) tAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGcc < 1:1998471/69‑1 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:1269724/1‑69 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:3192928/1‑69 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:2768371/1‑69 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:2537355/1‑69 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:82689/1‑69 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:1569268/1‑69 (MQ=255) cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGGTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc > 1:3220171/1‑69 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:2951729/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:1344329/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:1039022/44‑1 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:2732979/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:1481749/44‑1 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:2696627/44‑1 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:255887/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:1877697/44‑1 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:3301059/44‑1 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:2472639/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:2356273/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:827701/44‑1 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:2242668/1‑44 (MQ=255) gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc > 1:3404716/1‑44 (MQ=255) ccTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc < 1:2834453/43‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:828916/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:996661/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:568803/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:2860712/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:2803665/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:2711342/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:2477371/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:2059545/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:1695411/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:159905/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:1539465/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:1504128/69‑1 (MQ=255) tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt < 1:1469656/69‑1 (MQ=255) | AGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCT > minE/1610844‑1610937 |
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: ‑ |