-
Notifications
You must be signed in to change notification settings - Fork 0
/
wrapper_de_utr.sh
executable file
·37 lines (27 loc) · 1.01 KB
/
wrapper_de_utr.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
#!/bin/bash
set -e
# Example ./bo_utr.sh inputs/Arabidopsis_thaliana.TAIR10.31.dna_rm.toplevel.fa inputs/Arabidopsis_thaliana.TAIR10.31.gff3 ID=gene: 500 --No --UseUTR
#inputs
# $1 - File, e.g. PlantA.fa
# $2 - File, e.g. PlantA.gff3
# $3 - String, e.g. ID=gene:
# $4 - Integer, e.g. 500/2000/5000
# $5 - --NoOverlap (if checked)
# $6 - --UseUTR (if checked)
fileA=$1
fileB=$2
FILENAMEA=${fileA##*/} # remove the path and leave only the file name
FILENAMEB=${fileB##*/}
mv ${fileA} /apples/bin/utr_tool/
mv ${fileB} /apples/bin/utr_tool/
cd /apples/bin/utr_tool
./bo_utr.sh $1 $2 $3 $4 $5 $6
mv /apples/bin/utr_tool/works/promoters.fa /de-app-work/PlantA.fa
mv /apples/bin/utr_tool/works/promoters.bed /de-app-work/PlantA.bed
mv /apples/bin/utr_tool/works/utr3.bed /de-app-work/PlantA_utr3.bed
mv /apples/bin/utr_tool/works/utr5.bed /de-app-work/PlantA_utr5.bed
cd /de-app-work
# Make easier for the user to keep track of what input was used
touch "FA_"${FILENAMEA}
touch "GFF_"${FILENAMEB}
echo "${@:1}" > Input_Parameters.txt