File: spmop.sh

package info (click to toggle)
genometools 1.6.6%2Bds-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 50,576 kB
  • sloc: ansic: 271,876; ruby: 29,930; python: 5,106; sh: 3,083; makefile: 1,213; perl: 219; pascal: 159; haskell: 37; sed: 5
file content (21 lines) | stat: -rwxr-xr-x 315 bytes parent folder | download | duplicates (9) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
#!/bin/sh

set -e -x

if test $# -ne 1
then
  echo "Usage: $0 <filename>"
  exit 1
fi

runsfx()
{
  env -i bin/gt suffixerator -dna -suf -lcp -tis -spmopt 20 $*
}

cp $1 both.fna
reversecompl.pl 60 $1 >> both.fna

runsfx -indexname sfx1 -db $1 -mirrored
runsfx -indexname sfx2 -db both.fna
cmp -s sfx1.suf sfx2.suf