view VCF2Hapmap/vcf2FastaAndHapmap.sh @ 14:d15869b3731a draft

Deleted selected files
author dereeper
date Tue, 08 Jan 2019 08:46:57 -0500
parents 98c37a5d67f4
children
line wrap: on
line source

#!/bin/bash

tool_path=$(dirname $0)



filein=$1
fileout_label=$(date "+%Y%m%d%H%M%S")
fileout=$2
option=$3

option_text=''


if [ "$option" != "none" ]
then fileout_seq=$4 
	fileout_fa1=$5 
	filefasta=$6
	if [ "$option" == "fasta_gff" ]
		then filegff=$7
	fi
fi

if [ "$option" == "fasta" ]
then option_text="--reference $filefasta"
fi

if [ "$option" == "fasta_gff" ]
then option_text="--reference $filefasta --gff $filegff"
fi


perl $tool_path/VCF2FastaAndHapmap.pl --vcf $filein --out $fileout_label $option_text


cp  $fileout_label.hapmap $fileout ; rm $fileout_label.hapmap

if [ "$option" == "fasta_gff" ]
then cp  $fileout_label.flanking.txt $fileout_seq ; rm $fileout_label.flanking.txt ; cp  $fileout_label.gene_alignment.fas $fileout_fa1 ; rm $fileout_label.gene_alignment.fas ;
fi

if [ "$option" == "fasta" ]
then cp  $fileout_label.flanking.txt $fileout_seq ; rm $fileout_label.flanking.txt ;
fi