8.4. binReads

8.4.1. Purpose

For a single channel of input, creates wiggle files showing binned coverage, which can be visualized using the UCSC Genome Browser, IGV and other browsers. See http://genome.ucsc.edu/goldenPath/help/wiggle.html.

For multiple channels, generates tables of binned coverage which cannot be visualized.

8.4.2. Parameters

usage: usage="binReads [options] <alignment_file(s)>"

options:
  -h, --help            show this help message and exit
  -b BINSIZE, --binsize=BINSIZE
                        binsize [default=1000]
  -o OUTFILE, --outfile=OUTFILE
                        output filename [default=<uses base filename>]
  -v, --verbose         print debugging output [default=False]
  -s, --sort            sort files first [default=False]
  -z, --gzip            gzip resulting file [defualt=False]
  -m REFMAP, --refMap=REFMAP
                        transform refIDs to refnames in the map file (e.g. the
                        output of align2txt -map <file>)
  --merge               create a single binning based on multiple alignment
                        files
  --aggregate           creates a tab-deliminated (e.g. not WIGGLE) formated
                        tally with count columns for each input file for each
                        reference bin
  --zero                include bins with a zero count
  --minCount=MINCOUNT   report only bins with counts >= minCount  [default=1]
  --bufferSize=MEM      Buffer Size to allocate when sorting text files
                        (default=1G)
  -n TRACK_NAME, --track_name=TRACK_NAME
                        name to display on browser track
  -r, --reads_per_million
                        report bin counts in reads/million units
                        (default=False)