File: align.sh

package info (click to toggle)
squizz 0.99c%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 6,760 kB
  • ctags: 8,355
  • sloc: sh: 4,789; ansic: 2,640; lex: 1,999; yacc: 1,650; makefile: 119
file content (26 lines) | stat: -rwxr-xr-x 766 bytes parent folder | download | duplicates (3)
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
#! /bin/sh

## Inits
alifmt='CLUSTAL FASTA MSF NEXUS PHYLIPI PHYLIPS STOCKHOLM'

## Check sequences names with spaces
sed 's/\./ /' $srcdir/align/phylips.seq >_tmpali.seq
for f in $alifmt; do
  ../src/squizz -c $f _tmpali.seq >_tmpali2.seq 2>/dev/null || exit 1
  (../src/squizz _tmpali2.seq 2>&1 | grep -v $f) && exit 1
done

## Check for unexpected sequence name
sed -e '$d' $srcdir/align/clustal.seq | sed -e '$d' >_tmpali.seq
tail -n 2 $srcdir/align/clustal.seq | sed 's/_/-/' >>_tmpali.seq
(../src/squizz -A _tmpali.seq 2>&1 | grep -v UNKNOWN) && exit 1

## Check for duplicate sequence name
cat $srcdir/align/fasta.seq $srcdir/align/fasta.seq >_tmpali.seq
(../src/squizz -A _tmpali.seq 2>&1 | grep -v UNKNOWN) && exit 1

## Cleanup
rm -f _tmpali*.seq

exit 0