1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
#!/bin/sh
set -e -x
if test $# -ne 3
then
echo "Usage: $0 <Referencefile> <Queryfile> <threshold>"
exit 1
fi
reference=$1
query=$2
threshold=$3
# exec in gt root dir
# mkvtree.x -db ${reference} -indexname esa-mkv -dna -pl -lcp -suf -tis
# vstree2tex.x -tis -suf -s esa-mkv > tmp.tex
# gt suffixerator -tis -suf -dna -pl -db ${reference} -indexname esa-fwd
# gt tagerator -e 2 -maxocc 20 -rw -esa esa-fwd -q ${query}
gt packedindex mkindex -ssp -tis -indexname pck-rev -db ${reference} -sprank -dna -pl -bsize 10 -locfreq 32 -dir rev
gt prebwt -maxdepth 4 -pck pck-rev
gt tagerator -rw -best -e ${threshold} -pck pck-rev -q ${query}
|
