# HG changeset patch
# User nikos
# Date 1410272881 14400
# Node ID a47f03c50780abbd70f63e5ff8569d7a91a2704c
# Parent 0847315fa035c7ef942aea68e491a9bd9bd00c1b
Uploaded
diff -r 0847315fa035 -r a47f03c50780 bedClip.sh
--- a/bedClip.sh Tue Sep 09 10:22:40 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 0847315fa035 -r a47f03c50780 bedClip.xml
--- a/bedClip.xml Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,65 +0,0 @@
-
- Remove lines from bed file that refer to off-chromosome places.
-
- fetchChromSizes
- bedClip
-
-
- ## 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 0847315fa035 -r a47f03c50780 bedExtendRanges.xml
--- a/bedExtendRanges.xml Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,82 +0,0 @@
-
- Extend length of entries in bed 6+.
-
- bedExtendRanges
-
-
- ## 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.
-
-.. _bedExtendRanges: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bedExtendRanges
-
-**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 0847315fa035 -r a47f03c50780 bigWigSummary.xml
--- a/bigWigSummary.xml Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,43 +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 0847315fa035 -r a47f03c50780 bigwig2summary.sh
--- a/bigwig2summary.sh Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,131 +0,0 @@
-#!/bin/bash
-# $Id: bigwig2summary_stdout 34 2014-02-20 08:31:20Z jens $
-
-#USE IN GALAXY
-#This script extract summary values (mean, min, max, std or coverage) from a bigwig file for a number of equal sized bins across genomic regions given in bed file.
-#If the bed file contains 6 columns (or more), column 6 is expected to contain strand information. Summary values from a negative strand will be reversed.
-
-########################################
-#bigwigSummary has 3 non-standard outputs:
-#1) n/a #(no data in bin)
-#2) no data #(no data in entire region)
-#3) is not a valid option (typically if negative coordinate)
-
-#Default settings & input parsing. "" indicates required user input.
-nbins=1
-rm_header_line=0
-summary_type=mean
-
-#parse input
-while getopts hef:b:o:n:t: myarg
-do case "$myarg" in
- h) echo "Usage: bigwig2summary_stdout -f -b -o -n -t -e"
- echo "Corrects for strand if bed-file has 6 columns or more. Col 6 assumed to contain strand infomation"
- exit ;;
- f) bigwig_file="$OPTARG" ;;
- b) bed_file="$OPTARG" ;; #must be tab separated without header
- o) org_assembly="$OPTARG" ;;
- n) nbins="$OPTARG" ;;
- t) summary_type="$OPTARG" ;;
- e) rm_header_line=2 ;; #flag. if -e then first line is removed
- [?]) echo "Usage: bigwig2summary_stdout -f -b -o -n -t -e"
- echo "Corrects for strand if bed-file has 6 columns or more. Col 6 assumed to contain strand infomation"
- exit 1 ;;
- esac
-done
-
-###################################################
-###VALIDATE INPUT
-###################################################
-
-#get chromosome sizes from bigwig file. bigwig file does not contain name of genome assembly.
-org_assembly_file=`mktemp -u`
-fetchChromSizes $org_assembly 2>/dev/null > $org_assembly_file
-if [ $? -ne 0 ]; then
- echo "ERROR: Organism genome assembly does not exist"
- rm $org_assembly_file
- exit
-fi
-
-#check input bed_file. bedClip only checks first 3 columns!
-if [ `bedClip -verbose=2 <(tail -n +${rm_header_line} $bed_file) $org_assembly_file /dev/null 2>&1 | wc -l` -gt 0 ]; then
- echo -e "ERROR: Input bed file is not in proper format!\nTry 'bedClip' to find lines causing error"
- echo "Make sure that bigwig and bed files are using the same genome assembly"
- exit 1
-fi
-
-#make string of "nbins" 0's to insert in regions, where no reads are found
-if [ $nbins -gt 1 ]; then
- seq_string=`seq 1 $nbins`
- zero_string=`printf "0\t%.s" {$seq_string} | perl -pe "s/\t$//"`
-fi
-
-#make sure the given summary type exists
-if [ `echo $summary_type | egrep "(mean|max|min|std|coverage)" | wc -l` -ne 1 ]; then
- echo "ERROR: Summary type must be: mean, max, min, std or coverage. Default is 'mean'"
- exit 1
-fi
-
-#determine number of fields in bed_file
-if [ `tail -n +${rm_header_line} $bed_file | awk '{print NF}' | uniq | wc -l` -eq 1 ]; then
- nfields_bed=`tail -n +${rm_header_line} $bed_file | awk '{print NF}' | uniq`
-else
- echo "ERROR: Bed file does not have constant number of line columns"
- exit 1
-fi
-
-if [[ $nbins -gt 1 && $nfields_bed -ge 6 ]]; then
- strand_uniq_chars=`tail -n +${rm_header_line} $bed_file | cut -f6 | sort -u | perl -pe "s/\n//"`
- if [[ $strand_uniq_chars != "+-" && $strand_uniq_chars != "-+" ]] ; then
- echo "ERROR: Column 6 in bed file must only contain '+' or '-' characters"
- exit 1
- fi
-fi
-
-###################################################
-###EXTRACT DENSITIES FROM NORMALIZED BIGWIG FILE
-###################################################
-
-
-#if more than 1 bin AND >= 6 fields (i.e. has strand column)
-if [[ $nbins -gt 1 && $nfields_bed -ge 6 ]]; then
-
- #cut columns 1-3+6 | rm header if flag set
- cut -f1-3,6 $bed_file | tail -n +${rm_header_line} | while read -r line; do
-
- #read 4 fields into variables
- read -r cur_chr cur_start cur_end cur_strand <<<"$line"
- #run bigWigSummary. Combine stdout and stderr | treat exceptions and errors after 'done'
- bigWigSummary $bigwig_file $cur_chr $cur_start $cur_end $nbins -type=${summary_type} 2>&1 | perl -pe "s/no data.+$/${zero_string}/" | awk 'BEGIN{OFS="\t"}{ if("'"$cur_strand"'"=="-") { for (i=NF; i>0; i--) { printf("%s\t",$i) } printf("\n") } else { print $0 } }'
- done | perl -pe "s/n\/a/0/g" | perl -pe "s/^\t?0\t?$/${zero_string}/" | perl -pe "s/ +/\t/g" | perl -pe "s/\t$//" | sed '/^$/d'
-
-#if more than 1 bin AND less than 6 fields (i.e. no strand column)
-elif [[ $nbins -gt 1 && $nfields_bed -lt 6 ]]; then
-
- #cut columns 1-3 | rm header if flag set
- cut -f1-3 $bed_file | tail -n +${rm_header_line} | while read -r line; do
-
- #read 3 fields into variables
- read -r cur_chr cur_start cur_end <<<"$line"
- #run bigWigSummary. Combine stdout and stderr | treat exceptions and errors after 'done'
- bigWigSummary $bigwig_file $cur_chr $cur_start $cur_end $nbins -type=${summary_type} 2>&1
-
- done | perl -pe "s/n\/a/0/g" | perl -pe "s/no data.+$/${zero_string}/" | perl -pe "s/^\t?0\t?$/${zero_string}/" | perl -pe "s/ +/\t/g" | sed '/^$/d'
-
-
-#if 1 bin. Strand column irrelevant
-else
-
- cut -f1-3 $bed_file | tail -n +${rm_header_line} | while read -r line; do
-
- read -r cur_chr cur_start cur_end <<<"$line"
- bigWigSummary $bigwig_file $cur_chr $cur_start $cur_end 1 -type=${summary_type} 2>&1
-
- done | perl -pe "s/no data.+$/0/" | perl -pe "s/n\/a/0/g" | sed '/^$/d'
-
-fi
-
-
-
-
-
diff -r 0847315fa035 -r a47f03c50780 bigwig2summary.xml
--- a/bigwig2summary.xml Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,106 +0,0 @@
-
- Extract summary information from a bigWig file across multiple genomic regions specified by the user.
-
- bigWigSummary
- fetchChromSizes
- bedClip
-
-
- bigwig2summary.sh -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_, bedClip_ and fetchChromSizes_.
-
-.. _bigWigSummary: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bigWigSummary
-
-.. _bedClip: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bedClip
-
--- _fetchChromSizes: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/fetchChromSizes
-
------
-
-.. 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 0847315fa035 -r a47f03c50780 faCount.xml
--- a/faCount.xml Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,78 +0,0 @@
-
- Count base statistics and CpGs in FASTA files.
-
- faCount
-
-
- faCount $summary $dinuc $strands $input > $output
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-**What it does**
-
-faCount_ - Count base statistics and CpGs in FA files.
-
-.. _faCount: http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/faCount
-
-**Usage**
-
-faCount file(s).fa -summary -dinuc -strands
-
-**Examples**
-
-Example 1::
-
- faCount 1.fasta
-
- #seq len A C G T N cpg
- HSFAU1 515 125 138 146 106 0 23
- HSFAU2 514 124 138 146 106 0 25
- HSFAU3 518 125 139 149 105 0 25
- HSFAU4 524 128 142 148 106 0 26
- HSFAU5 518 124 138 147 109 0 25
- total 2589 626 695 736 532 0 124
-
-Example 2::
-
- faCount 1.fasta -summary
-
- #seq len A C G T N cpg
- total 2589 626 695 736 532 0 124
- prcnt 1.0 0.2418 0.2684 0.2843 0.2055 0.0000 0.0479
-
-Example 3::
-
- faCount 1.fasta -summary -strands
-
- #seq len A C G T N cpg
- total 5178 1158 1431 1431 1158 0 248
- prcnt 1.0 0.2236 0.2764 0.2764 0.2236 0.0000 0.0479
-
-
-
-
diff -r 0847315fa035 -r a47f03c50780 test-data/.DS_Store
Binary file test-data/.DS_Store has changed
diff -r 0847315fa035 -r a47f03c50780 test-data/._.DS_Store
Binary file test-data/._.DS_Store has changed
diff -r 0847315fa035 -r a47f03c50780 test-data/1.bed
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1.bed Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,5 @@
+chr19 50178708 50180708
+chr6 90348174 90350174
+chr16 58495848 58497848
+chr5 180580242 180582242
+chr9 120177017 120179017
diff -r 0847315fa035 -r a47f03c50780 test-data/1.bigwig
Binary file test-data/1.bigwig has changed
diff -r 0847315fa035 -r a47f03c50780 test-data/1.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1.fasta Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,50 @@
+>HSFAU1
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggccccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU2
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgcactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU3
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaagggggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU4
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgaaatagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttgcccgaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU5
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtaggccgcatgctttttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
\ No newline at end of file
diff -r 0847315fa035 -r a47f03c50780 test-data/1.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1.tabular Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,4 @@
+0.144886 0 0
+0.507327 1.14649 1.38456
+0.221471 0.144886 0.309857
+0.348944 0.426638 0.244495
diff -r 0847315fa035 -r a47f03c50780 test-data/1.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1.txt Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,1 @@
+0.796113 0.273952 0.437672
diff -r 0847315fa035 -r a47f03c50780 test-data/2.bed
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2.bed Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,5 @@
+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
diff -r 0847315fa035 -r a47f03c50780 test-data/2.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2.tabular Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,4 @@
+0 0 0 0.144886 0
+0.396878 0.970009 1.1091 1.17538 1.40612
+0.203829 0.21262 0.144886 0.409996 0.274558
+0.180424 0.430989 0.458898 0.35171 0.307157
diff -r 0847315fa035 -r a47f03c50780 test-data/3.bed
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/3.bed Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,7 @@
+chr19 50178708 50180708
+chr21 48119895 48139895
+chr6 90348174 90350174
+chr16 58495848 58497848
+chr5 180580242 180582242
+chr17 81194210 81196210
+chr9 120177017 120179017
diff -r 0847315fa035 -r a47f03c50780 test-data/3.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/3.tabular Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,7 @@
+#seq len A C G T N cpg
+HSFAU1 515 125 138 146 106 0 23
+HSFAU2 514 124 138 146 106 0 25
+HSFAU3 518 125 139 149 105 0 25
+HSFAU4 524 128 142 148 106 0 26
+HSFAU5 518 124 138 147 109 0 25
+total 2589 626 695 736 532 0 124
diff -r 0847315fa035 -r a47f03c50780 test-data/4.bed
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/4.bed Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,5 @@
+chr19 50178708 50183708 NM_198318 0 + PRMT1
+chr6 90345174 90350174 NM_020466 0 - LYRM2
+chr16 58495848 58500848 NM_020465 0 + NDRG4
+chr5 180580242 180585242 NM_206880 0 + OR2V2
+chr9 120174017 120179017 NM_014010 0 - ASTN2
diff -r 0847315fa035 -r a47f03c50780 test-data/4.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/4.tabular Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,3 @@
+#seq len A C G T N cpg AA AC AG AT CA CC CG CT GA GC GG GT TA TC TG TT
+total 5178 1158 1431 1431 1158 0 248 360 238 408 142 328 447 248 408 320 426 447 238 150 320 328 360
+prcnt 1.0 0.2236 0.2764 0.2764 0.2236 0.0000 0.0479 0.0695 0.0460 0.0788 0.0274 0.0633 0.0863 0.0479 0.0788 0.0618 0.0823 0.0863 0.0460 0.0290 0.0618 0.0633 0.0695
diff -r 0847315fa035 -r a47f03c50780 test-data/test-data/1.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test-data/1.fasta Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,50 @@
+>HSFAU1
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggccccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU2
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgcactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU3
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaagggggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU4
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgaaatagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtagcaggccgcatgcttgcccgaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
+>HSFAU5
+ttcctctttctcgactccatcttcgcggtagctgggaccgccgttcagtcgccaatatgc
+agctctttgtccgcgcccaggagctacacaccttcgaggtgaccggccaggaaacggtcg
+cccagatcaaggctcatgtagcctcactggagggcattgccccggaagatcaagtcgtgc
+tcctggcaggcgcgcccctggaggatgaggccactctgggccagtgcggggtggaggccc
+tgactaccctggaagtaggccgcatgctttttggaggtaaagttcatggttccctggccc
+gtgctggaaaagtgagaggtcagactcctaaggtggccaaacaggagaagaagaagaaga
+agacaggtcgggctaagcggcggatgcagtacaaccggcgctttgtcaacgttgtgccca
+cctttggcaagaagaagggccccaatgccaactcttaagtcttttgtaattctggctttc
+tctaataaaaaagccacttagttcagtcaaaaaaaaaa
\ No newline at end of file
diff -r 0847315fa035 -r a47f03c50780 test-data/test-data/3.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test-data/3.tabular Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,7 @@
+#seq len A C G T N cpg
+HSFAU1 515 125 138 146 106 0 23
+HSFAU2 514 124 138 146 106 0 25
+HSFAU3 518 125 139 149 105 0 25
+HSFAU4 524 128 142 148 106 0 26
+HSFAU5 518 124 138 147 109 0 25
+total 2589 626 695 736 532 0 124
diff -r 0847315fa035 -r a47f03c50780 test-data/test-data/4.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test-data/4.tabular Tue Sep 09 10:28:01 2014 -0400
@@ -0,0 +1,3 @@
+#seq len A C G T N cpg AA AC AG AT CA CC CG CT GA GC GG GT TA TC TG TT
+total 5178 1158 1431 1431 1158 0 248 360 238 408 142 328 447 248 408 320 426 447 238 150 320 328 360
+prcnt 1.0 0.2236 0.2764 0.2764 0.2236 0.0000 0.0479 0.0695 0.0460 0.0788 0.0274 0.0633 0.0863 0.0479 0.0788 0.0618 0.0823 0.0863 0.0460 0.0290 0.0618 0.0633 0.0695
diff -r 0847315fa035 -r a47f03c50780 tool_dependencies.xml
--- a/tool_dependencies.xml Tue Sep 09 10:22:40 2014 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,88 +0,0 @@
-
-
-
-
-
- http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/bigWigSummary
-
- .
- $INSTALL_DIR
-
- $INSTALL_DIR/bigWigSummary
-
- $INSTALL_DIR
-
-
-
-
-
-
-
-
-
- http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/bedClip
-
- .
- $INSTALL_DIR
-
- $INSTALL_DIR/bedClip
-
- $INSTALL_DIR
-
-
-
-
-
-
-
-
-
- http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/bedExtendRanges
-
- .
- $INSTALL_DIR
-
- $INSTALL_DIR/bedExtendRanges
-
- $INSTALL_DIR
-
-
-
-
-
-
-
-
-
- http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/faCount
-
- .
- $INSTALL_DIR
-
- $INSTALL_DIR/faCount
-
- $INSTALL_DIR
-
-
-
-
-
-
-
-
-
- http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/fetchChromSizes
-
- .
- $INSTALL_DIR
-
- $INSTALL_DIR/fetchChromSizes
-
- $INSTALL_DIR
-
-
-
-
-
-
-