http://www.gly.uga.edu/railsback/FieldImages.htmlMolecularMicro

PAVED- A Software suite for  the analysis of epigenome-derived next generation sequencing data

washUlogo





MAIN

INDEX


ANALYTICAL PIPELINE

CONTACT

SYSTEM REQUIREMENTS

PAVED Package  

Example Data

Tag low coverage regions

The utility filterOutLowReadDepth tags low coverage areas depending on the average coverage of that chromosome

Prerequisites

1) .depth file generated using  findReadDepth or  findFragmentDepth utility  in PAVED

2) Average coverage value per chromosome generated using findAverageCoveragePerChromosome utility in PAVED

How to run it?

Type java -jar PAVED.jar filterOutLowReadDepth -h to see list of parameters

filterOutLowReadDepth
 takes as input a .depth file and average coverage per chromosome file to tag low coverage regions

Run the utility as follows:

java -jar C:\Britta\manuscript\Analysis\PAVED.jar filterOutLowReadDepth -i C:\Britta\manuscript\Analysis\data\MNAse\normalizedData\MNAseRep2Chr5.depth -o C:\Britta\manuscript\Analysis\data\MNAse\normalizedDataMNAseRep2Chr5Filtered.depth -j  C:\Britta\manuscript\Analysis\data\MNAse\normalizedDataMNAseRep2Chr5Filtered.AvgCov -m 0.2


Here, "
C:\Britta\manuscript\Analysis\" is the location where the jar file is present on the local disk, -i is the input depth file, -j is the average coverage per chromosome file and -o is the output file with low coverage regions tagged and -m is the factor that determines low coverage regions (eg. -m 0.2 selects all regions with coverage less than 0.2*average coverage of that chromosome and tags them as low coverage regions)