diff count_ref_length.pl @ 50:7b5a48b972e9 draft

Uploaded

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/count_ref_length.pl	Fri Dec 05 00:11:02 2014 -0500
@@ -0,0 +1,58 @@
+#!/usr/bin/perl -w
+#Filename:
+#Author: Tian Dongmei
+#Email: [e-mail]
+#Date: 2009-05-06
+#Modified:
+#Description: ɾ³ýmatched reads 
+my $version=1.00;
+
+use strict;
+use Getopt::Long;
+
+my %opts;
+GetOptions(\%opts,"i=s","o=s","h");
+if (!(defined $opts{i} and defined $opts{o} ) || defined $opts{h}) { #necessary arguments
+&usage;
+}
+
+my $filein=$opts{'i'};
+my $fileout=$opts{'o'};
+
+open IN,"<$filein"; #input file  
+open OUT,">$fileout"; #output file  
+
+my ($name,$seq);
+while (my $aline=<IN>) {
+	chomp $aline;
+	if ($aline=~/^>(\S+)/) {
+		$name=$1;
+		while (my $new=<IN>) {
+			chomp $new;
+			if ($new=~/^>(\S+)/) {
+				print OUT $name,"\t",length($seq),"\n";
+				$seq ="";
+				$name=$1;
+				next;
+			}
+			else{$seq .=$new;}
+		}
+	}
+	print OUT $name,"\t",length($seq),"\n";
+}
+ 
+close IN;
+close OUT;
+sub usage{
+print <<"USAGE";
+Version $version
+Usage:
+$0 -i -o
+options:
+-i input file
+-o output file
+-h help
+USAGE
+exit(1);
+}
+