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
|
#!/bin/sh -e
pkg=last-align
if [ "$ADTTMP" = "" ] ; then
ADTTMP=$(mktemp -d /tmp/${pkg}-test.XXXXXX)
trap "rm -rf $ADTTMP" 0 INT QUIT ABRT PIPE TERM
fi
cd $ADTTMP
cp -a /usr/share/doc/${pkg}/examples/*.fa* .
gunzip -r *
# first example from tutorial
lastdb -cR01 humdb humanMito.fa
lastal humdb fuguMito.fa > myalns.maf
# multiMito example
lastdb -c humanMito humanMito.fa
lastal -e25 -j4 humanMito mouseMito.fa | last-split | maf-sort > hm.maf
lastal -e25 -j4 humanMito chickenMito.fa | last-split | maf-sort > hc.maf
lastal -e25 -j4 humanMito fuguMito.fa | last-split | maf-sort > hf.maf
maf-join hm.maf hc.maf hf.maf > out
|
