| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,299,448 | A→G | 10.5% | intergenic (‑373/+117) | 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,448 | 0 | A | G | 10.5% | 98.5 / 4.0 | 38 | intergenic (‑373/+117) | cobU/yeeX | bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (19/15); new base G (2/2); total (21/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01 | |||||||||||
TGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCA > minE/1299380‑1299517 | 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) tCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg > 1:1656376/1‑71 (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) 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) aTCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGa > 1:1018209/1‑71 (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) tAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGgatgat < 1:1762734/71‑1 (MQ=255) gTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACgg < 1:2348389/67‑1 (MQ=255) aTAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca > 1:3738676/1‑71 (MQ=255) aTAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca > 1:92975/1‑71 (MQ=255) tAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca > 1:3796255/1‑70 (MQ=255) aGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca > 1:2246264/1‑69 (MQ=255) aaTTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggatg < 1:1116850/53‑1 (MQ=255) aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggatg < 1:2820589/49‑1 (MQ=255) aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggatg < 1:3350800/49‑1 (MQ=255) aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGAt > 1:3014685/1‑68 (MQ=255) gAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATca < 1:2027807/71‑1 (MQ=255) | TGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCA > minE/1299380‑1299517 |
| 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: ‑ |