Retrieve NCBI GenBank records with a range of accession numbers
...es} >= $param{maxRetries} ) { die("Too many failures--giving up search."); }...>= $param{maxRetries} ) { die("Too many failures--giving up search."); }...2929 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 # Needleman-Wunsch global...2918 days ago
Count GC Content in nucleotide sequence with Perl
...Count\n"; $seq = ""; while () { chomp; if (/^>/) { #finish up previous line. if (length($...OUT "$id\t"; } else { $seq = $seq . $_; } } #finish up last line. &process_it; c...2918 days ago
2884 days ago
2706 days ago
Perl script to remove the duplicate sequences from multifasta file
use strict; use Bio::SeqIO; my %unique; my $file = "myseqs.fa"; my $seqio = Bio::SeqIO->new(-file => $file, -format => "fasta"); my $outseq = Bio::SeqIO->ne...2703 days ago
Perl script to insert the DNA string in genome
...ry_id eq $chr){ my $length = length($seq->seq); my $upstream=substr($seq->seq, 0, $p...seq, $pos,$length); my $seq_obj = Bio::Seq->new(-seq => $upstream.$seqI.$downstream,-disp...2672 days ago
Calculate ATGC percentage in parallel with perl
...es; # hash to resolve PID's back to child specific information my $pm = new Parallel::ForkManager($max_procs); # Setup a callback for when a child finishes up so we can # get it's exit...2669 days ago
2649 days ago
Genetic Algorithms demonstration with word DNA in Perl
...chomp; } # eliminate words over $max_entry_length letters, and uppercase them @dictionary =...onary; # build the letter frequencies hash (remember, all letters are uppercase) $freqs{$_}++ fore...2384 days ago