Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 899,672 | T→C | 8.3% | L15P (CTG→CCG) | pspE → | thiosulfate:cyanide sulfurtransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 899,672 | 0 | T | C | 8.3% | 194.9 / 9.3 | 72 | L15P (CTG→CCG) | pspE | thiosulfate:cyanide sulfurtransferase |
Reads supporting (aligned to +/- strand): ref base T (27/39); new base C (4/2); total (31/41) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.92e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01 |
TTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTT > minE/899606‑899739 | ttcattcaTAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc > 1:1984197/1‑71 (MQ=255) ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc > 1:1418498/1‑71 (MQ=255) ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc > 1:3079218/1‑71 (MQ=255) ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc > 1:3275012/1‑71 (MQ=255) cATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGt > 1:1828104/1‑71 (MQ=255) cATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc < 1:3005134/71‑1 (MQ=255) cATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc < 1:520483/71‑1 (MQ=255) tAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc < 1:3014207/69‑1 (MQ=255) gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgccg < 1:49542/71‑1 (MQ=255) gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgccg < 1:398955/71‑1 (MQ=255) aaGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgccg < 1:3491836/70‑1 (MQ=255) aCGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGtt < 1:661903/59‑1 (MQ=255) cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGaa > 1:3876223/1‑70 (MQ=255) cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGaa > 1:3659352/1‑70 (MQ=255) gCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGtttt > 1:3254336/1‑59 (MQ=255) tACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTg > 1:2800224/1‑71 (MQ=255) tACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTg > 1:3760399/1‑71 (MQ=255) aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTg > 1:3264434/1‑70 (MQ=255) aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTg > 1:61421/1‑70 (MQ=255) tttAAAAAATGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:3738658/71‑1 (MQ=255) tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTgg < 1:2966973/66‑1 (MQ=255) tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa > 1:2384868/1‑71 (MQ=255) tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:2547826/71‑1 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGAt > 1:2123320/1‑66 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGAt > 1:282300/1‑66 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:1930860/1‑71 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:209792/1‑71 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:3051776/1‑71 (MQ=255) aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGt < 1:1036882/70‑1 (MQ=255) aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGt < 1:1003257/70‑1 (MQ=255) aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGt < 1:168166/70‑1 (MQ=255) aaaaGGCTTACTTGCTCTGTCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTTGATCGATGTTc < 1:1717593/71‑1 (MQ=255) ttACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTc > 1:1688884/1‑70 (MQ=255) tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:3362055/1‑70 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTc > 1:3097274/1‑64 (MQ=255) gCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:3352755/1‑69 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:3894283/47‑1 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:1799197/47‑1 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:1803889/47‑1 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:1935866/47‑1 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:1280034/47‑1 (MQ=255) tGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGt < 1:3243145/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:486017/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:1312394/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:2936017/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:2918211/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt > 1:2877704/1‑69 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:1780944/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt > 1:3680247/1‑69 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCa < 1:1379736/71‑1 (MQ=255) cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:918905/44‑1 (MQ=255) cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:1539828/44‑1 (MQ=255) cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:3299230/44‑1 (MQ=255) cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:2597291/44‑1 (MQ=255) cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa < 1:2388946/44‑1 (MQ=255) gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCaa > 1:496611/1‑69 (MQ=255) tgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAg < 1:3828434/59‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCATAGCAGTATCAGc < 1:2605134/64‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAg < 1:3365759/57‑1 (MQ=255) tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGc < 1:2847840/63‑1 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCTATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg > 1:757099/1‑71 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg > 1:3905905/1‑44 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg > 1:3912010/1‑44 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg > 1:2279502/1‑44 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg > 1:440072/1‑71 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg < 1:3221734/70‑1 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg < 1:798017/70‑1 (MQ=255) cACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt > 1:420133/1‑71 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGTATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt < 1:1848761/70‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt < 1:3885477/70‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt < 1:2685603/70‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGAGTTCCAGAGCAGTATCAGCAAGAGCACGtt < 1:2317800/70‑1 (MQ=255) | TTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTT > minE/899606‑899739 |
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: ‑ |