view filter_vcf.sh @ 10:76f76eee9130 draft

Uploaded
author ulfschaefer
date Thu, 16 Jun 2016 08:17:09 -0400
parents e2ee6acf8752
children
line wrap: on
line source

#!/bin/bash

# echo $@

OUTPUT=$1
shift
VCF=$1
shift
REF=$1
shift
ONLYGOOD=$1
shift
FILTERS=$@

FILTERS=$(echo -n $FILTERS | sed 's/ /,/g')

CMD="phenix.py filter_vcf --vcf $VCF --filters $FILTERS --output $OUTPUT --reference $REF"

if [ $ONLYGOOD != "NOTTHERE" ]; then
    CMD="$CMD --only-good"
fi

echo "CMD: "$CMD
eval $CMD