diff Tools/CREF/Extract_motif_codes_galaxy.pl @ 0:229d36377838 draft

Uploaded
author amadeo
date Mon, 05 Sep 2016 05:53:08 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Tools/CREF/Extract_motif_codes_galaxy.pl	Mon Sep 05 05:53:08 2016 -0400
@@ -0,0 +1,54 @@
+#!/usr/bin/perl -w
+
+# The lines that start with # are comment lines that are not executed
+
+
+$|=1;
+use strict; 
+use warnings;
+
+
+#Script to take output from AME (part of memesuite-org) and extract a
+#list of the overrepresented motifs and print them to a new file
+#called ame-motif-id.list
+
+
+#Checking to see if the user has provided 1 argument - which is the
+#name of the AME results file
+
+if(@ARGV < 2){
+print "\nUsage: Extract_motif_codes.pl ame.txt ame-shorted.txt\n\n";
+exit(0);
+}
+
+#Declaring variables
+my @cols; #an array variable
+my $line; # a scalar varaible
+
+#Using a FIELHANDLE to open the input file
+open (INPUT, "<$ARGV[0]") || 
+   die "File '$ARGV[0]' not found\n" ;
+   
+open (OUTPUT, ">$ARGV[1]") || 
+  die "File '>$ARGV[1]' not found\n" ;
+
+#looping through each line of the file
+ while (<INPUT>){
+   #assigning line to variable $line 
+   #$_ is a special default variable that here holds the line contents
+   $line = $_;
+   #match lines that have Ranksum 
+   if ($line =~ /Ranksum/){
+    printf OUTPUT "%s\n", $line; 
+    #split the lines on white space, so each part of the line gets
+    #stored as an array element
+    @cols=split;
+    #Testing to see what line elements are stored in the array
+    #print "cols [0] is $cols[0] \n";
+    #print "cols [2] is $cols[2] \n\n";
+    
+    #Now see if you can print out the array elemnent that stores the
+    #motif ID to a new file called ame-motif-id.list.
+  }
+ }
+