#!/usr/bin/env perl

use strict;
use warnings;

my @k_vals = (75, 77, 79, 81, 83, 85, 87, 89);

print '#!/bin/bash', "\n\n", ;

foreach my $k_val (@k_vals) { 
    print "    velveth Csp9_Berkeley_19apr2012_k$k_val $k_val", ' -create_binary -fastq.gz -short unpaired/*.fq.gz -shortPaired paired/*.fq.gz ; ', "\n",
          "    program_done_e-ping.pl -p velveth_Csp9_Berkeley_17apr2012_k$k_val ; \n",
          "\n",
          "    velvetg Csp9_Berkeley_19apr2012_k$k_val -cov_cutoff auto -exp_cov auto -min_contig_lgth 200 -ins_length 335 ; \n",
          "    program_done_e-ping.pl -p velvet_assembly_Csp9_Berkeley_17apr2012_k$k_val ; \n",
          "\n",
          ;
}