2931 days ago
2931 days ago
Parse a genbank file using regular expressions
...*)/) { print "Locus: $2\n"; } elsif (/(VERSION.*GI:)(\d*)/) { print "GI: $2\n"; } elsif (/(DEFINITION\s*)(.*)(\.)/)...print "CDS: $cds_start - $cds_end\n"; } elsif (/(\/translation=")(.*)/)...2931 days ago
Perl script to count the number of Adenine, Thymine, Guanine and Cytosine in your DNA Sequence
...nt_of_T = 0; $errors = 0; foreach $base (@DNA) { if ( $base eq 'a' ) { ++$count_of_A; } elsif ( $base eq 'c' ) { ++$count_of_C; } elsif ( $base eq 'g' ) {...2930 days ago
Retrieve NCBI GenBank records with a range of accession numbers
...e => undef, url => 'http://www.ncbi.nlm.nih.gov/entrez/eutils', retries => 0,...); $count = $param{maxRecords}; } else { message( $par...2930 days ago
A multilayer perceptron (MLP) neural network in Perl
#!/usr/local/bin/perl -w #################################################### #MLP neural network in Perl Original source code by Phil Brierley #Translated into...2930 days ago
Blast result parser with Perl and Bioperl
...ery_length; # output "no hits found" if there is no hits if ( $result->num_hits == 0 ) { print OUT "\tNo hits found\n"; } else { my $count = 0;...2929 days ago
2929 days ago
Find the number of each 2 consecutive characters AA, AC,AG,AT,CC,CA... with Perl
...ubject =~ m/[ACTG][ATGC]/g) { # matched text = $& if(exists $results{$&}) { $results{$&}++ } else { $res...2924 days ago
Needleman-Wunsch Algorithm in Perl
#!/usr/bin/perl # USAGE: perl nw.pl HEAGAWGHEE PAWHEAE BLOSUM50.txt -8 # See: "Biological sequence anaysis" Durbin et al. ed. CUP 1998, Pg. 19 # Needlema...2920 days ago