Skip to content

Latest commit

 

History

History
53 lines (48 loc) · 2.27 KB

README.md

File metadata and controls

53 lines (48 loc) · 2.27 KB

bam_to_bigwig

Convert your alignment files in BAM format into coverage files in bigWig format

Usage

bam_to_bigwig.py BAM_FILE [BAM_FILE ...]
[-o, --bigwig_filename=<output file name>
 -t, --tempfile
 -k, --keep-tempfile
 -s, --ignore-secondary
 -q, --ignore-qc-fail
 -d, --ignore-optical-pcr-duplicate
 -u, --ignore-supplementary]

--bigwig_filename: if not given will save to <BAM file prefix>.bigwig.
--tempfile: If this is set, will use tempfile library to generate file names instead of using 
            <BAM file prefix>.wig and <BAM file prefix>.sizes.
--keep-tempfile: If this is set, will not delete the intermediate files <BAM file prefix>.wig 
                 and <BAM file prefix>.sizes.
--ignore-secondary: If this is set, rsem-bam2wig will ignore alignments with the 
                    "secondary alignment" flag bit 0x100 set.
--ignore-qc-fail: If this is set, rsem-bam2wig will ignore alignments with the 
                  "not passing quality controls" flag bit 0x200 set.
--ignore-optical-pcr-duplicate: If this is set, rsem-bam2wig will ignore alignments with the 
                                "PCR or optical duplicate" flag bit 0x400 set.
--ignore-supplementary: If this is set, rsem-bam2wig will ignore alignments with the 
                        "supplementary alignment" flag bit 0x800 set.

Requirements

Installation

  • pysam
    1. Install build dependencies
    • CentOS: yum install python-devel zlib-devel
    • Ubuntu: apt-get install python-dev zlib1g-dev
    1. Install using pip
    • pip install pysam
  • wigToBigWig
    1. Download binary: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/wigToBigWig
    2. Make sure wigToBigWig is in a directory in your $PATH
    • cp wigToBigWig /usr/local/bin
  • rsem-bam2wig
    1. Download the source code: https://github.com/deweylab/RSEM
    2. Compile
    • make
    1. Make sure rsem-bam2wig is in a directory in your $PATH
    • cp rsem-bam2wig /usr/local/bin