Perl script to count the number of Adenine, Thymine, Guanine and Cytosine in your DNA Sequence
...count_of_A = 0; $count_of_C = 0; $count_of_G = 0; $count_of_T = 0; $errors = 0; foreach $base (@DN...eq 'G' ) { ++$count_of_G; } else { print "Error - Unknown base: $base\n";...3054 days ago
Retrieve NCBI GenBank records with a range of accession numbers
...($esearch_result =~ m/(.*)/is) { die("ESearch returned an error: $1"); } me...} open( my $OUTFILE, ">" . $param{outputFile} ) or die("Error: Cannot open $param{outputFil...3054 days ago
A multilayer perceptron (MLP) neural network in Perl
...###### my $numInputs = 3; my $numPatterns = 4; ######################################## my $patNum; my $errThisPat; my $outPred; my $RMSerror; my @trainInputs; my @tra...3054 days ago
Perl program to implement sliding window !
#!/usr/bin/perl -w my $filename = 'data.txt'; open(my TR, '3053 days ago
3048 days ago
Perl script to extract lines with matching ids !!
#!/usr/bin/perl use strict; use warnings; my %patterns; #USAGE: perl extactByIds.pl Idsfile1 file2 > Result # Open file and get patterns to search for open(my $fh2,"3013 days ago
3010 days ago
Create genome scaffolding with Perl
...; } else { pod2usage({-exitVal => 1, -message => "Error: Unknown command-line option...ts->{"query"}){ pod2usage({-exitVal => 1, -message => "Error: No query assembly file provi...2447 days ago
Perl script to remove fasta sequences in multifasta file with certain length threshold
#!/usr/bin/perl use strict; use warnings; my $minlen = shift or die "Error: `minlen` parameter not provided\n"; { local $/=">"; while() { chomp; ne...2390 days ago
Plot custom gene density with R
library(karyoploteR) pp2388 days ago