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
|
#!/bin/sh -e
pkg=seqtk
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
cp -a /usr/share/doc/${pkg}/test_data/* .
gunzip -r *
seqtk seq -a in.fq > out.fa
seqtk seq -aQ64 -q20 in.fq > out.fa
seqtk seq -aQ64 -q20 -n N in.fq > out.fa
seqtk seq -Cl60 in.fa > out.fa
seqtk seq -l0 in.fq > out.fq
seqtk seq -r in.fq > out.fq
seqtk subseq in.fq name.lst > out.fq
seqtk sample -s100 in.fq 2 > sub.fq
seqtk trimfq in.fq > out.fq
seqtk trimfq -b 5 -e 10 in.fa > out.fa
|