Perl script to split fasta sequence / overlaps
#!/usr/bin/perl use strict; use warnings; my $len = 5000; my $over = 200; my $seq_id=$ARGV[0]; my $seqFile = $ARGV[1]; my $seq; open(my $fh, "2097 days ago
Perl script to create a consensus of nucleotide sequences !
use strict; use warnings; my @instances = qw ( AAAAA ATCGA ATAAA ); my @instances2 = qw ( AAAAA AACGA...AAA exit; sub consensus{ my @mi = @_; chomp(@mi); my...[] ); s/\s//g for @mi; my ($w) = sort {$b $a} map {len...sh @{$H{G}}, $g; } my @cons = (); my %prefOrde...2086 days ago
2058 days ago
Perl script to split fasta sequence and create overlaps
#!/usr/bin/perl use strict; use warnings; my $len = 5000; my $over = 200; my $seq_id=$ARGV[0]; my $seqFile = $ARGV[1]; my $seq; open(my $fh, "2026 days ago
Perl script to count occurrence of a character !
#!/usr/bin/env perl # -*- coding: utf-8 -*- #!/usr/bin/perl use strict; use warnings; my %count_of; while ( ) { my @val = split "\t", $_; #my ( $word) = m/(\w+)/; $count_...1954 days ago
1854 days ago
Perl script to run in parellel !
...kManager; use Bio::SeqIO; my ($sequence_data_ref) = parse_genome_files($ARGV[0]); my %genome=%{$sequence_data_ref}...er->new( $n_processes ); for my $i ( 1 .. $n_processes ) {...sub parse_genome_files { my $file=shift; my (%sequen...1740 days ago
Perl subroutine to creating kmer !
sub k_mers { my ($sequence, $k) = @_; my $len = length($sequence); my @result = (); for (my $i = 0; $i1621 days ago
Perl script to remove duplicated lines !
#!/usr/bin/perl use strict; use warnings; { $_ = ; my $next_line; while( $next_line = ) { #print "current line: $_ -- next line: $next_line$/"; print $...1614 days ago
1613 days ago