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 Bottle Brush Plots

This page walks you though various steps required to generate bottle brush plots using PAINT package.

Pipeline

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 an allele files using findAlleles utility of PAINT
4) Generate an list files using listAlleles utility of PAINT
5) Find parental SNVs using findSNPsSomy utility of PAINT. Extract Homozygous SNVs from the parental lines and place them in a directory.
6) Generate merged allele file using findAllelesSNPpositions utility of PAINT
7) Generate merged allele frequency file using findAlleleFreqsSNPpositions utility of PAINT
8) Divide the merged allele file into individual files using DivideMergedFiles2ChrFiles utility of PAINT.
9) Divide the merged allele frequency file into individual files using DivideMergedFiles2ChrFiles utility of PAINT.
10) Find parental allele frequencies using getParentAllelFrequencies utility of PAINT.
11) Find coverage using findReadDepth utility of PAINT
12) Divide the coverage files into individual chromosome coverage files using divideReadDepthFile utility of PAINT
13) Update all file locations and other information where it says " #UPDATE" in the script../DemoData/Scripts/DrawParFreqPloidyEachLocii.R
14) Use the script ../DemoData/Scripts/DrawParFreqPloidyEachLocii.R to generate the bottle brush plots.

Sample Output

bottle brush plot