Skip to content

Commit 67440aa

Browse files
committed
add stats
1 parent 6adc726 commit 67440aa

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

GetCoverageOfRegions.sh

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -2,10 +2,10 @@
22
bedFile=$1
33
covFile=$2
44
meanCov=`cat $3`
5-
5+
sdir=$4
66

77
while read line; do
88
rgn=`echo $line | awk '{ print $1":"$2"-"$3;}'`
9-
tabix $covFile $rgn | cut -f 4 | stats | awk -v m=$meanCov '{ print $1/m;}'
9+
tabix $covFile $rgn | cut -f 4 | $sdir/stats | awk -v m=$meanCov '{ print $1/m;}'
1010
done < $bedFile
1111

0 commit comments

Comments
 (0)