|
2
|
1 #!/bin/bash
|
|
|
2
|
|
|
3 tool_path=$(dirname $0)
|
|
|
4
|
|
|
5 filein=$1
|
|
3
|
6 fileout_label=$(date "+%Y%m%d%H%M%S")
|
|
|
7 fileout=$2
|
|
|
8 filelog=$3
|
|
|
9 export=$4
|
|
|
10 frequency=$5
|
|
|
11 max_freq=$6
|
|
|
12 allow_missing=$7
|
|
|
13 nb_alleles_min=$8
|
|
|
14 nb_alleles_max=9
|
|
|
15 type=${10}
|
|
|
16 bound_start=${11}
|
|
|
17 bound_end=${12}
|
|
2
|
18
|
|
|
19
|
|
4
|
20 if [ "${13}" != "None" ]
|
|
|
21 then samples="--samples ${13}"
|
|
2
|
22 fi
|
|
|
23
|
|
4
|
24 if [ "${14}" != "None" ]
|
|
|
25 then chromosomes="--chromosomes ${14}"
|
|
2
|
26 fi
|
|
|
27
|
|
|
28 if [ "$bound_start" -gt "$bound_end" ]
|
|
|
29 then tmp=$bound_start ; bound_start=$bound_end ; bound_end=$tmp ; echo "Warning : Lower bound must be lower than greater bound!" >&2
|
|
|
30 fi
|
|
|
31
|
|
|
32 if [ "$nb_alleles_min" -gt "$nb_alleles_max" ]
|
|
|
33 then tmp=$nb_alleles_min ; nb_alleles_min=$nb_alleles_max ; nb_alleles_max=$tmp ; echo "Warning : Minimum number of alleles must be lower than maximum number of allele!" >&2
|
|
|
34 fi
|
|
|
35
|
|
|
36 perl $tool_path/VCFToolsFilter.pl --input $filein --out $fileout_label --export $export --frequency $frequency --max_freq $max_freq --allow_missing $allow_missing --nb_alleles $nb_alleles_min','$nb_alleles_max --type $type --bounds $bound_start','$bound_end $samples $chromosomes
|
|
|
37
|
|
|
38 if [ "$export" = "VCF" ]
|
|
|
39 then cp $fileout_label.recode.vcf $fileout ; rm $fileout_label.recode.vcf
|
|
|
40 elif [ "$export" = "freq" ]
|
|
|
41 then cp $fileout_label.frq $fileout ; rm $fileout_label.frq
|
|
4
|
42 else cp $fileout_label.ped $fileout; cp $fileout_label.map ${15} ; rm $fileout_label.ped $fileout_label.map
|
|
2
|
43 fi
|
|
|
44
|
|
|
45 cp vcftools.log $filelog
|
|
|
46 rm vcftools.log
|
