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

Bio Epitranscriptomics Modification Visualization

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

Create metagene plots and browser tracks for RNA modification data using Guitar to visualize m6A distribution around genomic features.

About

Generates metagene plots and browser tracks for RNA modification data with the Guitar package. A developer uses it when visualizing m6A distribution patterns around features like stop codons.

  • Guitar metagene plots for RNA modifications
  • Browser tracks showing m6A distribution

Bio Epitranscriptomics Modification Visualization 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-modification-visualization

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

Create metagene plots and browser tracks for RNA modification data using Guitar to visualize m6A distribution around genomic features.

Files

SKILL.mdMarkdownGitHub ↗

Modification Visualization

Metagene Plots with Guitar

library(Guitar)
library(TxDb.Hsapiens.UCSC.hg38.knownGene)

# Load m6A peaks
peaks <- import('m6a_peaks.bed')

# Create metagene plot
# Shows distribution relative to transcript features
GuitarPlot(
    peaks,
    txdb = TxDb.Hsapiens.UCSC.hg38.knownGene,
    saveToPDFprefix = 'm6a_metagene'
)

Custom Metagene with deepTools

# Create bigWig from IP/Input ratio
bamCompare -b1 IP.bam -b2 Input.bam \
    --scaleFactors 1:1 \
    --ratio log2 \
    -o IP_over_Input.bw

# Metagene around stop codons
computeMatrix scale-regions \
    -S IP_over_Input.bw \
    -R genes.bed \
    --regionBodyLength 2000 \
    -a 500 -b 500 \
    -o matrix.gz

plotProfile -m matrix.gz -o metagene.pdf

Browser Tracks

# Create normalized bigWig for genome browser
bamCoverage -b IP.bam \
    --normalizeUsing CPM \
    -o IP_normalized.bw

# Peak BED to bigBed
bedToBigBed m6a_peaks.bed chrom.sizes m6a_peaks.bb

Heatmaps

library(ComplexHeatmap)

# m6A signal around peaks
Heatmap(
    signal_matrix,
    name = 'm6A signal',
    cluster_rows = TRUE,
    show_row_names = FALSE
)

Related Skills

  • epitranscriptomics/m6a-peak-calling - Generate peaks for visualization
  • data-visualization/genome-tracks - IGV, UCSC integration
  • chip-seq/chipseq-visualization - Similar techniques

Related skills

This week in AI coding

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

unsubscribe anytime.