Skip to content

A tool to compute various coverage statistics given a .bed file (i.e. a CLIP or ChIP-Seq result) and a reference genome.

Notifications You must be signed in to change notification settings

alexamlie/bed_statistics

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

38 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

genewise_bed_coverage.py - compute coverage of a .bed file over all genes
entrywise_bed_coverage.py - from the other direction, given a .bed file and 4 .bed files
		         containing promoters, exons, introns, and repeat regions, calculate
		         genomic partition of each entry in the .bed file (no UTRs)
			 

UPDATE THIS

bed_gene_feature_coverage.py is the main script here. It is used to compute various coverage statistics given
a .bed file and a reference genome, which can be in any format as long as you give a header file
with the format for it. This can be generated by downloading a table from the UCSC genome browser
and then putting the top line into the header file, like:
head -n 1 ref_file.txt > ref_header.txt

The script is run as follows:

python bed_gene_feature_coverage.py reffile ref_format bedfile output log

Where reffile and ref_format are the reference and header files described above, bedfile is the
.bed file you want to compute statistics on, output is the path to the output file, and log is
the path to a log file, which contains details of peak checking for debugging purposes. 

About

A tool to compute various coverage statistics given a .bed file (i.e. a CLIP or ChIP-Seq result) and a reference genome.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published