Retrieve NCBI GenBank records with a range of accession numbers
...else { message( $param{verbose}, "Retrieving $count records out of $count available records.\n" ); } my $retmax = 500; if ( $retmax > $count...2930 days ago
A multilayer perceptron (MLP) neural network in Perl
...#To include an input bias create an #extra input in the training data #and set to 1 ################ User settings ######### my $numEpochs = 500; my $numHidden = 4; my $LR...2930 days ago
Needleman-Wunsch Algorithm in Perl
...# USAGE: perl nw.pl HEAGAWGHEE PAWHEAE BLOSUM50.txt -8 # See: "Biologi...SMATCH scores we will use values read from BLOSUM50) my $MATCH = 1; # +1 for...y $GAP = $gapcost; # for any gap my %BLOSUM50 = (); my @aalist = (); #...2920 days ago
2719 days ago
2708 days ago
picard tools command to get some insert statistics
...e=8g -jar $PICARD_HOME/picard.jar CollectInsertSizeMetrics MINIMUM_PCT=0 HISTOGRAM_FILE=.pdf INPUT= OUTPUT=.sorted.collectInseSize HISTOGRAM_WIDTH=5001352 days ago
2509 days ago
Create genome scaffolding with Perl
...ontig file for query sequences "prefix" => "psl_scaffold_", # prefix for contig names "pid" => 90, # percent ID threshold "trimlimit" => 50, # max number of overlapping...2323 days ago
Plot the clock using Lastz -gerenal outfile
use strict; use warnings; use Statistics::R ; use List::Util qw(sum); #Usage perl clockPlot.pl Palindrome.palfc 1500 my $R = Statistics::R->new() ; $R->startR ; my $file...2313 days ago
Plot the density of genes in R
#column1 = chromosome name and column2 = start position of the gene # check if ggplot2 is installed, if so, load it, # if not, install and load it if("ggplot2"...2299 days ago