# HG changeset patch # User nikos # Date 1409673604 14400 # Node ID eebdd6da935398810a14516a3fdf0c35b2efcc0b # Parent 558f7f9baa78b3bb34796c4f687a65610b98554b Deleted selected files diff -r 558f7f9baa78 -r eebdd6da9353 bedClip.sh --- a/bedClip.sh Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -#!/bin/bash - -bedClip $1 <(fetchChromSizes $2 2> /dev/null) $3 diff -r 558f7f9baa78 -r eebdd6da9353 bedClip.xml --- a/bedClip.xml Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,53 +0,0 @@ - - Remove lines from bed file that refer to off-chromosome places. - - ## Set genome assembly - - #set $Genome = str( $genome_cond.genome ) - #if str( $genome_cond ) == 'OTHER': - #set $Genome = str( $genome_cond.genome_other ) - #end if - - bedClip.sh $input $Genome $output - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -**What it does** -bedClip - Remove lines from bed file that refer to off-chromosome places. - -**Usage** - - bedClip input.bed chrom.sizes output.bed - - - - diff -r 558f7f9baa78 -r eebdd6da9353 bedExtendRanges.xml --- a/bedExtendRanges.xml Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,68 +0,0 @@ - - Extend length of entries in bed 6+. - - ## Set genome assembly - - #set $Genome = str( $genome_cond.genome ) - #if str( $genome_cond ) == 'OTHER': - #set $Genome = str( $genome_cond.genome_other ) - #end if - - bedExtendRanges -user=genome -host=genome-mysql.cse.ucsc.edu $Genome $length $input 2> /dev/null > $output - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -**What it does** - -bedExtendRanges - extend length of entries in bed 6+ data to be at least the given length, -taking strand directionality into account. - -**Usage** - - bedExtendRanges database length files(s) - -**Example** - - * bedExtendRanges -user=genome -host=genome-mysql.cse.ucsc.edu hg18 250 stdin - - will transform: - chr1 500 525 . 100 + - chr1 1000 1025 . 100 - - to: - chr1 500 750 . 100 + - chr1 775 1025 . 100 - - - - - diff -r 558f7f9baa78 -r eebdd6da9353 bigWigSummary.xml --- a/bigWigSummary.xml Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,33 +0,0 @@ - - Extract summary information from a bigWig file. - - bigWigSummary $input $chrom $start $end $dataPoints -type=$type> $output - - - - - - - - - - - - - - - - - - - - - - -**Usage** - * bigWigSummary file.bigWig chrom start end dataPoints - * Get summary data from bigWig for indicated region, broken into dataPoints equal parts. (Use dataPoints=1 for simple summary.) - - - - diff -r 558f7f9baa78 -r eebdd6da9353 bigwig2summary.xml --- a/bigwig2summary.xml Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,85 +0,0 @@ - - Extract summary information from a bigWig file across multiple genomic regions specified by the user. - - bigWigSummary - bigwig2summary_stdout - bedClip - - - /home/jens/bin/bigwig2summary_stdout -f $input_bw -b $input_bed -n $bins -o $assembly -t $type $header > $output - - - - - - - - - - - - - - - - - - - - - -This tool extracts summary values (mean, min, max, std or coverage) from a **bigWig** file for a number of equal sized bins across genomic regions given in an a "BED-like" file. - -The script this tool is based on is written by Jens Vilstrup Johansen and uses bigWigSummary_ and bedClip_. - -.. _bigWigSummary: https://github.com/adamlabadorf/ucsc_tools/blob/master/executables/bigWigSummary - -.. _bedClip: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bedClip - ------ - -.. class:: infomark - -The file contaning the genomic region must be TAB-delimited with at list 3 columns representing Chromosome, ChrStart and ChrEnd. If the file contains 6 columns (or more), column 6 is expected to contain strand information. Summary values from a negative strand will be reversed. - ------ - -**Example 1** - -Input BED file:: - - chr19 50178708 50180708 - chr6 90348174 90350174 - chr16 58495848 58497848 - chr5 180580242 180582242 - chr9 120177017 120179017 - -Extract summary (*#* of bins = 3):: - - 0 0 0 - 0.144886 0 0 - 0.507327 1.14649 1.38456 - 0.221471 0.144886 0.309857 - 0.348944 0.426638 0.244495 - -**Example 2** - -Input BED file (with strand information):: - - chr19 50178708 50180708 NM_198318 0 + PRMT1 - chr6 90348174 90350174 NM_020466 0 - LYRM2 - chr16 58495848 58497848 NM_020465 0 + NDRG4 - chr5 180580242 180582242 NM_206880 0 + OR2V2 - chr9 120177017 120179017 NM_014010 0 - ASTN2 - -Extract summary (*#* of bins = 3):: - - 0 0 0 - 0 0 0.144886 - 0.507327 1.14649 1.38456 - 0.221471 0.144886 0.309857 - 0.244495 0.426638 0.348944 - - - - diff -r 558f7f9baa78 -r eebdd6da9353 make_ucsc_track.sh --- a/make_ucsc_track.sh Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,59 +0,0 @@ -#!/bin/bash - -# This tool creates a UCSC track line given a bigWig or bedGraph input file, a track name and genome assembly. - -## Define URL prefix - -URL="http://data.bric.dk/galaxy/UCSC" - -#parse input -while getopts hi:t:n:o: myarg -do case "$myarg" in - h) echo "Usage: create UCSC track line -i -t -n -o " - exit ;; - i) input_file="$OPTARG" ;; - t) file_type="$OPTARG" ;; - n) track_name="$OPTARG" ;; - o) org_assembly="$OPTARG" ;; - [?]) echo "Usage: create UCSC track line -i -t -n -o " - - exit 1 ;; - esac -done - -#################################################################### -# Check if a link to that file exists in the UCSC webfolder -flag=0 - -for link in `ls -1 $UCSCPATH` -do - if [ `readlink -f $input_file` == `readlink -f $UCSCPATH/$link` ] ; then - link_path=/$link # use already existing link - flag=1 - fi -done - -# If link to input file doesn't exist, create anonymous soft-link for input file in a folder accessible from the web. - -if [ $flag == 0 ]; then - - link_path=`mktemp -u --tmpdir=/` - ln -s $input_file $UCSCPATH$link_path -fi - -#################################################################### -##Convert 4th character of file_type to uppercase (e.g. bigwig -> bigWig). UCSC browser doesn't recognize the file type otherwise. - -##Split file_type string to prefix (big) and suffix (wig, bed) and uppercase the first character of the latter. -prefix=`echo ${file_type:0:3}` -suffix=`echo ${file_type:3} | sed -e "s/\b\(.\)/\u\1/g"` - -file_type=`echo ${prefix}${suffix}` - -#################################################################### -#Contruct URL -bigURL=`echo ${URL}${link_path}` - -#################################################################### -## Print track line -echo "track name='$track_name' db=$org_assembly type=$file_type visibility=full alwaysZero=ON bigDataUrl=$bigURL" diff -r 558f7f9baa78 -r eebdd6da9353 make_ucsc_track.xml --- a/make_ucsc_track.xml Tue Jul 15 07:19:46 2014 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,34 +0,0 @@ - - Create a track line to import a bigWig or bigBed file in the UCSC genome browser - - make_ucsc_track.sh -i $input -t ${input.ext} -n '$name' -o $assembly > $output - - - - - - - - - - - - - - - - - - -This tool creates a UCSC track line for a bigWig or bigBed dataset in Galaxy. The text output can be pasted in the "Paste URLs or data" form in the `Add Custom Tracks`_ section of the UCSC Genome Browser. - -.. _Add Custom Tracks: http://genome.ucsc.edu/cgi-bin/hgCustom - - - -