diff --git a/test/getfasta/test-getfasta.sh b/test/getfasta/test-getfasta.sh index 64dbcd55995dbeee46b07be902bfbdbbb040f012..fcc06c45d5b7c4055468db4b436c0f25afd477d4 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 f95782a9a3853a0cf2fd0e68a530a58f70da5fdc..424bd4df75fe030c569e7ae408ca01d1cc8498b2 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"