bio-long-read-sequencing-nanopore-methylation — independently scanned and version-tracked by SaferSkills.
SaferSkills independently audited bio-long-read-sequencing-nanopore-methylation (Agent Skill) and scored it 100/100 (green). The audit ran 55 deterministic rules across Security, Supply Chain, Maintenance, Transparency, and Community; it found 0 high-severity and 0 lower-severity findings. The full rule-by-rule trace and per-finding evidence are below. Free, methodology-open.
Findings & checks · 0 flagged
Every scanned point with the score it earned and what moved between them.
First recorded scan — no prior version to compare against.
The primary manifest — the file an agent reads to learn what this artifact does.
Reference examples tested with: methylKit 1.28+, minimap2 2.26+, samtools 1.19+
Before using code patterns, verify installed versions match. If versions differ:
<tool> --version then <tool> --help to confirm flagsIf code throws ImportError, AttributeError, or TypeError, introspect the installed package and adapt the example to match the actual API rather than retrying.
"Call methylation from my Nanopore reads" → Extract 5mC/6mA modification probabilities from basecalled reads and summarize per-site methylation frequencies.
modkit pileup aligned.bam methylation.bed --ref ref.faONT's modkit is the recommended tool for methylation analysis from basecalled data.
# Assumes BAM has MM/ML tags from dorado basecalling
modkit pileup input.bam methylation.bed \
--ref reference.fa \
--cpg \
--combine-strands# bedMethyl format
chr1 1000 1001 . 10 + 1000 1001 0,0,0 10 80.5
# Columns: chrom, start, end, name, score, strand, thickStart, thickEnd,
# itemRgb, coverage, percent_modified# Dorado basecalling with 5mC model
dorado basecaller [email protected] \
pod5_dir/ \
--modified-bases 5mCG \
> calls.bam
# Index and align
samtools fastq calls.bam | \
minimap2 -ax map-ont -y reference.fa - | \
samtools sort -o aligned.bam
samtools index aligned.bam# CpG islands only
modkit pileup aligned.bam cpg_islands.bed \
--ref reference.fa \
--cpg \
--include-bed cpg_islands.bed
# Promoter regions
modkit pileup aligned.bam promoters.bed \
--ref reference.fa \
--cpg \
--include-bed promoters.bed# Get modification summary statistics
modkit summary aligned.bam
# Output includes:
# - Total reads with modifications
# - Modification types detected
# - Fraction modified per type# Create BED files for each sample
modkit pileup sample1.bam sample1.bed --ref ref.fa --cpg
modkit pileup sample2.bam sample2.bed --ref ref.fa --cpg
# Compare with methylKit or DSS in R~30 seconds. Free. No account. Every finding cites a rule and a line of evidence.