Skip to content
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
15 changes: 12 additions & 3 deletions book/src/intro_dmr_isoform.md
Original file line number Diff line number Diff line change
Expand Up @@ -25,11 +25,15 @@ To generate an isoform SVG plot like the one above, add the following arguments:
modkit dmr isoform \
${bedmethyl} \
isoform_dmr_${gene_name}.bed \
--plot path/to/plot/dir \ # <-- add this option
--gene-name ${gene_name} \ # <-- requires a gene name (or gene id to plot)
--plot path/to/plot/dir \
--gene-name ${gene_name} \
--gtf ${gtf}
```

> **Important:** `--plot` requires either `--gene-name` or `--gene-id` and **cannot** be used
> in whole-transcriptome mode (i.e. without one of those arguments). Attempting to use `--plot`
> without a gene specifier will result in an error.

*Note* that exons without marks in certain transcripts indicates that there is no data in the input bedMethyl.
The command `modkit bedmethyl map-to-genome` can map a transcript-aligned bedMethyl to genome coordinates.

Expand Down Expand Up @@ -63,6 +67,11 @@ Through alternative splicing gene exons can be combined in multiple permutations

## Filtering the number of methylation marks
Some long genes may have many modified positions and drawing them all can get crowded.
To only plot positions with a maximum p-value or minimum score use the `--max-pvalue` or `--min-score` arguments, respectively.
To only plot positions with a maximum p-value or minimum score use the `--max-pval` or
`--min-score` arguments, respectively.

> **Note:** `--max-pval` and `--min-score` are **plot-only** options. They require `--plot`
> (and therefore also `--gene-name` or `--gene-id`) and have no effect outside of single-gene
> mode. These two flags are **mutually exclusive** — only one may be provided at a time.
Comment thread
SuhasSrinivasan marked this conversation as resolved.