|
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")
|
|
8
|
7 fileout_het=$2
|
|
|
8 fileout_imiss=$3
|
|
|
9 fileout_sum=$4
|
|
|
10 filelog=$5
|
|
2
|
11
|
|
|
12
|
|
|
13
|
|
|
14 perl $tool_path/VCFToolsStats.pl --input $filein --out $fileout_label
|
|
|
15
|
|
|
16 cp $fileout_label.het $fileout_het ; rm $fileout_label.het
|
|
|
17 cp $fileout_label.imiss $fileout_imiss ; rm $fileout_label.imiss
|
|
|
18 cp $fileout_label.TsTv.summary $fileout_sum ; rm $fileout_label.TsTv.summary
|
|
|
19
|
|
|
20 cp vcftools.log $filelog
|
|
|
21 rm vcftools.log
|
