// Run SnpEff, a variant annotation and effect prediction tool for genomic variants. Use when the user wants to annotate VCF files, build genome databases, or use SnpEff utilities.
| name | snpeff |
| description | Run SnpEff, a variant annotation and effect prediction tool for genomic variants. Use when the user wants to annotate VCF files, build genome databases, or use SnpEff utilities. |
| allowed-tools | Bash, Read, Grep, Glob |
SnpEff is a variant annotation and effect prediction tool. It analyzes VCF files and predicts the functional effects of genetic variants (such as amino acid changes) on known genes.
Full documentation: https://pcingola.github.io/SnpEff/
Use the wrapper script at .claude/skills/snpeff/snpeff.sh. It handles JVM memory defaults and argument passthrough. The JAR file is expected at $HOME/snpEff/snpEff.jar.
IMPORTANT: SnpEff output is typically very large (thousands to millions of lines). ALWAYS redirect output to a file. NEVER let output print to stdout, as it will fill the context window and make the conversation unusable.
# Correct: redirect to file
.claude/skills/snpeff/snpeff.sh <command> [options] [arguments] > output.vcf 2> snpeff.log
# WRONG: never do this
.claude/skills/snpeff/snpeff.sh <command> [options] [arguments]
| Command | Description | Docs |
|---|---|---|
ann / eff | Annotate variants (default command) | commandline, running, input/output |
build | Build a SnpEff database from reference genome files | build_db, build_db_gff_gtf |
buildNextProt | Build NextProt database from XML files | commands |
cds | Compare CDS sequences (database check) | commands |
closest | Annotate closest genomic region | commands |
count | Count reads/bases overlapping genomic intervals | commands |
databases | List available databases | commands |
download | Download a pre-built database | commands |
dump | Dump database contents (BED/TXT) | commands |
genes2bed | Create BED file from gene list | commands |
len | Calculate genomic length per marker type | commands |
pdb | Build interaction database from PDB/AlphaFold data | build_pdb |
protein | Compare protein sequences (database check) | commands |
seq | Translate DNA sequence to protein | commands |
show | Show gene/transcript text representation | commands |
translocReport | Create translocation report with SVG | commands |
| Topic | File |
|---|---|
| Introduction | introduction |
| Additional annotations | additionalann |
| Building regulation databases | build_reg |
| Cancer samples | cansersamples |
| FAQ | faq |
| Human genomes | human_genomes |
| Integration (GATK, Galaxy) | integration |
| Output summary | outputsummary |
| Troubleshooting | troubleshooting |
| Download & Install | download |
| Examples | examples |
| Help | help |