view filter_vcf.sh @ 15:48c41c6ad84e draft

Uploaded

#!/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