forester/ruby/scripts/preprocess.rb

   1 #!/usr/local/bin/ruby -w
   2
   3
   4 input = ARGV[ 0 ]
   5
   6 system( "hmmscan --nobias --domtblout " + input + "_hmmscan_260_10 -E 10 Pfam-A.hmm " + input + ".ni.fasta" )
   7
   8 system( "hsp " + input + "_hmmscan_260_10 " + input + "_hmmscan_260_10_domain_table" )
   9
  10 system( "d2f " + input + "_hmmscan_260_10_domain_table " + input + ".ni.fasta " + input + "_hmmscan_260_10.dff" )