1 |
gpertea |
23 |
#!/bin/tcsh |
2 |
|
|
# if ($1'x' == 'x') then |
3 |
|
|
# echo "Usage: ghmm_remap.csh <ghmm..." |
4 |
|
|
# echo |
5 |
|
|
# exit 1 |
6 |
|
|
# endif |
7 |
|
|
|
8 |
|
|
|
9 |
|
|
foreach f ( *.isoFa.[1-9] *.isoFa.[1-9][0-9] ) |
10 |
|
|
set ghpd = $f:s/.isoFa./.isoFa.ghmm.pdom./ |
11 |
|
|
set gh = $f:s/.isoFa./.isoFa.ghmm./ |
12 |
|
|
set ghopd = $f:r |
13 |
|
|
set ghopd = $ghopd:s/.isoFa/.glimmerhmm_pdom.gff/ |
14 |
|
|
set gho = $f:r |
15 |
|
|
set gho = $gho:s/.isoFa/.glimmerhmm.gff/ |
16 |
|
|
set ofs = `head -1 $f | cut -f2 -d ' '` |
17 |
|
|
mapgffback.pl $ghpd $ofs >> $ghopd |
18 |
|
|
mapgffback.pl $gh $ofs >> $gho |
19 |
|
|
#exit |
20 |
|
|
end |