Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,299,417 | A→G | 13.8% | intergenic (‑342/+148) | 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,417 | 0 | A | G | 13.8% | 145.1 / 14.3 | 58 | intergenic (‑342/+148) | cobU/yeeX | bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (25/25); new base G (2/6); total (27/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01 |
AAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAG > minE/1299348‑1299486 | aaaCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAAc < 1:2359040/71‑1 (MQ=255) aaCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACt < 1:358104/71‑1 (MQ=255) aCGGGCATTGTGTTCTGAAATCGAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACtt < 1:3273185/71‑1 (MQ=255) gCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaa < 1:1654604/71‑1 (MQ=255) ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACg > 1:753072/1‑70 (MQ=255) ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGt < 1:3902214/71‑1 (MQ=255) tCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAAc > 1:2582390/1‑56 (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:2454941/1‑49 (MQ=255) cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa > 1:2355526/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:3785583/71‑2 (MQ=255) gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatg < 1:3822283/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) | AAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAG > minE/1299348‑1299486 |
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: ‑ |