bio-epitranscriptomics-m6a-peak-calling
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
More from gptomics/bioskills
bioskills
Installs 425 bioinformatics skills covering sequence analysis, RNA-seq, single-cell, variant calling, metagenomics, structural biology, and 56 more categories. Use when setting up bioinformatics capabilities or when a bioinformatics task requires specialized skills not yet installed.
100bio-single-cell-batch-integration
Integrate multiple scRNA-seq samples/batches using Harmony, scVI, Seurat anchors, and fastMNN. Remove technical variation while preserving biological differences. Use when integrating multiple scRNA-seq batches or datasets.
5bio-epitranscriptomics-merip-preprocessing
Align and QC MeRIP-seq IP and input samples for m6A analysis. Use when preparing MeRIP-seq data for peak calling or differential methylation analysis.
5bio-data-visualization-multipanel-figures
Combine multiple plots into publication-ready multi-panel figures using patchwork, cowplot, or matplotlib GridSpec with shared legends and panel labels. Use when combining multiple plots into publication figures.
5bio-data-visualization-specialized-omics-plots
Reusable plotting functions for common omics visualizations. Custom ggplot2/matplotlib implementations of volcano, MA, PCA, enrichment dotplots, boxplots, and survival curves. Use when creating volcano, MA, or enrichment plots.
5bio-read-qc-fastp-workflow
All-in-one read preprocessing with fastp including adapter trimming, quality filtering, deduplication, base correction, and HTML report generation. Use when preprocessing Illumina data and wanting a single fast tool instead of separate Cutadapt, Trimmomatic, and FastQC steps.
5