Now liveThe Skillselion MCP - thousands of ranked skills, loaded into your agent mid-task. No install.Get it →
gptomics avatar

Bio Epitranscriptomics M6a Peak Calling

  • 4 installs
  • 1.1k repo stars
  • Updated July 25, 2026
  • gptomics/bioskills

Call m6A peaks from MeRIP-seq IP versus input comparisons using exomePeak2 to identify modification sites.

About

Identifies m6A modification sites by calling peaks from MeRIP-seq IP versus input samples with exomePeak2. A developer uses it when locating m6A sites from methylated RNA immunoprecipitation data.

  • exomePeak2 IP-vs-input peak calling
  • Identifies m6A modification sites

Bio Epitranscriptomics M6a Peak Calling by the numbers

  • 4 all-time installs (skills.sh)
  • Ranked #1,625 of 2,064 Data Science & ML skills by installs in the Skillselion catalog
  • Data as of Aug 5, 2026 (Skillselion catalog sync)
npx skills add https://github.com/gptomics/bioskills --skill bio-epitranscriptomics-m6a-peak-calling

Add your badge

Show developers this skill is listed on Skillselion. Paste this into your README.

Listed on Skillselion
Installs4
repo stars1.1k
Last updatedJuly 25, 2026
Repositorygptomics/bioskills

What it does

Call m6A peaks from MeRIP-seq IP versus input comparisons using exomePeak2 to identify modification sites.

Files

SKILL.mdMarkdownGitHub ↗

m6A Peak Calling

exomePeak2 (Recommended)

library(exomePeak2)

# Peak calling with biological replicates
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
)

# Export peaks
exportResults(result, format = 'BED')

MACS3 Alternative

# 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

MeTPeak

library(MeTPeak)

# 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'
)

Peak Filtering

# Filter by fold enrichment and q-value
# FC > 2, q < 0.05 typical thresholds
awk '$7 > 2 && $9 < 0.05' peaks.xls > filtered_peaks.bed

Related Skills

  • merip-preprocessing - Prepare data for peak calling
  • m6a-differential - Compare peaks between conditions
  • chip-seq/peak-calling - Similar concepts

Related skills

Data Science & MLpipelinesanalytics

This week in AI coding

Five minutes, every Monday - the tools, releases and tactics for developers.

unsubscribe anytime.