| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,015,838 | G→C | 24.8% | P281P (CCG→CCC) | 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,838 | 0 | G | C | 24.8% | 200.8 / 67.6 | 125 | P281P (CCG→CCC) | metC | cystathionine beta‑lyase, PLP‑dependent |
| Reads supporting (aligned to +/- strand): ref base G (14/80); new base C (19/12); total (33/92) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.57e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 | |||||||||||
Only 100 of 125 total aligned reads displayed.
AGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAA > minE/2015772‑2015897 | aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa < 1:2229207/69‑1 (MQ=255)aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa > 1:1704759/1‑69 (MQ=255) tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg < 1:60189/69‑1 (MQ=255) gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt < 1:500730/68‑1 (MQ=255) gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt < 1:443798/68‑1 (MQ=255) gTCAACATCATGAAAGCAGTATGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt < 1:3399390/68‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:1086066/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:2915920/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:1208876/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:3042189/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:3059581/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:2509861/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:136949/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:1630664/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac < 1:2361095/69‑1 (MQ=255) cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCCAGTTAac < 1:1038301/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:1519952/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:708550/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:2424879/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:3399263/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:3321296/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:1236079/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac < 1:990212/69‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:3429090/68‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:359847/68‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:2551652/68‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:2176321/68‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:2923755/68‑1 (MQ=255) catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:488700/68‑1 (MQ=255) atGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTg > 1:408373/1‑69 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTTGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:1294104/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:2971152/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:628064/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:2121903/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:1635444/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:3183957/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:1414412/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc < 1:2479816/68‑1 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGtt > 1:2487544/1‑57 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTg > 1:2303020/1‑67 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:1094467/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:283834/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:3011618/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:2083284/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:1305925/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:1340347/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:134188/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:2484602/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:758152/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:2315036/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:1655359/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:1692589/1‑68 (MQ=255) gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc > 1:2203668/1‑68 (MQ=255) cTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca < 1:1854868/53‑1 (MQ=255) tGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg > 1:836475/1‑43 (MQ=255) aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa < 1:1491898/69‑1 (MQ=255) aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa < 1:3318333/69‑1 (MQ=255) aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa < 1:2706342/69‑1 (MQ=255) aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa < 1:546135/69‑1 (MQ=255) gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa < 1:1060870/68‑1 (MQ=255) gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa < 1:2437436/68‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:2931938/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:1162778/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:982790/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:745524/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:737135/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:637286/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:2266201/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:322791/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:3157357/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:2928959/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:3069843/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:2970741/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:2911939/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:2751695/49‑1 (MQ=255) ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt < 1:273948/49‑1 (MQ=255) gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt > 1:3300975/1‑69 (MQ=255) gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt > 1:3131526/1‑69 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:446231/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:746635/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:1587551/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:561828/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:2038803/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:2216875/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:230316/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:2641345/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGAGTTAAACACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:2519435/67‑1 (MQ=255) gctggcAGAACATCCCCAAGTAGCGCGACTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa < 1:2296054/67‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:1247180/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:1889297/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:963995/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:1867486/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:1754557/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:3040452/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa < 1:157391/69‑1 (MQ=255) aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGa < 1:2009794/68‑1 (MQ=255) aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa > 1:444014/1‑66 (MQ=255) aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa > 1:2231118/1‑66 (MQ=255) aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa > 1:2502209/1‑66 (MQ=255) aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa > 1:2524745/1‑66 (MQ=255) | AGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAA > minE/2015772‑2015897 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |