Raj at 11:12, 20 May 2016

2016-05-20T11:12:04Z

← Older revision		Revision as of 11:12, 20 May 2016
Line 1:		Line 1:
	<~~syntaxhighlight~~ lang="perl">		<source lang="perl">

	#! /usr/bin/perl -w		#! /usr/bin/perl -w
Line 87:		Line 87:

	# END		# END
	</~~syntaxhighlight~~>		</source>

Raj at 07:01, 23 April 2011

2011-04-23T07:01:12Z

← Older revision		Revision as of 07:01, 23 April 2011
Line 1:		Line 1:
	<~~source~~ lang=perl>		<syntaxhighlight lang="perl">

	#! /usr/bin/perl -w		#! /usr/bin/perl -w

Line 86:		Line 87:

	# END		# END
	</~~source~~>		</syntaxhighlight>

Raj at 14:17, 26 December 2009

2009-12-26T14:17:48Z

← Older revision		Revision as of 14:17, 26 December 2009
Line 83:		Line 83:
	}		}

	print "\~~nDone~~ All Man\n\n";		print "\n\tDone All Man\nOutfile written into $resfile\n\n";
	~~#Removing Duplicates~~


	# END		# END
	</source>		</source>

Raj: New page: #! /usr/bin/perl -w # ************************************************************* # This Programme does the following in sequence # -> Open the misa results # -> Par...

2009-12-26T13:43:07Z

New page: <source lang=perl> #! /usr/bin/perl -w # ************************************************************* # This Programme does the following in sequence # -> Open the misa results # -> Par...

New page

<source lang=perl>
#! /usr/bin/perl -w

# *************************************************************
# This Programme does the following in sequence
# -> Open the misa results
# -> Parse the files
# -> compare the original EST file and pick the sequence that have repeats
# -> Outputs as fasta
# Author: Rajkumar (itc@rajkumar.in)

# *************************************************************

# Libraries
use strict;
use Bio::SeqIO;
use Bio::Perl;

# Variables
my $path = "/home/raj/bio";
my $ssrdir = "$path/euca/ssr";
my $datadir = "$ssrdir/data";
my $resdir = "$ssrdir/result";

my $est = "$datadir/globulus_est_trim.fa";
my $repeat = "$resdir/globulus_est_trim.fa.misa";
my $resfile = "$datadir/globulus_est_repeat.fa";

my ($i, $ssr_id, $total, $math ) = $_;
my (@repeat_id, @identifier) = @_;

# Function

# Lets Clean it
if (-e "$resfile") {`rm $resfile`;}

# start your hack

open(CSV, "$repeat") or die "Idiot! Where is that Misa Outfile ? 'CSV': $!\n";
<CSV>; # skip header line
while (my $line = <CSV>) {

my ($ssr_id) = split('\t', $line);
next if ($ssr_id =~ m/^ID/);

push (@repeat_id, "$ssr_id");
}

# Making them Unique
my %seen;
my @unique_repeat_id = grep !$seen{$_}++, @repeat_id;

# Getting some statistics
my $count_unique_repeat_id = @unique_repeat_id;
$i = 0;

# Lets Grab the sequence
foreach my $repeat_id(@unique_repeat_id) {

$i++;
my $in = Bio::SeqIO->new(-file => "$est", -format => 'Fasta');
my $identifier;
my $sequence;

while ( my $seq = $in->next_seq() ) {
$identifier = $seq->id;
$sequence = $seq->seq;

next if ($repeat_id ne $identifier);

if ($repeat_id eq $identifier) {

my $fasta = ">$identifier\n$sequence\n";
open (FH, ">>".$resfile);
print FH $fasta ;
close FH ;
print "\r\tWritten $identifier. Still $i of $count_unique_repeat_id to be written!";
}
}
}

print "\nDone All Man\n\n";
#Removing Duplicates

# END
</source>

Parse misa.pl - Revision history

Raj at 11:12, 20 May 2016

Raj at 07:01, 23 April 2011

Raj at 14:17, 26 December 2009

Raj: New page: #! /usr/bin/perl -w # ************************************************************* # This Programme does the following in sequence # -> Open the misa results # -> Par...