The blasr from nine428

blasr's Introduction

Installation.

1. Requirements.

  To build BLASR, you must have hdf 1.8.0 or above installed and
  configured with c++ support (you should have the library
  libhdf5_cpp.a).  If you are intalling the entire PacBio secondary
  analysis software suite, appropriate hdf libraries are already
  distributed and no configuration is necessary.  Otherwise, it is
  necessary to point two environment variables, HDF5INCLUDEDIR and
  HDF5LIBDIR to the locations of the HDF5 libraries.

For example:

> export HDF5INCLUDEDIR=/usr/include/hdf
> export HDF5LIBDIR=/usr/lib/hdf

2. Build the source tree.

> make


3. The executables will be in cpp/alignment/bin/blasr

Running BLASR


Typing blasr -h or blasr -help on the command line will give you a
list of options.  At the least, provide a fasta,fastq, or bas.h5 file,
and a genome.  

Some typical use cases:

# Align reads from reads.bas.h5 to ecoli_K12 genome, and output in SAM
#  format.

> blasr reads.bas.h5  ecoli_K12.fasta -sam


# Same as above, but with soft clipping

> blasr reads.bas.h5  ecoli_K12.fasta -sam -clipping soft

# Use multiple threads

> blasr reads.bas.h5  ecoli_K12.fasta -sam -clipping soft -out  alignments.sam -nproc 16

# Include a larger minimal match, for faster but less sensitive
# alignments

> blasr reads.bas.h5  ecoli_K12.fasta -sam -clipping soft -minMatch 15


# Produce alignments in a pairwise human readable format

> blasr reads.bas.h5  ecoli_K12.fasta -m 0


# Use a precomputed suffix array for faster startup

# first precompute the suffix array
> sawriter hg19.fasta.sa hg19.fasta  
> blasr reads.bas.h5 hg19.fasta -sa hg19.fasta.sa


# Use a precomputed BWT-FM index for smaller runtime memory footprint,
# but slower alignments.
> sa2bwt hg19.fasta hg19.fasta.sa hg19.fasta.bwt
> blasr reads.bas.h5 hg19.fasta -bwt hg19.fasta.bwt

Recommend Projects