AGELESSMolecularMicro

PAINT: Leishmania Sexual Reproductive Strategies As Resolved Through Compuational Methods Designed for Aneuploid Genomes

washUlogo
Jahangheer S. Shaik, Deborah E. Dobson, David L. Sacks and Stephen M. Beverley




MAIN

INDEX


ANALYTICAL PIPELINE

CONTACT

SYSTEM REQUIREMENTS

PAINT Package  

Example Data

Generate Heatmap of ploidies

This page walks you though various steps required to generate a heatmap of ploidies using PAINT package.

Prerequisites

1) Align the fastq files to the reference genome of interest using your favourite aligner (Bowtie, BWA, Novoalign etc.)
2) Sort the generated Binary alignment map (BAM) file using samtools 'sort' function
3) Generate read depth (coverage) for each of the samples using findReadDepth utility of PAINT
4) Generate concatenated ploidy matrix using ConcatenatedPloidyMatrix utility of PAINT

How to Run it?

The R script to generate the heat map is in the directory ../DemoData/scripts/heatmap.R . The location of the concatenated ploidy matrix needs to be updated in the R script. It is marked using "#UPDATE" in the script. Also, it requires "gplots" package in R to draw the heatmap. 

Output

The output of the listAlleles program looks as follows:
heatmap
The chromosomes are the columns and samples are the rows. color scheme with histogram is depicted to the topleft. The values are overlaid on the heatmap. Edit the R script for font sizes, titles etc. as needed.