| 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 |
