An amplicon sequencing based translocation detection tool from NGS.
16
Detect chromosomal translocations from amplicon sequencing data — packaged as a ready-to-run Docker image. No local Python, conda, Java, or samtools installation required. Everything the tool needs is baked into the image.
PASTA (Primer Anchored Statistical Translocation Analysis) takes demultiplexed amplicon sequencing reads (uBAM files) and:
It is designed for primer-anchored assays where P5 and P7 indices are both fixed or random.
Note: "Amplicons" here means two primers per target site. PASTA does not support current AMP-seq technology.
docker pull idtdna/pasta:v1.1.0
The image's entrypoint is the PASTA orchestrator, so you can pass --help directly:
docker run --rm idtdna/pasta:v1.1.0 --help
PASTA reads all sample metadata and file paths from a single YAML configuration file. Because your data (uBAMs, reference genome, BED files) and the config live on your machine, you must mount them into the container and make sure the paths inside the config match the mounted locations.
docker run --rm \
-v /path/to/your/data:/data \
-v /path/to/results:/results \
idtdna/pasta:v1.1.0 \
--config /data/pasta_config.yaml \
--threads 8 \
--tool-threads 4 \
--outdir /results
-v /path/to/your/data:/data — mounts your local data folder (config, uBAMs, genome, BED
files) into the container at /data.-v /path/to/results:/results — mounts a local folder where outputs will be written.Important: The file paths inside your YAML config (
uBAM,Amplicon_Sequences,Genome.hg38.Sequence,TargetBed, etc.) must be the paths as seen inside the container (e.g./data/...), not the paths on your host machine.
| Option | Default | Description |
|---|---|---|
--config | (required) | Path to the YAML configuration file describing samples and parameters. |
--threads | 8 | Number of concurrent tasks (workers) the orchestrator runs in parallel. |
--tool-threads | 32 | Number of threads handed to each underlying tool (e.g. primer identification). |
--outdir | results/ | Root directory for all outputs. |
--skip-primer-id | (off) | Skip primer approximation/identification and reuse existing *_OUTPUT_PRIMERS_ID.txt files. |
--primers-dir | (none) | Directory of precomputed *_OUTPUT_PRIMERS_ID.txt files (only used with --skip-primer-id). |
Performance guidance: PASTA manages the whole workflow within a single process (no external
workflow manager). For best results use one large machine and scale --threads cautiously —
some steps are memory intensive and too many workers can exhaust available memory.
Create a YAML file (for example pasta_config.yaml) describing each sample. A minimal example:
Samples:
Sample_Control:
uBAM: /data/uBAM/control.unmapped.bam
ComparisonGroup: ['1'] # groups this sample participates in
SampleType: 'Control' # 'Control' or 'Treatment'
Description: Control_Rep1
Amplicon_Sequences: /data/primers/amplicons_6col_amp.bed
pasta_downsample: TRUE # 'TRUE'/'FALSE' — enable per-sample downsampling
pasta_downsample_value: 100000 # reads to downsample to when enabled
pasta_pvalue: 0.01 # significance threshold
pasta_signal: 0.01 # minimum % translocation signal
pasta_method: 'deseq2' # 'hypergeom' or 'deseq2'
pasta_min_cov: 100 # minimum coverage before flagging
pasta_calc: 'mean' # translocation frequency calc method
SummaryGraphs: true # optional — render summary graphs
Genome:
hg38:
Sequence: /data/genome/hg38_primary.fa
CustomGenomeFasta: /data/genome/hg38_primary.fa
TargetBed: /data/primers/amplicons.bed
SnpEffGenome: hg38
Sample_Treatment:
uBAM: /data/uBAM/treatment.unmapped.bam
ComparisonGroup: ['1']
SampleType: 'Treatment'
Description: Treatment_Rep1
Amplicon_Sequences: /data/primers/amplicons_6col_amp.bed
pasta_method: 'hypergeom'
pasta_pvalue: 0.01
pasta_signal: 0.01
pasta_min_cov: 100
pasta_calc: 'mean'
Genome:
hg38:
Sequence: /data/genome/hg38_primary.fa
CustomGenomeFasta: /data/genome/hg38_primary.fa
TargetBed: /data/primers/amplicons.bed
SnpEffGenome: hg38
Requirements for a valid comparison: each ComparisonGroup must contain at least one
Control and one Treatment sample. Groups that don't meet this are skipped.
| Parameter | Meaning |
|---|---|
SampleType | Control or Treatment. |
ComparisonGroup | List of group IDs this sample belongs to. |
uBAM | Path (inside the container) to the demultiplexed unmapped BAM. |
Amplicon_Sequences / TargetBed | BED file describing the expected amplicon bounds. |
Genome.hg38.Sequence | Reference genome FASTA used to reconstruct primers. |
pasta_method | deseq2 or hypergeom. |
pasta_pvalue | p-value significance threshold. |
pasta_signal | Minimum % translocation signal to keep a call. |
pasta_min_cov | Minimum coverage before a call is flagged. |
pasta_downsample / pasta_downsample_value | Enable and set per-sample read downsampling. |
Point --outdir at a mounted folder. PASTA writes:
results/
├── bamfiles/ # BAMs produced during primer identification
├── primers/ # Primer inputs, identified primers, and amplicon files (troubleshooting)
├── group_translocations/ # Per-comparison translocation CSVs + summaries (key deliverable)
├── sample_translocations/ # Legacy directory (slated for deprecation)
├── graphs/ # Summary graphs (when SummaryGraphs is enabled)
└── all_identified_translocations.csv # Final aggregated results (also printed to the console)
# Layout on the host:
# /home/me/pasta_run/
# ├── pasta_config.yaml
# ├── uBAM/ (your *.unmapped.bam files)
# ├── genome/ (hg38 FASTA)
# └── primers/ (amplicon BED files)
docker run --rm \
-v /home/me/pasta_run:/data \
-v /home/me/pasta_run/results:/results \
idtdna/pasta:latest \
--config /data/pasta_config.yaml \
--threads 8 \
--tool-threads 4 \
--outdir /results
When it finishes, your calls are in /home/me/pasta_run/results/ and the aggregated table is
printed to the terminal.
/data/...).--threads / --tool-threads and scale up as resources allow.--user "$(id -u):$(id -g)" to the docker run command if you want output files owned by
your host user instead of root.PASTA is made available for academic, educational, governmental, non-profit, and internal commercial research purposes.
Commercial organizations may use PASTA for internal research and development activities, including the discovery, development, testing, and evaluation of therapeutic, agricultural, industrial, and research applications.
Without obtaining a separate commercial license from the copyright holder, PASTA may not be used by any person or organization whose primary business includes the commercial development, manufacture, marketing, sale, licensing, or provision of genome editing products, reagents, platforms, software, or services.
This restriction includes, but is not limited to, entities that:
The following are examples of permitted uses under this license:
For the avoidance of doubt, the development of therapeutic products, diagnostic products, agricultural products, or other end-use applications is considered a permitted use under this license, provided that the user is not commercially offering genome editing reagents, genome editing platforms, genome editing software products, or genome editing analysis services to third parties.
Organizations developing drugs, biologics, cell therapies, gene therapies, or other therapeutic products for their own pipelines may use PASTA for internal research and development activities, even when those activities involve genome editing technologies.
The restrictions in this license are intended to apply primarily to organizations whose principal commercial activities involve the sale, licensing, manufacture, or provision of genome editing reagents, genome editing products, genome editing software, or genome editing services, and are not intended to restrict the internal use of PASTA by therapeutic developers, pharmaceutical companies, biotechnology companies, or similar organizations developing end-use products.
By downloading, accessing, or using PASTA, you agree to the terms described above. Organizations wishing to engage in a restricted commercial use must obtain a separate written commercial license from the copyright holder before using PASTA for such purposes.
Content type
Image
Digest
sha256:43705e278…
Size
1.4 GB
Last updated
3 months ago
docker pull idtdna/pasta:v1.1.0