File: extract-seqs.sh

package info (click to toggle)
genometools 1.6.1%2Bds-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 50,412 kB
  • sloc: ansic: 271,241; ruby: 30,339; python: 4,880; sh: 3,193; makefile: 1,194; perl: 219; pascal: 159; haskell: 37; sed: 5
file content (14 lines) | stat: -rwxr-xr-x 219 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
#!/bin/sh

if test $# -lt 2
then
  echo "Usage: $0 <fastafile> seqnum1 [seqnum2 ...]"
  exit 1
fi

inputfile=$1
shift
for seqnum in $*
do
  env -i bin/gt seq -width 70 -showseqnum `expr ${seqnum} + 1` ${inputfile}
done