git://source.jalview.org / jalview.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (from parent 1: 75f9103)
in progress
authorcmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Wed, 21 Mar 2012 04:23:35 +0000 (04:23 +0000)
committercmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Wed, 21 Mar 2012 04:23:35 +0000 (04:23 +0000)
forester/ruby/scripts/preprocess.rb patch | blob | history

diff --git a/forester/ruby/scripts/preprocess.rb b/forester/ruby/scripts/preprocess.rb
index f3447ef..8829eb6 100644 (file)
--- a/forester/ruby/scripts/preprocess.rb
+++ b/forester/ruby/scripts/preprocess.rb
@@ -1,9 +1,16 @@
 #!/usr/local/bin/ruby -w
 
 
+if ( ARGV == nil || ARGV.length != 3 )
+  puts( "usage: preprocess.rb <input> <path to Pfam A HMMs>" )         
+  exit( -1 )
+end 
+
 input = ARGV[ 0 ]
 
-system( "hmmscan --nobias --domtblout " + input + "_hmmscan_260_10 -E 10 Pfam-A.hmm " + input + ".ni.fasta" )
+pfam = ARGV[ 1 ]
+
+system( "hmmscan --nobias --domtblout " + input + "_hmmscan_260_10 -E 10 " + pfam +  " + input + ".ni.fasta" )
 
 system( "hsp " + input + "_hmmscan_260_10 " + input + "_hmmscan_260_10_domain_table" )