0
|
1 #!/bin/bash
|
|
2
|
|
3 if [ $# != 4 ]; then
|
|
4 echo " ==== ERROR ... you called this script inappropriately."
|
|
5 echo ""
|
9
|
6 echo " usage: $0 fragsFILE.frags/.csv fastaX fastaY alignments.txt"
|
0
|
7 echo ""
|
|
8 exit -1
|
|
9 fi
|
|
10
|
|
11 FRAGS=$1
|
|
12 FASTAX=$2
|
|
13 FASTAY=$3
|
|
14 ALIGN=$4
|
|
15
|
|
16
|
|
17
|
|
18 BINDIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"
|
|
19
|
|
20 # generate indices
|
|
21
|
9
|
22 $BINDIR/indexmaker $FASTAX $FASTAX.idx
|
|
23 $BINDIR/indexmaker $FASTAY $FASTAY.idx
|
0
|
24
|
9
|
25 $BINDIR/reverseComplement $FASTAY $FASTAY.rev
|
|
26 $BINDIR/indexmaker $FASTAY.rev $FASTAY.rev.idx
|
0
|
27
|
|
28 # extract frags
|
|
29
|
9
|
30 if [ ${FRAGS: -4} == ".csv" ]
|
|
31 then
|
|
32
|
|
33 echo "....using csv"
|
|
34 sed 's/,/ /g' $FRAGS > $FRAGS.fix
|
|
35 $BINDIR/csvFrags2text $FRAGS.fix $FASTAX $FASTAY $FASTAY.rev $FASTAX.idx $FASTAY.idx $FASTAY.rev.idx $ALIGN
|
|
36 rm $FRAGS.fix
|
0
|
37
|
9
|
38 else
|
|
39 echo "....using frags"
|
|
40 $BINDIR/frags2text $FRAGS $FASTAX $FASTAY $FASTAY.rev $FASTAX.idx $FASTAY.idx $FASTAY.rev.idx $ALIGN
|
|
41 fi
|
|
42
|
0
|
43
|
9
|
44 rm $FASTAX.idx
|
|
45 rm $FASTAY.idx
|
|
46 rm $FASTAY.rev
|
|
47 rm $FASTAY.rev.idx
|
|
48
|