bisnake2001/local-methylation-profile
22.local-methylation-profile/SKILL.md
This skill analyzes the local DNA methylation profiles around target genomic regions provide by user. Use this skill when you want to vasulize the average methylation profile around target regions (e.g. TSS, CTCF peak or other target regions).
$ install --global
skillsauthnpx skillsauth add bisnake2001/chromskills local-methylation-profileInstall this skill globally with one command. Works with Claude Code, Cursor, and Windsurf.
Security Scan Results
3 of 9 scanners reported clean
Some scanners were skipped, did not run, or reported a non-clean status. Review each row below.
3
Scanners Passed
9
Scanners in report
Clean
TrivyContainer and dependency vulnerability scanner
95%
Clean
SemgrepStatic code analysis for vulnerabilities
95%
Clean
mcp-scan (Snyk)Model Context Protocol security validation
95%
Skipped
Snyk (dep)Open source security scanning
50%
Skipped
Socket.devSupply chain security analysis
50%
Skipped
VirusTotalMulti-engine malware detection
50%
Skipped
CrowdStrikeAdvanced threat intelligence
50%
Skipped
OSV-ScannerOpen Source Vulnerability database check
50%
Skipped
OWASP Dep-Check
50%
Last scanned: Apr 3, 2026, 7:26 PM4.3s1 file scanned
SKILL.md
- name:
- local-methylation-profile
- description:
- This skill analyzes the local DNA methylation profiles around target genomic regions provide by user. Use this skill when you want to vasulize the average methylation profile around target regions (e.g. TSS, CTCF peak or other target regions).
Local Methylation Profile Analysis
Overview
- Always prompt user for which columns in the BED files are methylation fraction/percent. Never decide by yourself.
- Generat profile: Bin methylation around regions (±flank, fixed bin size), aggregate mean±SE.
- Visualize: Plot mean profile with ribbon and center line.
Inputs & Outputs
Inputs
methylation.bed
target_regions.bed
Outputs
local_methyl_profile/
stats/
CpG_around_target.tsv
plots/
CpG_around_target.pdf
temp/
... # other temp file generated
Decision Tree
Step 1: Preprocess input → 5-column BED (for methylKit), and 3-column BED (for target regions)
awk -F'\t' 'BEGIN {OFS="\t"} {print $1, $2, $3, $<i_methylation>}, $<i_coverage>}' methylation.bed # n is provide by user, *100 if is fraction
awk -F'\t' 'BEGIN {OFS="\t"} {print $1, $2, $3}' target_regions.bed
Step 2: Build methylation profiles around regions
Call:
mcp__methyl-tools__build_local_methylation_profile
with:
methyl_bed_path: 5-column BED-like file from preprocess_methylation.regions_bed_path: 3-column BED-like file from preprocess_regions.output_profile_tsv_path: path for aggregated profile table (TSV).flank_size: flank size in bp around region center (default 2000).bin_size: bin size in bp (default 50).min_coverage: minimum coverage threshold for CpGs (default 10).
Step 3: Visualization
Call:
mcp__methyl-tools__plot_profile
with:
profile_tsv_path: TSV from build_methylation_profile.output_plot_path: output figure path (PNG/PDF; format inferred from extension).title: plot title (optional).
Parameter Guidelines
| Context | Flank | Bin | Min cov | |-----------|-------|------|---------| | TF peaks | ±2 kb | 50bp | 10x | | Promoters | ±1 kb | 50bp | 10x | | Enhancers | ±5 kb | 100bp| 5x | | Motifs | ±0.5kb| 10–20| 10x |
Notes
- Snippets are usage hints and must be adapted to your paths and column indices.
Related Skills
bisnake2001/reads-mapping
development
VerifiedTrustedCommunity
Align ChIP-seq or ATAC-seq FASTQ files to a reference genome using Bowtie2, with strict input validation, library layout detection, output organization and logging. Use it when raw sequencing reads must be converted into sorted/indexed BAM files before downstream QC, peak calling, or footprinting.
5SKILL.mdUpdated Apr 29, 2026
bisnake2001/dna-methylation-alignment-bismark
development
VerifiedTrustedCommunity
Align bisulfite sequencing DNA methylation reads using Bismark only, with explicit validation of reference preparation, library layout detection, output organization, logging, and alignment QC. Use it for WGBS, RRBS, or other bisulfite-converted DNA methylation sequencing data when raw FASTQ files must be aligned before methylation extraction and downstream analysis.
5SKILL.mdUpdated Apr 26, 2026
bisnake2001/peak-calling
data-ai
VerifiedTrustedCommunity
Perform peak calling for ChIP-seq or ATAC-seq data using MACS3, with intelligent parameter detection from user feedback. Use it when you want to call peaks for ChIP-seq data or ATAC-seq data.
4SKILL.mdUpdated Apr 23, 2026
bisnake2001/TF-differential-binding
devops
VerifiedTrustedCommunity
The TF-differential-binding pipeline performs differential transcription factor (TF) binding analysis from ChIP-seq datasets (TF peaks) using the DiffBind package in R. It identifies genomic regions where TF binding intensity significantly differs between experimental conditions (e.g., treatment vs. control, mutant vs. wild-type). Use the TF-differential-binding pipeline when you need to analyze the different function of the same TF across two or more biological conditions, cell types, or treatments using ChIP-seq data or TF binding peaks. This pipeline is ideal for studying regulatory mechanisms that underlie transcriptional differences or epigenetic responses to perturbations.
4SKILL.mdUpdated Apr 3, 2026