3838 const path_to_root = "" ;
3939 const default_light_theme = "light" ;
4040 const default_dark_theme = "navy" ;
41- window . path_to_searchindex_js = "searchindex-b1d0ec24 .js" ;
41+ window . path_to_searchindex_js = "searchindex-ec6434a5 .js" ;
4242 </ script >
4343 <!-- Start loading toc.js asap -->
44- < script src ="toc-a172ac56 .js "> </ script >
44+ < script src ="toc-4ae5f900 .js "> </ script >
4545 </ head >
4646 < body >
4747 < div id ="mdbook-help-container ">
@@ -2219,14 +2219,17 @@ <h2 id="extract-calls"><a class="header" href="#extract-calls">extract calls</a>
22192219< h2 id ="extract-read-stats "> < a class ="header " href ="#extract-read-stats "> extract read-stats</ a > </ h2 >
22202220< pre > < code class ="language-text "> Produce a table where modification counts are summarized on the read level. This
22212221table will have one record per valid read and count the number of modified and
2222- unmodified bases for each base moficiation requested
2222+ unmodified bases for each base modification requested
22232223
22242224Usage: modkit extract read-stats [OPTIONS] <IN_BAM> <OUT_PATH>
22252225
22262226Arguments:
2227- <IN_BAM> Input modBAM. Can be a path to a file or "-"/"stdin" for standard
2228- input
2229- <OUT_PATH> Path to file for output, default will be to stdout
2227+ <IN_BAM>
2228+ Input modBAM. Can be a path to a file or "-"/"stdin" for standard
2229+ input
2230+
2231+ <OUT_PATH>
2232+ Path to file for output, default will be to stdout
22302233
22312234Options:
22322235 --filter-threshold <FILTER_THRESHOLD>
@@ -2240,26 +2243,37 @@ <h2 id="extract-read-stats"><a class="header" href="#extract-read-stats">extract
22402243 default for all other bases with: --filter-threshold A:0.70
22412244 --filter-threshold 0.9 will specify a threshold value of 0.70 for
22422245 adenine and 0.9 for all other base modification calls
2246+
22432247 --mod-thresholds <MOD_THRESHOLDS>
22442248 Specify a passing threshold to use for a base modification,
22452249 independent of the threshold for the primary sequence base or the
22462250 default. For example, to set the pass threshold for 5hmC to 0.8 use
22472251 `--mod-threshold h:0.8`. The pass threshold will still be estimated as
22482252 usual and used for canonical cytosine and other modifications unless
22492253 the `--filter-threshold` option is also passed
2254+
22502255 --mod-codes <MOD_CODES>...
22512256 Specify which modcodes to tabulate, should be
22522257 <primary_base>:<mod_code>
2258+
22532259 --queue-size <QUEUE_SIZE>
22542260 Queue size, maximum number of reads to be held in memory waiting for
2255- records to be written [default: 10000]
2261+ records to be written
2262+
2263+ [default: 10000]
2264+
22562265 --io-threads <IO_THREADS>
22572266 [default: 4]
2267+
22582268 -h, --help
2259- Print help
2269+ Print help (see a summary with '-h')
22602270
22612271Logging Options:
2262- --log-filepath <LOG_FILEPATH> Path to file to write run log
2272+ --log-filepath <LOG_FILEPATH>
2273+ Path to file to write run log
2274+
2275+ --suppress-progress
2276+ Hide the progress bar
22632277</ code > </ pre >
22642278< h2 id ="motif-bed "> < a class ="header " href ="#motif-bed "> motif bed</ a > </ h2 >
22652279< pre > < code class ="language-text "> Create BED file with all locations of a sequence motif. Example: modkit motif
@@ -2900,6 +2914,186 @@ <h2 id="dmr-multi"><a class="header" href="#dmr-multi">dmr multi</a></h2>
29002914 --io-threads <IO_THREADS>
29012915 Number of threads to use when for decompression [default: 4]
29022916</ code > </ pre >
2917+ < h2 id ="dmr-isoform "> < a class ="header " href ="#dmr-isoform "> dmr isoform</ a > </ h2 >
2918+ < pre > < code class ="language-text "> Use a transcriptome-aligned bedMethyl table and GTF file to compare methylation
2919+ across gene transcript isoforms. Run analysis on the entire transcriptome or on
2920+ a single gene. Optionally also make plots for single genes
2921+
2922+ Usage: modkit dmr isoform [OPTIONS] --gtf <GTF> <INPUT_BEDMETHYL> <OUTPUT_TABLE>
2923+
2924+ Arguments:
2925+ <INPUT_BEDMETHYL>
2926+ Input bedmethyl table. Expected to be bgzip-compressed and have a .tbi
2927+ tabix index. This bedMethyl table must be generated from a
2928+ transcriptome-aligned modBAM, NOT a genome-aligned modBAM. Expected
2929+ that this table was generated by pileup with `--modified-bases`
2930+ argument
2931+
2932+ <OUTPUT_TABLE>
2933+ Path to output BED file or '-' for standard output
2934+
2935+ Options:
2936+ --gtf <GTF>
2937+ Path to GTF file to use for transcript models. May be compressed
2938+
2939+ --min-valid-coverage <MIN_VALID_COVERAGE>
2940+ Minimum valid coverage required to use an entry from a bedMethyl. See
2941+ the help for pileup for the specification and description of valid
2942+ coverage
2943+
2944+ [default: 1]
2945+
2946+ --single-code <SINGLE_MODIFICATION_CODE>
2947+ Calculate differential methylation for this modification code in
2948+ isolation. Unmodified and other modification counts will be combined
2949+
2950+ --full
2951+ Include per-modification combined proportions per position and
2952+ per-isoform proportions in output. Adding this flag will make the
2953+ output file substantially larger
2954+
2955+ --log-filepath <LOG_FILEPATH>
2956+ Path to optional debug log (recommended)
2957+
2958+ -g, --gene-id <GENE_ID>
2959+ Only process this Gene Id
2960+
2961+ -n, --gene-name <GENE_NAME>
2962+ Only process this Gene Name
2963+
2964+ --suppress-progress
2965+ Don't show the progress bar
2966+
2967+ -t, --threads <THREADS>
2968+ Number of concurrent processors to use. Only used when running on all
2969+ genes
2970+
2971+ [default: 4]
2972+
2973+ --plot <PLOT>
2974+ Generate an SVG plot displaying the methylation positions on each
2975+ isoform (filtered by --min-score or --max-pval) and the common exons
2976+ of the gene. Also displays negative log of the p-value of the DMR
2977+ metric. See the online documentation for an example
2978+
2979+ --marker-positions <MARKER_POSITIONS>...
2980+ Comma-separated list of genomic positions to draw a vertical dotted
2981+ line and a label
2982+
2983+ --width <WIDTH>
2984+ SVG width in px
2985+
2986+ [default: 1600]
2987+
2988+ --track-height <TRACK_HEIGHT>
2989+ Height per transcript track in px
2990+
2991+ [default: 26]
2992+
2993+ --track-gap <TRACK_GAP>
2994+ Gap between transcript tracks in px
2995+
2996+ [default: 18]
2997+
2998+ --left-margin <LEFT_MARGIN>
2999+ Left margin in px
3000+
3001+ [default: 220]
3002+
3003+ --right-margin <RIGHT_MARGIN>
3004+ Right margin in px
3005+
3006+ [default: 40]
3007+
3008+ --top-margin <TOP_MARGIN>
3009+ Top margin in px
3010+
3011+ [default: 210]
3012+
3013+ --bottom-margin <BOTTOM_MARGIN>
3014+ Bottom margin in px
3015+
3016+ [default: 50]
3017+
3018+ --intron-px <INTRON_PX>
3019+ Length, in px, of the compressed intron segments
3020+
3021+ [default: 24]
3022+
3023+ --top-bar-height <TOP_BAR_HEIGHT>
3024+ Height in px of the top bar
3025+
3026+ [default: 90]
3027+
3028+ --min-score <MIN_SCORE>
3029+ Only plot methylation markers with negative log p-value greater than
3030+ this value
3031+
3032+ --max-pval <MAX_PVAL>
3033+ Only plot methylation markers with p-value less than this value
3034+
3035+ -h, --help
3036+ Print help (see a summary with '-h')
3037+ </ code > </ pre >
3038+ < h2 id ="dmr-compare-tx-sites "> < a class ="header " href ="#dmr-compare-tx-sites "> dmr compare-tx-sites</ a > </ h2 >
3039+ < pre > < code class ="language-text "> Uses two transcriptome-aligned bedmMethyl tables to compare methylation for two
3040+ conditions at transcriptome sites aggregated at the gene level. This command
3041+ will combine the methylation counts across isoforms into a single summary
3042+ statistic per-site for the gene, then compare those counts between conditions at
3043+ the gene level and report results on genomic coordinates. The transcript models
3044+ are provided by a GTF file. Optionally make a volcano plot of the top-k
3045+ differentially methylated sites per gene
3046+
3047+ Usage: modkit dmr compare-tx-sites [OPTIONS] -a <COND_A> -b <COND_B> --out <OUTPUT_TABLE> --gtf <GTF>
3048+
3049+ Options:
3050+ -a <COND_A>
3051+ Input bedmethyl table for first condition (condition 'A'). Expected to
3052+ be bgzip-compressed and have a .tbi tabix index. This bedMethyl table
3053+ must be generated from a transcriptome-aligned modBAM, NOT a
3054+ genome-aligned modBAM. Expected that this table was generated by
3055+ pileup with `--modified-bases` argument
3056+ -b <COND_B>
3057+ Input bedmethyl table for first condition (condition 'B')
3058+ -o, --out <OUTPUT_TABLE>
3059+ Path to output BED file or '-' for stdout
3060+ --full
3061+ Include combined proportions and counts for each condition. This will
3062+ increase the size of the output file
3063+ --gtf <GTF>
3064+ Path to GTF file to use for transcript models. Can be compressed
3065+ --min-valid-coverage <MIN_VALID_COVERAGE>
3066+ Minimum valid coverage required to use an entry from a bedMethyl. See
3067+ the help for pileup for the specification and description of valid
3068+ coverage [default: 1]
3069+ --single-code <SINGLE_MODIFICATION_CODE>
3070+ Calculate differential methylation for this modification code in
3071+ isolation. Unmodified and other modification counts will be combined
3072+ --log-filepath <LOG_FILEPATH>
3073+ Path to optional debug log (recommended)
3074+ --suppress-progress
3075+ Don't show the progress bar
3076+ -t, --threads <THREADS>
3077+ Number of concurrent processors to use [default: 4]
3078+ --plot <PLOT>
3079+ Path to file to create an SVG volcano plot where the x-axis is log2
3080+ fold change and the y-axis is negative log p-value. See the online
3081+ documentation for and example
3082+ --top-k <TOP_K>
3083+ For each gene, sort the positions in decreasing negative log p-value
3084+ and take this many for plotting. Setting a large number will increase
3085+ the size of the SVG plot [default: 1]
3086+ --min-effect-size <MIN_EFFECT_SIZE>
3087+ Discard positions with an effect size less than this before
3088+ considering them for plotting [default: 0.1]
3089+ --plot-title <PLOT_TITLE>
3090+ Add a title to the plot
3091+ --gene-labels <GENE_LABELS>
3092+ Path to plain text file of a list of gene names, one per line, to mark
3093+ in the volcano plot
3094+ -h, --help
3095+ Print help
3096+ </ code > </ pre >
29033097< h2 id ="bedmethyl-merge "> < a class ="header " href ="#bedmethyl-merge "> bedmethyl merge</ a > </ h2 >
29043098< pre > < code class ="language-text "> Perform an outer join on two or more bedMethyl files, summing their counts for
29053099records that overlap
@@ -3051,6 +3245,22 @@ <h2 id="bedmethyl-tobigwig"><a class="header" href="#bedmethyl-tobigwig">bedmeth
30513245 --suppress-progress
30523246 Hide the progress bar
30533247</ code > </ pre >
3248+ < h2 id ="bedmethyl-map-to-genome "> < a class ="header " href ="#bedmethyl-map-to-genome "> bedmethyl map-to-genome</ a > </ h2 >
3249+ < pre > < code class ="language-text "> Map a transcriptome-aligned bedMethyl file to genome-coordinates based on a GTF
3250+
3251+ Usage: modkit bedmethyl map-to-genome [OPTIONS] --transcript-id <TRANSCRIPT_ID> --gtf <GTF> <INPUT_BEDMETHYL> <OUTPUT_BEDMETHYL>
3252+
3253+ Arguments:
3254+ <INPUT_BEDMETHYL>
3255+ <OUTPUT_BEDMETHYL>
3256+
3257+ Options:
3258+ -t, --transcript-id <TRANSCRIPT_ID>
3259+ --transcript-version <TRANSCRIPT_VERSION>
3260+ --gtf <GTF>
3261+ --header
3262+ -h, --help Print help
3263+ </ code > </ pre >
30543264< h2 id ="modbam-adjust-mods "> < a class ="header " href ="#modbam-adjust-mods "> modbam adjust-mods</ a > </ h2 >
30553265< pre > < code class ="language-text "> Performs various operations on BAM files containing base modification
30563266information, such as converting base modification codes and ignoring
0 commit comments