Data & ReproBioinformatics & GenomicsFreedomIntelligence/OpenClaw-Medical-SkillsData & Reproduction

bio-epitranscriptomics-m6a-peak-calling

Maintainer FreedomIntelligence · Last updated April 1, 2026

Call m6A peaks from MeRIP-seq with exomePeak2 or MACS2.

OpenClawNanoClawAnalysisReproductionbio-epitranscriptomics-m6a-peak-calling🧠 bioos extended suitebioos extended bioinformatics suitecall

Original source

FreedomIntelligence/OpenClaw-Medical-Skills

https://github.com/FreedomIntelligence/OpenClaw-Medical-Skills/tree/main/skills/bio-epitranscriptomics-m6a-peak-calling

Maintainer: FreedomIntelligence
License: MIT
Last updated: April 1, 2026

Skill Snapshot

Key Details From SKILL.md

2 min

Key Notes

result <- exomePeak2( bam_ip = c('IP_rep1.bam', 'IP_rep2.bam'), bam_input = c('Input_rep1.bam', 'Input_rep2.bam'), gff = 'genes.gtf', genome = 'hg38', paired_end = TRUE ).
exportResults(result, format = 'BED').

Source Doc

Excerpt From SKILL.md

exomePeak2 (Recommended)

library(exomePeak2)

## Call peaks treating input as control

macs3 callpeak \
    -t IP_rep1.bam IP_rep2.bam \
    -c Input_rep1.bam Input_rep2.bam \
    -f BAMPE \
    -g hs \
    -n m6a_peaks \
    --nomodel \
    --extsize 150 \
    -q 0.05

GTF-aware peak calling

metpeak( IP_BAM = c('IP_rep1.bam', 'IP_rep2.bam'), INPUT_BAM = c('Input_rep1.bam', 'Input_rep2.bam'), GENE_ANNO_GTF = 'genes.gtf', OUTPUT_DIR = 'metpeak_output' )

Use cases

Use bio-epitranscriptomics-m6a-peak-calling for genomics and bioinformatics workflows.
Apply bio-epitranscriptomics-m6a-peak-calling to sequencing, variant, or omics analysis tasks.