| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,015,787 | T→C | 35.8% | R264R (CGT→CGC) | metC → | cystathionine beta‑lyase, PLP‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,015,787 | 0 | T | C | 35.8% | 47.0 / 46.3 | 78 | R264R (CGT→CGC) | metC | cystathionine beta‑lyase, PLP‑dependent |
| Reads supporting (aligned to +/- strand): ref base T (34/16); new base C (17/11); total (51/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.21e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01 | |||||||||||
AGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTT > minE/2015720‑2015853 | aGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTc < 1:1907134/69‑1 (MQ=255) tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCaa > 1:3196763/1‑68 (MQ=255) tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCaa > 1:941310/1‑68 (MQ=255) tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAAc > 1:3410918/1‑69 (MQ=255) tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAAc > 1:548854/1‑69 (MQ=255) aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa < 1:2061819/69‑1 (MQ=255) aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa < 1:2121986/69‑1 (MQ=255) aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa < 1:2456814/69‑1 (MQ=255) aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa < 1:47791/69‑1 (MQ=255) aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa < 1:3034561/69‑1 (MQ=255) aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa < 1:2849034/69‑1 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:2601667/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:2522895/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:555232/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:2113035/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:44302/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:2982367/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:33794/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:3127082/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:637813/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:1323027/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:1483998/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:1669815/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat > 1:324815/1‑60 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatca > 1:3319813/1‑59 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatca > 1:3063512/1‑59 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAgata > 1:2691615/1‑65 (MQ=255) aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTAGCGCCAAcatca > 1:1604232/1‑59 (MQ=255) ccTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGa < 1:2458232/69‑1 (MQ=255) taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg > 1:373631/1‑69 (MQ=255) taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg > 1:2870823/1‑69 (MQ=255) taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg > 1:2868453/1‑69 (MQ=255) taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg > 1:2114939/1‑69 (MQ=255) cGTGGCCTGCGCACATT‑GGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:17255/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:2752883/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:687848/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:3277069/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:452136/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:1800645/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:920105/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:592044/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:1738206/69‑1 (MQ=255) gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:3203269/69‑1 (MQ=255) tGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc < 1:2578563/68‑1 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:3092125/1‑69 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:995444/1‑69 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:2790791/1‑69 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:2407151/1‑69 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:239642/1‑69 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:219115/1‑69 (MQ=255) cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa > 1:2089824/1‑69 (MQ=255) tAGGTGTGCGTTTGGGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:1856463/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggct > 1:1761841/1‑55 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACGTcc > 1:3348394/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:802292/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:949873/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:2850928/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:1043008/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:1490955/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:1825043/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:2024544/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:3172030/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:1383641/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:2687241/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:10284/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:1265260/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:2864416/1‑67 (MQ=255) tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc > 1:2892792/1‑67 (MQ=255) aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa < 1:2229207/69‑1 (MQ=255) aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa < 1:712547/69‑1 (MQ=255) aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa > 1:1704759/1‑69 (MQ=255) tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg < 1:2629259/69‑1 (MQ=255) tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg < 1:60189/69‑1 (MQ=255) tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg < 1:3002648/69‑1 (MQ=255) tGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGa > 1:2247121/1‑68 (MQ=255) gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt < 1:500730/68‑1 (MQ=255) gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt < 1:443798/68‑1 (MQ=255) gTCAACATCATGAAAGCAGTATGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt < 1:3399390/68‑1 (MQ=255) | AGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTT > minE/2015720‑2015853 |
| 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: ‑ |