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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
#!/bin/sh
# Author: Lance Lin <LQi254@protonmail.com>
# Date: 18 February 2022
# File: run-unit-test
# Purpose: run a simple-unit test for the lagan package
#
# Details:
# Use a publicly accessible .fasta file format example. There are others, but they quite large (>30MB).
# Compare it with the included sample.fasta file in the lagan package.
set -e
if [ "${AUTOPKGTEST_TMP}" = "" ] ; then
AUTOPKGTEST_TMP=$(mktemp -d /tmp/${pkg}-test.XXXXXX)
trap "rm -rf ${AUTOPKGTEST_TMP}" 0 INT QUIT ABRT PIPE TERM
fi
cd "${AUTOPKGTEST_TMP}"
# This is needed for lagan to run
export LAGAN_DIR=/usr/lib/lagan
echo "run-unit-test: generating sample data"
# FASTA format file from "UT Southwestern Medical Center"
# URL: http://prodata.swmed.edu/promals/info/fasta_format_file_example.htm
# Website title: Fasta format file example
# Date accessed: 18 February 2022
cat <<EOF > testfile.fasta
>seq0
FQTWEEFSRAAEKLYLADPMKVRVVLKYRHVDGNLCIKVTDDLVCLVYRTDQAQDVKKIEKF
>seq1
KYRTWEEFTRAAEKLYQADPMKVRVVLKYRHCDGNLCIKVTDDVVCLLYRTDQAQDVKKIEKFHSQLMRLME LKVTDNKECLKFKTDQAQEAKKMEKLNNIFFTLM
>seq2
EEYQTWEEFARAAEKLYLTDPMKVRVVLKYRHCDGNLCMKVTDDAVCLQYKTDQAQDVKKVEKLHGK
>seq3
MYQVWEEFSRAVEKLYLTDPMKVRVVLKYRHCDGNLCIKVTDNSVCLQYKTDQAQDVK
>seq4
EEFSRAVEKLYLTDPMKVRVVLKYRHCDGNLCIKVTDNSVVSYEMRLFGVQKDNFALEHSLL
>seq5
SWEEFAKAAEVLYLEDPMKCRMCTKYRHVDHKLVVKLTDNHTVLKYVTDMAQDVKKIEKLTTLLMR
>seq6
FTNWEEFAKAAERLHSANPEKCRFVTKYNHTKGELVLKLTDDVVCLQYSTNQLQDVKKLEKLSSTLLRSI
>seq7
SWEEFVERSVQLFRGDPNATRYVMKYRHCEGKLVLKVTDDRECLKFKTDQAQDAKKMEKLNNIFF
>seq8
SWDEFVDRSVQLFRADPESTRYVMKYRHCDGKLVLKVTDNKECLKFKTDQAQEAKKMEKLNNIFFTLM
>seq9
KNWEDFEIAAENMYMANPQNCRYTMKYVHSKGHILLKMSDNVKCVQYRAENMPDLKK
>seq10
FDSWDEFVSKSVELFRNHPDTTRYVVKYRHCEGKLVLKVTDNHECLKFKTDQAQDAKKMEK
EOF
# This is sample.fasta from the lagan package
cat <<EOF > sample.fasta
>sample1
GGCATGTCCAGAAAATCCAAGTGCCTCTTCCTCTTGATCTTCTCCAACGATGTCCAGA
AAATCCAAGTGCCTCATTCCTCTTGATCTTCTCCAGGCATGTCCAGAAAATCCAAGTG
CCTCTTCCTCTCTGATCTTCTCCTCGGTTGGTCCAGAAAATCCAAGTGCCTCTTCCTC
TTGATCTTCTCCAGAAATGTCCAGAAAATCCAAGTAGCCTCTTCCTCTTGATCGGCTC
CAGAAATGTCCAGAAAAATCCAAGTGCCTCTTCCTCTTGATCGGCTCCATAAATGTCC
AGAAAATCCAACGTGCCTCTTCCTCTTGATCGGCTCCAGAAATGTCCAGAAATATCCA
AGTGCCTCTTCCTCTTGATCGGCTCCTTA
>sample2
CGATCCCAAATCCAAGTGCCTCAGAGTCTACTTGATCTTCAATTCAGATCCCAAATCC
AAGTGCCTCAGAGTCTACTTGAATCTTCTATCGGGTCCCAAATCCAAGTGCCTCAGAG
TCTACTTGATCTTCTCTCTCGATCCCATATCCAAGTGCCTCCTAGAGTCTACTTGATC
TTCTCGATAACCAAAATCCAAGTGCCTCAGAGTCTACTTCACTCTTCTCGACTAACCC
AAATCCAAGTGCCTCAGATGAGTCTACTTCCTCTTCTCATAACTCAAATCCAAGTGCC
TCAGAGTCTAACTTCCTCTTCTCGAATAACCCAAATCCAAGTGCCTCAGAGTGTCTAC
TTCCTCTTCTCG
>sample3
TACCCAAATCCAAGTGCCTCAGCGTCTAATAAAACAAGTCTTGATCTTCAACTCCTCC
CAAATCCAAGTGCAACCTCAGCCGCTAATAAAAAGTCTTGATCTTCTCGCGTCCGGCA
AATCCAAGTGCCTCAGCGCTAATAAAAAGTCTTGATCTTCTCGGGAGTCCCAAATCCA
AGTGCCTCAGCGCTAATAAAAAGTCTTGATCTTCTCGGAGGAACAACAAATCCAAGTG
CCTCAAGCGCTAATAAAAAGTCCCGATCTTCTCGTGACAATACAAATCCAAGTGCCTC
AGCGCTAATAAAAAGTCCCGATCTTCTCCCGTGTAAACAAATCCAAGTGCCTCAGCGC
TAATAAAAAGTCCCGATCTTCTCTGGTAACACAACAAATCCAAGTCACGCCTCAGATA
CGCTAATAAAAAGTCCCGATCTTCTC
EOF
echo "run-unit-test: generating md5 sums for verification"
echo "bb24ed1e1fb416d76a27580cd874c1f7 output.txt" >> checksums
echo "run-unit-test: running lagan/order"
order testfile.fasta sample.fasta > output.txt
echo "run-unit-test: verifying md5 checksums"
md5sum --check checksums
|
