From d745b05e0f0c7bc712fc7e10e3dd44a76df88e0e Mon Sep 17 00:00:00 2001 From: Charles Plessy <plessy@debian.org> Date: Sun, 4 Nov 2012 19:09:10 +0900 Subject: [PATCH] Call bedtools via a shell variable called BT. --- test/getfasta/test-getfasta.sh | 14 ++++++++------ test/groupBy/test-groupby.sh | 14 ++++++++------ 2 files changed, 16 insertions(+), 12 deletions(-) diff --git a/test/getfasta/test-getfasta.sh b/test/getfasta/test-getfasta.sh index 64dbcd55..fcc06c45 100644 --- a/test/getfasta/test-getfasta.sh +++ b/test/getfasta/test-getfasta.sh @@ -6,6 +6,8 @@ # aggggggggg # cggggggggg +BT=../../bin/bedtools + check() { if $1 = $2; then @@ -16,7 +18,7 @@ check() } echo " getfasta.t01...\c" -LINES=$(echo $'chr1\t1\t10' | ../../bin/bedtools getfasta -fi t.fa -bed stdin -fo - | awk 'END{ print NR }') +LINES=$(echo $'chr1\t1\t10' | $BT getfasta -fi t.fa -bed stdin -fo - | awk 'END{ print NR }') if [ "$LINES" != "2" ]; then echo fail else @@ -24,7 +26,7 @@ else fi echo " getfasta.t02...\c" -LEN=$(../../bin/bedtools getfasta -split -fi t.fa -bed blocks.bed -fo stdout | awk '(NR == 2){ print length($0) }') +LEN=$($BT getfasta -split -fi t.fa -bed blocks.bed -fo stdout | awk '(NR == 2){ print length($0) }') if [ "$LINES" != "2" ]; then echo fail else @@ -32,7 +34,7 @@ else fi echo " getfasta.t03...\c" -SEQ=$(../../bin/bedtools getfasta -split -fi t.fa -bed blocks.bed -fo stdout | awk '(NR == 4){ print $0 }') +SEQ=$($BT getfasta -split -fi t.fa -bed blocks.bed -fo stdout | awk '(NR == 4){ print $0 }') if [ "$SEQ" != "cta" ]; then echo fail else @@ -41,7 +43,7 @@ fi # test -fo - echo " getfasta.t04...\c" -SEQ=$(../../bin/bedtools getfasta -split -fi t.fa -bed blocks.bed -fo - | awk '(NR == 4){ print $0 }') +SEQ=$($BT getfasta -split -fi t.fa -bed blocks.bed -fo - | awk '(NR == 4){ print $0 }') if [ "$SEQ" != "cta" ]; then echo fail else @@ -51,9 +53,9 @@ fi # test -split with -s - echo " getfasta.t05...\c" -SEQ=$(../../bin/bedtools getfasta -split -s -fi t.fa -bed blocks.bed -fo - | awk '(NR == 4){ print $0 }') +SEQ=$($BT getfasta -split -s -fi t.fa -bed blocks.bed -fo - | awk '(NR == 4){ print $0 }') if [ "$SEQ" != "tag" ]; then echo fail else echo ok -fi \ No newline at end of file +fi diff --git a/test/groupBy/test-groupby.sh b/test/groupBy/test-groupby.sh index f95782a9..424bd4df 100644 --- a/test/groupBy/test-groupby.sh +++ b/test/groupBy/test-groupby.sh @@ -1,7 +1,9 @@ -lines_a=$(../../bin/groupBy -g 3-1 -o collapse -c 4 -i ../map/values3.bed | wc -l) -lines_b=$(../../bin/groupBy -g 1-3 -o collapse -c 4 -i ../map/values3.bed | wc -l) -lines_c=$(../../bin/groupBy -g 1,2,3 -o collapse -c 4 -i ../map/values3.bed | wc -l) -lines_d=$(../../bin/groupBy -g 1-2,3 -o collapse -c 4 -i ../map/values3.bed | wc -l) +BT=../../bin/bedtools + +lines_a=$($BT groupby -g 3-1 -o collapse -c 4 -i ../map/values3.bed | wc -l) +lines_b=$($BT groupby -g 1-3 -o collapse -c 4 -i ../map/values3.bed | wc -l) +lines_c=$($BT groupby -g 1,2,3 -o collapse -c 4 -i ../map/values3.bed | wc -l) +lines_d=$($BT groupby -g 1-2,3 -o collapse -c 4 -i ../map/values3.bed | wc -l) check(){ if [ "$1" != "$2" ]; then @@ -15,13 +17,13 @@ check $lines_a $lines_d H=$(head -n 1 values3.header.bed) -A=$(../../bin/bedtools groupby -i values3.header.bed -g 1,2,3 -c 4 -o concat -inheader | head -n 1) +A=$($BT groupby -i values3.header.bed -g 1,2,3 -c 4 -o concat -inheader | head -n 1) if [ "$A" != $'chr1\t0\t10\ta1' ]; then echo "fail groupby" fi -B=$(../../bin/bedtools groupby -i values3.header.bed -g 1,2,3 -c 4 -o concat -header | head -n 1) +B=$($BT groupby -i values3.header.bed -g 1,2,3 -c 4 -o concat -header | head -n 1) if [ "$B" != $'#chrom\tstart\tend\tconcat(A)' ]; then echo "fail groupby" -- GitLab