Codebase list cd-hit / 5ba8578
Merge pull request #86 from tseemann/fix_perl_hashbangs Set all hashbangs to /usr/bin/env perl Weizhong Li authored 4 years ago GitHub committed 4 years ago
59 changed file(s) with 61 addition(s) and 61 deletion(s). Raw diff Collapse all Expand all
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Storable;
33 use strict;
0 #!/usr/bin/perl -w
0 #!/usr/bin/env perl
11 # =============================================================================
22 # CD-HIT
33 # http://cd-hit.org/
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my $script_name = $0;
33 my $script_dir = $0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 #
22
33 my $rep;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #not like cd-hit-div, this script do not sort input
33 #or throw away seq
0 #!/usr/bin/perl -w
0 #!/usr/bin/env perl
11 # =============================================================================
22 # CD-HIT
33 # http://cd-hit.org/
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 $clstr = shift;
33 $fr = shift; # for nr80.clstr $fr = 0.8
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my $no = 0;
33 my $clstr_no = "";
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #usage: clstr_xml.pl [-len|-size] level1.clstr [level2.clstr level3.clstr ...]
33 #purpose: to create xml file from cd-hit or hierarchical cd-hit(h-cd-hit) results
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #keep only top $no proteins in cluster
33
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Storable;
33 use strict;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Storable;
33 use strict;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 # the order of clusters need to be identical
33 my ($master_clstr, @clstr) = @ARGV;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 # order of clusters don't need to be the same
33 # but then I have to read everything into memory
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 ## calculate the sensitivity and specificity of clusters
33 ## if the input fasta file has pre-defined classification term
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 ## calculate the sensitivity and specificity of clusters
33 ## if the input fasta file has pre-defined classification term
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22
33 $file90 = shift;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 $no = 0;
22 while($ll=<>){
33 if ($ll =~ /^>Cluster (\d+)/) {
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 $rep = "";
33 $no = 0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 # output single fasta file
22 # for each cluster output at least $cutoff seqs
33
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 # if nr90 from nr100 and
22 # nr80 from nr90, so I have nr90.clstr and nr80.clstr
33 # but, in nr80.clstr, some gi numbers whose from nr100 are there
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #my $by = shift;
33 my $min;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #my $by = shift;
33 my $min;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 if(@ARGV==0){
33 print "Usage:\n\tclstr_size_histogram.pl [-bin N] clstr_file\n";
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 if(@ARGV==0){
33 print "Usage:\n\tclstr_size_stat.pl clstr_file\n";
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my $sort_by_what = shift;
33 $sort_by_what = "no" unless $sort_by_what;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my $sort_by = shift;
33 $sort_by = "len" unless ($sort_by);
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 if(@ARGV==0){
33 print "Usage:\n\tclstr_sql_tbl.pl clstr_file tbl_file\n";
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 if(@ARGV==0){
33 print "Usage:\n\tclstr_sql_tbl_sort.pl table_file level\n";
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #note you have to use "-d 0" in the cd-hit run
33 #note you better to use "-g 1" in the cd-hit run
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Image::Magick;
33
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 $file90 = shift;
33 $segs = shift;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 #not like cd-hit-div, this script do not sort input
33 #or throw away seq
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my $by = shift;
33 my $min;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my $by = shift;
33 my $min;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my ($gi_file, $seq_file) = @ARGV;
33
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 my ($gi_file, $seq_file) = @ARGV;
33
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22
33 my $script_name = $0;
0 #!/usr/bin/perl -w
0 #!/usr/bin/env perl
11 ################################################################################
22 ######### PSI-cd-hit written by Weizhong Li at http://cd-hit.org
33 ################################################################################
11371137
11381138 open(REPERL, "> $remote_perl_script") || die;
11391139 print REPERL <<EOD;
1140 #!/usr/bin/perl
1140 #!/usr/bin/env perl
11411141 \$host = shift;
11421142 \$arg = shift;
11431143
0 #!/usr/bin/perl -w
0 #!/usr/bin/env perl
11 ################################################################################
22 ######### PSI-cd-hit written by Weizhong Li at http://cd-hit.org
33 ################################################################################
13711371
13721372 open(REPERL, "> $remote_perl_script") || die;
13731373 print REPERL <<EOD;
1374 #!/usr/bin/perl
1374 #!/usr/bin/env perl
13751375 \$host = shift;
13761376 \$instance = shift;
13771377 \$arg = shift;
0 #!/usr/bin/perl -w
0 #!/usr/bin/env perl
11 ################################################################################
22 ######### PSI-cd-hit written by Weizhong Li at http://cd-hit.org
33 ################################################################################
0 #!/usr/bin/perl -w
0 #!/usr/bin/env perl
11 ################################################################################
22 ######### PSI-cd-hit written by Weizhong Li at http://cd-hit.org
33 ################################################################################
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ## =========================== NGS tools ==========================================
22 ## NGS tools for metagenomic sequence analysis
33 ## May also be used for other type NGS data analysis
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ################################################################################
22 # NGS workflow by Weizhong Li, http://weizhongli-lab.org
33 ################################################################################
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 # =============================== NG-Omics-WF ==================================
22 # _ _ _____ ____ _ __ ________
33 # | \ | |/ ____| / __ \ (_) \ \ / / ____|
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ## =========================== NGS tools ==========================================
22 ## NGS tools for metagenomic sequence analysis
33 ## May also be used for other type NGS data analysis
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Getopt::Std;
33 my $script_name = $0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 #
22 use Getopt::Std;
33 getopts("i:s:S:o:f:j:b:h:",\%opts);
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Getopt::Std;
33 my $script_name = $0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Getopt::Std;
33 my $script_name = $0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Getopt::Std;
33 my $script_name = $0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Getopt::Std;
33 my $script_name = $0;
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ## =========================== NGS tools ==========================================
22 ## NGS tools for metagenomic sequence analysis
33 ## May also be used for other type NGS data analysis
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 #
22 use Getopt::Std;
33 getopts("s:S:o:f:j:",\%opts);
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ## =========================== NGS tools ==========================================
22 ## NGS tools for metagenomic sequence analysis
33 ## May also be used for other type NGS data analysis
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ## =========================== NGS tools ==========================================
22 ## NGS tools for metagenomic sequence analysis
33 ## May also be used for other type NGS data analysis
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 ################################################################################
22 # NGS workflow by Weizhong Li, http://weizhongli-lab.org
33 ################################################################################
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11 #
22 use Getopt::Std;
33 getopts("i:s:S:o:f:j:",\%opts);
0 #!/usr/bin/perl
0 #!/usr/bin/env perl
11
22 use Getopt::Std;
33 my $script_name = $0;