1
|
1 #!/bin/bash
|
|
2
|
|
3 FL=1000 # frequency limit
|
|
4
|
|
5 if [ $# != 7 ]; then
|
|
6 echo " ==== ERROR ... you called this script inappropriately."
|
|
7 echo ""
|
|
8 echo " usage: $0 seqXName seqYName lenght similarity WL fixedL strand"
|
|
9 echo ""
|
|
10 exit -1
|
|
11 fi
|
|
12
|
|
13 seqXName=$(basename "$1")
|
|
14 extensionX="${seqXName##*.}"
|
|
15 seqXName="${seqXName%.*}"
|
|
16
|
|
17 seqYName=$(basename "$2")
|
|
18 extensionY="${seqYName##*.}"
|
|
19 seqYName="${seqYName%.*}"
|
|
20
|
|
21 #seqXName=`basename $1 .fasta`
|
|
22 #seqYName=`basename $2 .fasta`
|
|
23
|
|
24 BINDIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"
|
|
25
|
|
26 length=${3}
|
|
27 similarity=${4}
|
|
28 WL=${5} # wordSize
|
|
29 fixedL=${6}
|
|
30 strand=${7}
|
|
31 distance=$((4*${WL}))
|
|
32
|
|
33 if [[ ! -f ../hits/${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered ]]; then
|
|
34 echo "${BINDIR}/hits ${seqXName} ${seqYName} ${seqXName}-${seqYName}-K${WL}.hits ${FL} ${WL}"
|
|
35 ${BINDIR}/hits ${seqXName} ${seqYName} ${seqXName}-${seqYName}-K${WL}.hits ${FL} ${WL}
|
|
36
|
|
37 echo "${BINDIR}/sortHits 10000000 32 ${seqXName}-${seqYName}-K${WL}.hits ${seqXName}-${seqYName}-K${WL}.hits.sorted"
|
|
38 ${BINDIR}/sortHits 10000000 32 ${seqXName}-${seqYName}-K${WL}.hits ${seqXName}-${seqYName}-K${WL}.hits.sorted
|
|
39
|
|
40 echo "${BINDIR}/filterHits ${seqXName}-${seqYName}-K${WL}.hits.sorted ${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered ${WL}"
|
|
41 ${BINDIR}/filterHits ${seqXName}-${seqYName}-K${WL}.hits.sorted ${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered ${WL}
|
|
42
|
|
43 mv ${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered ../hits/
|
|
44 fi
|
|
45
|
|
46 ln -s ../hits/${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered .
|
|
47
|
|
48 echo "${BINDIR}/FragHits $1 $2 ${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered ${seqXName}-${seqYName}-s${strand}.frags ${length} ${similarity} ${WL} ${fixedL} ${strand}"
|
|
49 ${BINDIR}/FragHits $1 $2 ${seqXName}-${seqYName}-K${WL}.hits.sorted.filtered ${seqXName}-${seqYName}-s${strand}.frags ${length} ${similarity} ${WL} ${fixedL} ${strand}
|
|
50
|
|
51 echo "--------------------DONE------------------"
|
|
52
|