For studies using high-throughput sequencing (HTS) data, simulations can be vital for planning sampling design and testing bioinformatic tools. However, most HTS sequencing tools provide only very simple ways of adding deviations from a reference genome. For HTS studies that focus on patterns of genomic variation among individuals, populations, or species, having a tool that can simulate realistic patterns of molecular evolution and generate HTS data from those simulations would be quite useful.
jackalope simply and efficiently simulates (i) haplotypes from reference genomes and (ii) reads from both Illumina and Pacific Biosciences (PacBio) platforms. It can either read reference genomes from FASTA files or simulate new ones. Variant haplotypes can be simulated using summary statistics, phylogenies, Variant Call Format (VCF) files, and coalescent simulationsâthe latter of which can include selection, recombination, and demographic fluctuations. jackalope can simulate single, paired-end, or mate-pair Illumina reads, as well as reads from Pacific Biosciences These simulations include sequencing errors, mapping qualities, multiplexing, and optical/PCR duplicates. All outputs can be written to standard file formats.
Before installing jackalope, you should update the packages Rhtslib and zlibbioc. Since both of these are on Bioconductor, you should update BiocManager, too.
To use multithreading in jackalope, youâll need to compile it from source using the proper flags. If youâve enabled OpenMP properly, running jackalope:::using_openmp() in R should return TRUE.
The first step is to add the following to the .R/Makevars (.R/Makevars.win on Windows) file inside the home directory:
Then, you should be able to install jackalope by running the following in R:
Follow the directions here to install R compiler tools: https://thecoatlessprofessor.com/programming/cpp/r-compiler-tools-for-rcpp-on-macos/.
Check your version of gcc using gcc --version in the Terminal. Then, check the table at https://mac.r-project.org/openmp/ to see which version of the runtime OpenMP downloads you need. For LLVM version 9.0.1, you run the following in the Terminal:
For the next step of actually installing jackalope, one option is to add the following to your ~/.R/Makevars file:
⦠then installing jackalope by running install.packages("jackalope", type = "source") or remotes::install_github("lucasnell/jackalope") in R.
This might not be desirable since it affects all package installations. An alternative method is to use the package withr:
Add the following to the .R/Makevars file inside the home directory:
Next, go to https://cran.r-project.org/bin/macosx/tools/ and download the newest versions of (1) the clang compiler (version 8 at the time of writing) and (2) GNU Fortran (version 6.1 at the time of writing). The downloads will have the .pkg extension. Next, install clang and gfortran by opening these .pkg files and following the directions.
After this, add the following to your ~/.R/Makevars file (replacing clang8 with your version of the clang compiler):
Now you should be able to install jackalope by running install.packages("jackalope", type = "source") in R.
For more information, please see https://thecoatlessprofessor.com/programming/openmp-in-r-on-os-x/.
UsageBelow shows how to simulate a 10kb genome, then create haplotypes from that genome using a phylogenetic tree:
Below simulates 500 million paired-end, 100 bp reads from the haplotypes:
Below simulates 500 thousand PacBio reads from the reference genome:
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4