| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,299,421 | G→A | 13.8% | intergenic (‑346/+144) | cobU ← / ← yeeX | bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,299,421 | 0 | G | A | 13.8% | 152.0 / 12.7 | 58 | intergenic (‑346/+144) | cobU/yeeX | bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (26/24); new base A (2/6); total (28/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.56e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGA > minE/1299354‑1299488 | gCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaa < 1:1654604/71‑1 (MQ=255) ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACg > 1:753072/1‑70 (MQ=255) ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGt < 1:3902214/71‑1 (MQ=255) cTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTaa > 1:3636786/1‑71 (MQ=255) cTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTaa > 1:1557396/1‑71 (MQ=255) gAAATCCAGGCTATTGATTAAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTaa > 1:2106251/1‑69 (MQ=255) aaTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAAt > 1:3110752/1‑68 (MQ=255) aTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGaa > 1:2986970/1‑71 (MQ=255) tCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAAt < 1:1399187/71‑1 (MQ=255) ccAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAaac > 1:3733922/1‑57 (MQ=255) cAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa < 1:1853309/62‑1 (MQ=255) aGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaa < 1:14849/52‑1 (MQ=255) gCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATcc < 1:2646719/71‑1 (MQ=255) cTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGa < 1:3770518/63‑1 (MQ=255) aTTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCgg > 1:2517659/1‑70 (MQ=255) ttGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGt < 1:450261/44‑1 (MQ=255) ttGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCg < 1:3277954/68‑1 (MQ=255) tGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAt > 1:1415213/1‑71 (MQ=255) ggttCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGAtg > 1:513489/3‑71 (MQ=255) gATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtg > 1:387079/1‑71 (MQ=255) gATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAt < 1:1958453/70‑1 (MQ=255) aTTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAg < 1:3181737/56‑1 (MQ=255) ttCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaaa < 1:156618/43‑1 (MQ=255) ttCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTAt < 1:2362783/61‑1 (MQ=255) tCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg > 1:1656376/1‑71 (MQ=255) cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa > 1:2355526/1‑49 (MQ=255) cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa > 1:2454941/1‑49 (MQ=255) cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTAt > 1:3089553/1‑59 (MQ=255) cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg < 1:299404/70‑1 (MQ=255) agcGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGc < 1:1515239/69‑1 (MQ=255) aaCGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggt > 1:2573165/1‑68 (MQ=255) aaCGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGc > 1:1999654/1‑71 (MQ=255) cGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGc > 1:2213787/1‑69 (MQ=255) tCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAAt < 1:546458/46‑1 (MQ=255) tCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACt < 1:176595/71‑1 (MQ=255) tCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGCGagc < 1:2394593/71‑3 (MQ=255) gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatg < 1:3822283/71‑2 (MQ=255) gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatg < 1:3785583/71‑2 (MQ=255) cGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaa > 1:1357433/1‑71 (MQ=255) gACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatgca > 1:1481928/1‑68 (MQ=255) gATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGa < 1:3212400/50‑1 (MQ=255) gtCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAg < 1:490054/37‑1 (MQ=12) gtCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAg < 1:1636649/37‑1 (MQ=12) aTCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAg < 1:3282274/38‑1 (MQ=255) aTCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGa > 1:1018209/1‑71 (MQ=255) tCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATcc > 1:33409/1‑45 (MQ=255) tcgcTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCgg < 1:1036078/43‑1 (MQ=255) tcgcTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCgg < 1:1508678/43‑1 (MQ=255) cgcTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGTAATGGTGGCGACTATAATTTCGATTg < 1:1071593/69‑1 (MQ=255) cgcTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGAt < 1:3443648/67‑1 (MQ=255) cTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTtc > 1:321380/1‑71 (MQ=255) aaaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCg > 1:433460/1‑70 (MQ=255) aaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtgg > 1:2604181/1‑44 (MQ=255) gAACTTGTAAAACAACGTAATAGAATTATCCGGATTGGTGGCGACTATAATTTCGATTGGTTCTCggat > 1:3381622/1‑69 (MQ=255) aaCTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTgg < 1:3322195/59‑1 (MQ=255) aaCTTGTAAAACAACGTAATAGAATTATCCGGAATGGAGGCGACTATAATTTCGATTGGTTCTCggatgag > 1:1299343/1‑71 (MQ=255) cTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg > 1:3760524/1‑37 (MQ=255) cTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgagga < 1:3250968/71‑1 (MQ=255) | GCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGA > minE/1299354‑1299488 |
| 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: ‑ |