Skip to content

mcvickerlab/WASP2

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

391 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

WASP2 - Allele-Specific Analysis Pipeline

PyPI Bioconda CI Documentation

Installation

Bioconda

mamba install -c conda-forge -c bioconda wasp2

PyPI

pip install wasp2

The PyPI package does not install external tools such as samtools, bcftools, or bedtools; install those separately.

Docker

docker pull ghcr.io/mcvickerlab/wasp2:1.4.0
docker run --rm ghcr.io/mcvickerlab/wasp2:1.4.0 wasp2-count --help

Singularity/Apptainer

singularity pull wasp2.sif docker://ghcr.io/mcvickerlab/wasp2:1.4.0
singularity exec wasp2.sif wasp2-count --help

CLI Tools

WASP2 installs four command-line entry points:

  • wasp2-map
  • wasp2-count
  • wasp2-analyze
  • wasp2-ipscore

Quick Start

1. Correct mapping bias

wasp2-map make-reads input.bam variants.vcf.gz -s sample1 -o remap_dir

# Realign remap_dir/*_swapped_alleles_r1.fq and r2.fq with the same aligner
# and settings used for the original BAM, then:

wasp2-map filter-remapped remapped.bam \
  -j remap_dir/input_wasp_data_files.json \
  -o filtered.bam

2. Count alleles

wasp2-count count-variants filtered.bam variants.vcf.gz -s sample1 -o counts.tsv

3. Test for imbalance

wasp2-analyze find-imbalance counts.tsv -o ai_results.tsv

Single-Cell Example

wasp2-count count-variants-sc \
  cellranger.bam \
  variants.vcf.gz \
  barcodes.tsv \
  --samples sample1 \
  --feature genes.gtf \
  --out_file allele_counts.h5ad

wasp2-analyze find-imbalance-sc \
  allele_counts.h5ad \
  barcode_groups.tsv \
  --sample sample1 \
  -o ai_results.tsv

iPSCORE Utilities

wasp2-ipscore inventory --output inventory.tsv
wasp2-ipscore manifest --output manifest.csv
wasp2-ipscore validate

See the documentation for complete usage, tutorials, and API details.

About

WASP2: High-performance allele-specific analysis pipeline for next-generation sequencing data

Topics

Resources

License

Code of conduct

Contributing

Stars

Watchers

Forks

Packages

 
 
 

Contributors