bio-epitranscriptomics-m6a-differential
Identify differential m6A methylation between conditions from MeRIP-seq. Use when comparing epitranscriptomic changes between treatment groups or cell states.
npx skills add majiayu000/claude-skill-registry --skill m6a-differential --agent claude-code
Same command for any agent — swap --agent for codex, cursor, copilot.
Weekly change comes from our own snapshots, not the repository page — it measures attention, not adoption.
# Differential m6A Analysis ## exomePeak2 Differential Analysis ```r library(exomePeak2) # Define sample design # condition: factor for comparison design <- data.frame( condition = factor(c('ctrl', 'ctrl', 'treat', 'treat')) ) # Differential peak calling result <- exomePeak2( bam_ip = c('ctrl_IP1.bam', 'ctrl_IP2.bam', 'treat_IP1.bam', 'treat_IP2.bam'), bam_input = c('ctrl_Input1.bam', 'ctrl_Input2.bam', 'treat_Input1.bam', 'treat_Input2.bam'), gff = 'genes.gtf', genome = 'hg38', experiment_design = design ) # Get differential sites diff_sites <- results(result, contrast = c('condition', 'treat', 'ctrl')) ``` ## QNB for Differential Methylation ```r library(QNB) # Requires count matrices from peak regions # IP and input counts per sample qnb_result <- qnbtest( IP_count_matrix, Input_count_matrix, group = c(1, 1, 2, 2) # 1=ctrl, 2=treat ) # Filter significant # padj < 0.05, |log2FC| > 1 sig <- qnb_result[qnb_result$padj < 0.05 & abs(qnb_result$log2FC) > 1, ] ``` ## Visualization ```r library(ggplot2) # Volcano plot ggplot(diff_sites, aes(x = log2FoldChange, y = -log10(padj))) + geom_point(aes(color = padj < 0.05 & abs(log2FoldChange) > 1)) + geom_hline(yintercept = -log10(0.05), line
- exomePeak2 Differential Analysis
- QNB for Differential Methylation
- Visualization
- Related Skills
What does the bio-epitranscriptomics-m6a-differential skill do?
Identify differential m6A methylation between conditions from MeRIP-seq. Use when comparing epitranscriptomic changes between treatment groups or cell states.
How do I install it?
Run `npx skills add majiayu000/claude-skill-registry --skill m6a-differential --agent claude-code` — it drops the skill into your project so the agent can pick it up. Swap the --agent value for codex, cursor or copilot if you use one of those.
Where does this skill come from?
From majiayu000/claude-skill-registry, a repository with 534 stars. We read it straight from the repository tree rather than a submitted listing, so what you see here is what is actually published.
Is a popular skill a good skill?
Not necessarily. Stars measure attention, not adoption — a repository can trend for a week and be abandoned. That is why we show the weekly change from our own snapshots next to the total, instead of a single flattering number.