|
5
|
1
|
|
0
|
2 #!/bin/bash
|
|
|
3 dir="$(cd "$(dirname "$0")" && pwd)"
|
|
|
4
|
|
|
5 testID=$1
|
|
|
6 species=$2
|
|
|
7 substitutionModel=$3
|
|
|
8 mutabilityModel=$4
|
|
|
9 clonal=$5
|
|
|
10 fixIndels=$6
|
|
|
11 region=$7
|
|
2
|
12 inputs=$8
|
|
|
13 inputs=($inputs)
|
|
|
14 IDs=$9
|
|
|
15 IDs=($IDs)
|
|
|
16 ref=${10}
|
|
|
17 output=${11}
|
|
3
|
18 selection=${12}
|
|
14
|
19 output_table=${13}
|
|
2
|
20 outID="result"
|
|
|
21
|
|
14
|
22 echo "$PWD"
|
|
|
23
|
|
2
|
24 echo "testID = $testID"
|
|
|
25 echo "species = $species"
|
|
|
26 echo "substitutionModel = $substitutionModel"
|
|
|
27 echo "mutabilityModel = $mutabilityModel"
|
|
|
28 echo "clonal = $clonal"
|
|
|
29 echo "fixIndels = $fixIndels"
|
|
|
30 echo "region = $region"
|
|
|
31 echo "inputs = ${inputs[@]}"
|
|
|
32 echo "IDs = ${IDs[@]}"
|
|
|
33 echo "ref = $ref"
|
|
|
34 echo "output = $output"
|
|
|
35 echo "outID = $outID"
|
|
|
36
|
|
|
37 fasta="$PWD/baseline.fasta"
|
|
|
38
|
|
|
39
|
|
|
40 count=0
|
|
|
41 for current in ${inputs[@]}
|
|
|
42 do
|
|
|
43 f=$(file $current)
|
|
|
44 zipType="Zip archive"
|
|
11
|
45 if [[ "$f" == *"$zipType"* ]] || [[ "$f" == *"XZ compressed data"* ]]
|
|
2
|
46 then
|
|
|
47 id=${IDs[$count]}
|
|
|
48 echo "id=$id"
|
|
11
|
49 if [[ "$f" == *"Zip archive"* ]] ; then
|
|
|
50 echo "Zip archive"
|
|
|
51 echo "unzip $input -d $PWD/files/"
|
|
|
52 unzip $current -d "$PWD/$id/"
|
|
|
53 elif [[ "$f" == *"XZ compressed data"* ]] ; then
|
|
|
54 echo "ZX archive"
|
|
|
55 echo "tar -xJf $input -C $PWD/files/"
|
|
|
56 mkdir -p "$PWD/$id/files"
|
|
|
57 tar -xJf $current -C "$PWD/$id/files/"
|
|
|
58 fi
|
|
3
|
59 summaryfile="$PWD/summary_${id}.txt"
|
|
|
60 gappedfile="$PWD/gappednt_${id}.txt"
|
|
|
61 filtered="$PWD/filtered_${id}.txt"
|
|
6
|
62 filecount=`ls -l $PWD/$id/ | wc -l`
|
|
7
|
63 if [[ "$filecount" -eq "2" ]]
|
|
6
|
64 then
|
|
|
65 cat $PWD/$id/*/1_* > $summaryfile
|
|
|
66 cat $PWD/$id/*/2_* > $gappedfile
|
|
|
67 else
|
|
|
68 cat $PWD/$id/1_* > $summaryfile
|
|
|
69 cat $PWD/$id/2_* > $gappedfile
|
|
|
70 fi
|
|
5
|
71 Rscript $dir/filter.r $summaryfile $gappedfile "$selection" $filtered 2>&1
|
|
3
|
72
|
|
|
73 final="$PWD/final_${id}.txt"
|
|
|
74 cat $filtered | cut -f2,4,7 > $final
|
|
|
75 python $dir/script_imgt.py --input $final --ref $ref --output $fasta --id $id
|
|
2
|
76 else
|
|
|
77 python $dir/script_xlsx.py --input $current --ref $ref --output $fasta
|
|
|
78 fi
|
|
|
79 count=$((count+1))
|
|
|
80 done
|
|
|
81 workdir="$PWD"
|
|
|
82 cd $dir
|
|
8
|
83 echo "file: ${inputs[0]}"
|
|
|
84 #Rscript --verbose $dir/Baseline_Main.r $testID $species $substitutionModel $mutabilityModel $clonal $fixIndels $region ${inputs[0]} $workdir/ $outID 2>&1
|
|
2
|
85 Rscript --verbose $dir/Baseline_Main.r $testID $species $substitutionModel $mutabilityModel $clonal $fixIndels $region $fasta $workdir/ $outID 2>&1
|
|
|
86
|
|
|
87 echo "$workdir/${outID}.txt"
|
|
|
88
|
|
5
|
89 rows=`tail -n +2 $workdir/${outID}.txt | grep -v "All sequences combined" | grep -n 'Group' | grep -Eoh '^[0-9]+' | tr '\n' ' '`
|
|
2
|
90 rows=($rows)
|
|
5
|
91 #unset rows[${#rows[@]}-1]
|
|
0
|
92
|
|
|
93 cd $dir
|
|
2
|
94 Rscript --verbose $dir/comparePDFs.r $workdir/${outID}.RData $output ${rows[@]} 2>&1
|
|
14
|
95 cp $workdir/result.txt ${output_table}
|
|
0
|
96
|
|
|
97
|
|
2
|
98
|
|
0
|
99
|
