The benefit of pyRAD over most alternative methods for analyzing RADseq-like data comes in its use of an alignment-clustering method (vsearch) that allows for the inclusion of indel variation which improves identification of homology across highly divergent samples. For this reason pyRAD is commonly employed for RADseq studies at deeper phylogenetic scales, however, it works equally well at shallow scales.
pyRAD is intended for use with any type of restriction-site associated DNA. It currently supports RAD, ddRAD, PE-ddRAD, GBS, PE-GBS, EzRAD, PE-EzRAD, 2B-RAD, nextRAD, and can be extended to other types. Below is the download link as well as a number of tutorials. There is a “general use” tutorial which explains how to install and setup input files, and also tutorials with example data sets for different data types and analyses. The software was initially described in the following publication (preprint, journal)
Downloads (now on github)
You can download a stable release version, or if you are comfortable with git you may clone the repository. Whichever you choose, please check for updates frequently, as I generally make bug fixes or updates weekly (see changelog).
General Use Tutorial
Recommended example Tutorial
- SE ddRAD Tutorial v.3.0.4
- PE ddRAD Tutorial v.3.0.4
- PE ddRAD w/ merged reads Tutorial v.3.0.4
- SE GBS (or ezRAD) Tutorial v.3.0.4
- PE GBS (or ezRAD) w/ merged reads Tutorial v.3.0.4
- D-statistics tutorial()
(For now, contact the google group forum for additional supported datatypes: 2B-RAD, nextRAD, etc.)