Mercurial > repos > davidvanzessen > change_o

view create_germlines.sh @ 1:326165da9ece draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Thu, 04 Aug 2016 08:59:47 -0400
parents dda9b2e72e2b
children 6f57d43a1f4c
line wrap: on
line source

#!/bin/bash
dir="$(cd "$(dirname "$0")" && pwd)"

input=$1
type=$2
cloned=$3
output=$4

cp $input $PWD/input.tab #file has to have a ".tab" extension

if [ "true" == "$cloned" ] ; then
	cloned="--cloned"
else
	cloned=""
fi

mkdir $PWD/outdir

/home/galaxy/anaconda3/bin/python $dir/CreateGermlines.py -d $PWD/input.tab -r $dir/IMGT_Human_IGH[VDJ].fasta --outdir $PWD/outdir --outname output -g $type $cloned
#/data/users/david/anaconda3/bin/python $dir/CreateGermlines.py -d $PWD/input.tab -r $dir/IMGT_Human_IGH[VDJ].fasta --outdir $PWD/outdir --outname output -g $type $cloned

mv $PWD/outdir/output_germ-pass.tab $output