annotate mapper_wrapper.pl @ 0:17e442abb3de draft default tip

Uploaded
author david-hoover
date Wed, 23 Jul 2014 10:27:07 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
1 #!/usr/bin/perl
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
2
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
3 $ENV{PATH} .= ":/usr/local/apps/mirdeep2/2.0.0.5:/usr/local/apps/bowtie/0.12.9:/usr/local/apps/viennarna/current/bin:/usr/local/randfold-2.0/bin";
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
4
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
5 $bowtie_index_name = $ARGV[0];
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
6 $output_arf = $ARGV[1];
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
7 $bowtie_files_path = $ARGV[2];
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
8 $options = join (" ", @ARGV[3..$#ARGV]);
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
9
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
10 $basename = `basename $bowtie_index_name`;
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
11 chomp $basename;
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
12
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
13 system ("mapper.pl $options -p $bowtie_files_path/$basename -t $output_arf 2> /dev/null");
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
14
17e442abb3de Uploaded
david-hoover
parents:
diff changeset
15 if (-s $output_arf == 0) {die "No reads aligned to the reference.";}