VAriation with LOng Range information
This repository includes code for VALOR2. The original version 1 is deprecated but also available in the rasekh-etal directory for archival purposes.
VALOR2: characterization of large-scale structural variants using linked-reads.
Fatih Karaoglanoglu, Camir Ricketts, Ezgi Ebren, Marzieh Eslami Rasekh,
Iman Hajirasouliha, Can Alkan. Genome Biology, 21 (1): 72, 2020.
- gcc > 5.4
- zlib (http://www.zlib.net)
- htslib (included; http://htslib.org/)
- sonic (included; https://github.com/calkan/sonic)
git clone https://github.com/BilkentCompGen/valor.git --recursive
cd valor
make libs
make
cp valor /path/to/your/favorite/binaries
Required Parameters:
-i, --input [BAM files] : Input files in sorted BAM format.
-o, --out [output prefix] : Prefix for the output files
-s, --sonic [sonic file] : Sonic file. Check: https://github.com/calkan/sonic.
-f, --svs_to_find [sv type] : Comma separated list of SVs (i.e. INV,DUP,IDUP,TRA,ITRA,DEL)
Optional Parameters:
-l, --log_file [logfile name]: default is [output prefix]-valor.log
-t, --threads [thread count]: sets number of threads to be used (default is 1)
-p, --ploidy [Number of chromosome sets] : Default is 2
-y, --single-copy-chr [Chromosome name] : sets copy number of a chromosome to 1 (example -y X -y Y for male humans).
Help Parameters:
-v, --version : Print version and exit.
-h, --help : Print help screen and exit.
Example:
valor -i input.bam -s hg19.sonic -o output_prefix -f DUP,IDUP,DEL
SONIC files are available under https://github.com/BilkentCompGen/sonic-prebuilt/
-
human_g1k_v37.sonic: SONIC file for Human Reference Genome GRCh37 (1000 Genomes Project version)
-
ucsc_hg19.sonic: SONIC file for the human reference genome, UCSC version build hg19.
-
GRCh38.sonic: SONIC file for the human reference genome build 38.
Make sure that the same reference was used to align the reads beforehand (BAM file) and to create the SONIC file. The SONIC files and the reference FASTA files linked above are compatible.
Please refer to the SONIC development repository: https://github.com/calkan/sonic/
in [output prefix]-predicted_inversions.bedpe
Chromosome-name BP1-start BP1-end Chromosome-name BP2-start BP2-end SV_TYPE 10XG-Support Read_Pair_Support Molecule_Depth
- Read pair Support: Number of read-pairs that support these breakpoints
- 10XG Support: Number of 10XG Molecule-pairs that support these breakpoints
| chr1 | start1 | end1 | chr2 | start2 | end2 | type | Support info |
|---|---|---|---|---|---|---|---|
| chr | bp1-start | bp1-end | chr | bp2-start | bp-end | inversion | ... |
| chr | molecule-start | start | chr | end | molecule-end | deletion | |
| chr | source-start | source-end | chr | target-start | target-end | duplication | |
| source-chr | source-start | source-end | target-chr | target-start | target-end | translocation |
To build a valor Docker image
cd docker
docker build . -t valor:latest
Your image named "valor" should be ready. You can run valor using this image by
docker run --user=$UID -v /path/to/inputs:/input -v /path/to/outputdir:/output valor [args]
[args]are usual arguments you would pass to valor executable. Be careful about mapping. You need to specify folders respective to container directory structure.- You need to map host machine input and output directory to responding volume directories inside the container. These options are specified by '-v' argment.
- Docker works with root user by default. "--user" option saves your outputs.