Blast result parser with Perl and Bioperl
...in->next_result ) { # the name of the query sequence print OUT $result->query.... "\t"; # the length of the query sequence print OUT $result->quer...# get the start and the end of the query sequence in the alignment print O...2907 days ago
2906 days ago
Perl to print indivisual nucleotide from a sequence!
#!/usr/bin/perl use strict; use warnings; my $string = "ATGCTTGCGT?AAATG??CT?GCGTA"; my @chars = split("", $string); print "First character: $chars[0]\n";2906 days ago
2902 days ago
Needleman-Wunsch Algorithm in Perl
...erl nw.pl HEAGAWGHEE PAWHEAE BLOSUM50.txt -8 # See: "Biological sequence anaysis" Durbin et al. ed. CU...# usage statement die "usage: $0 \n" unless @ARGV == 4; # get sequences, matrix and gapcost from com...2898 days ago
Count GC Content in nucleotide sequence with Perl
#!/usr/bin/perl -w ### Usage: get_gc_content.pl ### #------------------------------------------------------...2898 days ago
Implementation of biological random mutation with Perl
#!/usr/bin/perl -w use strict; use warnings; #sequence for a better recognition my $DNA="AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA \ n"; my $I; my $...2894 days ago
Perl script to generate a random psuedo DNA sequence !
#!/usr/bin/perl print "Enter a number of nucleotides: \n"; chomp ($N = ); @b=qw/A T G C/;print ">Genome\n";while($l2875 days ago
Perl script to extract fasta sequence by matching name/ids !!
#!/usr/bin/perl use strict; use warnings; use Text::Trim qw(trim); #Usage perl extractSeqbyID.pl ids.txt seq.fasta Result.fasta $ARGV[2] or die "use extractSe...2867 days ago
Perl script introduces control structures, arrays and hashes.
..., 'Number of seqs' => 2 ); print $sequence{'DNA'}, "\n"; # Control Loop: foreach foreach my $key (sort (keys %sequence)) { print "Key of hash i...2854 days ago