Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 1,761,776 | C→T | 71.9% | L154L (CTG→CTA) | sufS ← | selenocysteine lyase, PLP‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 1,761,776 | 0 | C | T | 71.9% | 45.2 / 21.7 | 32 | L154L (CTG→CTA) | sufS | selenocysteine lyase, PLP‑dependent |
Reads supporting (aligned to +/- strand): ref base C (4/5); new base T (12/11); total (16/16) | |||||||||||
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.07e-01 |
CGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCAT > W3110S.gb/1761712‑1761844 | cGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGt > 1:373250/1‑68 (MQ=255) cGTTGGAGACATGAGTAATAGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGt > 1:2365458/1‑68 (MQ=255) gaCATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACg < 1:282646/68‑1 (MQ=255) gaCATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACg < 1:2693419/68‑1 (MQ=255) gaCATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGt > 1:110948/1‑69 (MQ=255) aCATGAGTTATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTa < 1:25451/69‑1 (MQ=255) aCATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTa < 1:2882570/69‑1 (MQ=255) cATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTAc < 1:3563199/69‑1 (MQ=255) cATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTAc < 1:3235738/69‑1 (MQ=255) aTTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGa > 1:326050/1‑69 (MQ=255) gCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCggg > 1:442089/1‑65 (MQ=255) cagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCgg > 1:827450/1‑69 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCggg > 1:830757/1‑59 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCggg > 1:419928/1‑59 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCggg > 1:2716152/1‑59 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCggg > 1:45985/1‑59 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCggg < 1:1353794/69‑1 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCggg < 1:481982/69‑1 (MQ=255) agGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGATGCAACGTACCATCGGGATTGAGCggg < 1:2253321/69‑1 (MQ=255) gCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCGGGAt < 1:1497171/69‑1 (MQ=255) gAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCgg < 1:1637493/64‑1 (MQ=255) gAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCGGGATca > 1:1894869/1‑69 (MQ=255) gAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCGGGATca > 1:3022181/1‑69 (MQ=255) caAACAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGctc < 1:1078918/67‑1 (MQ=255) aaaCAGCGTAGGCAGCGTCTCTAGTTGCAACGTAc > 1:1728436/1‑35 (MQ=255) cAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCgg > 1:2965962/1‑38 (MQ=255) cAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCgg > 1:2525599/1‑38 (MQ=255) cAGCGTAGGCAGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGcc < 1:3257063/69‑1 (MQ=255) gTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAAcg < 1:2040109/69‑1 (MQ=255) aGCGTCTCTAGTTGCAACGTACCATCGGGATTGAGCGGGATca > 1:848499/1‑43 (MQ=255) aGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTgcg < 1:1424774/69‑1 (MQ=255) cGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCa < 1:346074/69‑1 (MQ=255) cAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCAt > 1:3559664/1‑69 (MQ=255) cAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCAt > 1:2041273/1‑69 (MQ=255) | CGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCAT > W3110S.gb/1761712‑1761844 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |