draw.target.coverage.plot: Draw target coverage plot

In ProfessionalFarmer/loonR: Common function in bioinformatics analysis

Draw target coverage plot

Description

Draw target coverage plot

Usage

draw.target.coverage.plot(dir)

Arguments

dir	bedtools coverage output file should be in the same directory. File should suffix with "hist.all.txt$"

Examples

Clone from https://gist.github.com/stephenturner/9396409   Thanks
input file from: bedtools coverage -hist -b  samp.01.bam -a target_regions.bed | grep ^all > samp.01.bam.hist.all.txt  or bedtools genomecov -ibam ../merge.bam -g genome.file.txt | grep ^genome > genome.cov
bedtools -g and -sorted should be noticed
Ref: http://www.gettinggeneticsdone.com/2014/03/visualize-coverage-exome-targeted-ngs-bedtools.html

ProfessionalFarmer/loonR documentation built on Oct. 9, 2024, 9:56 p.m.