The unique stacks program will take as input a set of short-read sequences and align them into exactly-matching stacks (or putative alleles). Comparing the stacks it will form a set of putative loci and detect SNPs at each locus using a maximum likelihood framework1.

1Hohenlohe PA, Bassham S, Etter PD, Stiffler N, Johnson EA, Cresko WA. 2010. Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags. PLoS Genetics, 6(2):e1000862.

Program Options

ustacks -f file_path -o path [-M max_dist] [-m min_cov] [-p num_threads]

Stack assembly options:

Gapped assembly options:

Model Options

For the SNP or Bounded SNP model:

For the Bounded SNP model:

For the Fixed model:

Example Usage

  1. Here we run ustacks independently on four samples from a genetic cross, two parents and two progeny. We specify the parameters for creating putative alleles (-m) and merging alleles into putative loci (-M). We speed up the matching process by specifying 16 parallel threads.

    % ustacks -f ./samples/f0_male.fq.gz -o ./stacks -m 3 -M 4 -p 16 % ustacks -f ./samples/f0_female.fq.gz -o ./stacks -m 3 -M 4 -p 16 % ustacks -f ./samples/progeny_01.fq.gz -o ./stacks -m 3 -M 4 -p 16 % ustacks -f ./samples/progeny_02.fq.gz -o ./stacks -m 3 -M 4 -p 16

  2. Here we run ustacks against three samples from a population. Gapped alignments will be made between alleles when forming putative loci, by default.

    % ustacks -f ./samples/sample_39-1.fq.gz -o ./stacks -M 6 -p 8 % ustacks -f ./samples/sample_40-2.fq.gz -o ./stacks -M 6 -p 8 % ustacks -f ./samples/sample_41-1.fq.gz -o ./stacks -M 6 -p 8

Other Pipeline Programs

Raw reads


Execution control

Utility programs