Misa run.pl: Difference between revisions

From RAJ INFO
Jump to navigation Jump to search
m New page: <source lang=perl> #! /usr/bin/perl -w # ************************************************************* # This Programme does the following in sequence # -> Runs MISA # -> Produces 2 outfi...
 
mNo edit summary
 
Line 20: Line 20:
my $infile = "$path/euca/data/globulus/globulus_est_trim.fa";
my $infile = "$path/euca/data/globulus/globulus_est_trim.fa";
my $ssrdir  = "$path/euca/ssr";
my $ssrdir  = "$path/euca/ssr";
my $resdir  = "$ssrdir/result";


# Running the File
# Running the File

Latest revision as of 13:21, 26 December 2009

#! /usr/bin/perl -w

# *************************************************************
# This Programme does the following in sequence
# -> Runs MISA
# -> Produces 2 outfiles, one ssr and other statistics
# Author: Author: Thomas Thiel ; Modified by Rajkumar (itc@rajkumar.in)
# Release: JUL 2008
# *************************************************************

# Libraries
	use strict;

# At the top: Time Start
	my ($sec,$min,$hour,$mday,$mon,$year,$wday,$yday,$isdst) = localtime(time);

# Variables
	my $path    = "/home/raj/bio";
	my $infile = "$path/euca/data/globulus/globulus_est_trim.fa";
	my $ssrdir  = "$path/euca/ssr";

# Running the File
	chdir $ssrdir;
	`./misa.pl $infile`;


# At the End: Time Calculation
	my ($sec_,$min_,$hour_,$mday_,$mon_,$year_,$wday_,$yday_,$isdst_) = localtime(time);
	my $days = $mday_-$mday;
	my $hours = $hour_-$hour;
	my $mins = $min_-$min;
	my $secs = $sec_-$sec;
	print "\n------------------------\nDuration: ";
	if ($days>0) { print "$days days";}
	if ($hours>0) { print ": $hours hours";}
	if ($mins>0) { print ": $mins mins";}
	if ($secs>0) { print ": $secs secs";}
	print "\n------------------------\n";

# END