BOL: Related items

Bioinformatics 101 - Running BLAST

Tue, 03 Sep 2013 14:59:50 -0500

How to format the database for BLAST, run the command, view the output file, and use BioPerl and Perl to parse the output. By David Francis, Ohio State University. Delivered live at the Tomato Disease Workshop 2010. For more information, please visit http://www.extension.org/pages/32521/bioinformatics-101-video.

Bioinformatics Scripts

Jit — Thu, 22 Jan 2015 22:29:39 -0600

Some of the useful bioinformatics scripts.

For example ... contig-stats.pl is a Perl script that will automatically describe features of a sequence assembly.

http://milkweedgenome.org/?q=scripts

Address of the bookmark: http://milkweedgenome.org/?q=scripts

Reverse Complement Problem Solved with Perl

Jit — Tue, 09 Jun 2015 23:37:23 -0500

Question at http://rosalind.info/problems/1b/

#Find the reverse complement of a DNA string.
#Given: A DNA string Pattern.
#Return: Pattern, the reverse complement of Pattern.

use strict;
use warnings;

my $string="AAAACCCGGT";
my $finalString="";
my %hash = (
   "C" => "G",
   "A" => "T",
   "T" => "A",
   "G" => "C",
);

for (my $aa=0; $aa<=(length($string)-1); $aa++) {
   my $char=substr $string, $aa, 1;
   #print $hash{$char};
   $finalString="$hash{$char}"."$finalString";
}

print $finalString;
print "\n";

Clump Finding Problem Solved with Perl

Jit — Wed, 10 Jun 2015 00:17:17 -0500

The question at http://rosalind.info/problems/1d/

Script are moved to http://bioinformaticsonline.com/snippets/view/34633/clump-finding-problem-solved-with-perl

CSBB-v1.0

Neel — Wed, 29 Jun 2016 07:33:05 -0500

CSBB is a command line based bioinformatics suite to analyze biological data acquired through varied avenues of biological experiments. CSBB is implemented in Perl, while it also leverages the use of R and python in background for specific modules. Major focus of CSBB is to allow users from biology and bioinformatics community, to get benefited by performing down-stream analysis tasks while eliminating the need to write programming code. CSBB is currently available on Linux, UNIX, MAC OS and Windows platforms.

Currently CSBB provides 13 modules focused on analytical tasks like performing upper-quantile normalization on expression data or convert genome wide gene expression to z-scores when comparing expression data from different platforms.

More at https://github.com/skygenomics/CSBB-v1.0

Address of the bookmark: https://github.com/skygenomics/CSBB-v1.0

SNPGenie

Jit — Thu, 30 Mar 2017 17:38:02 -0500

SNPGenie is a Perl script for estimating evolutionary parameters, mainly from pooled next-generation sequencing (NGS) single-nucleotide polymorphism (SNP) variant data. SNP reports (acceptable in a variety of formats) much each correspond to a single population, with variants called relative to a single reference sequence (one sequence in one FASTA file). Just run the main script, snpgenie.pl, in a directory containing the necessary input files, and we take care of the rest! For the earlier version, see Hughes Lab Bioinformatics Resource.

Address of the bookmark: https://github.com/hugheslab/snpgenie

RestrictionDigest: A powerful Perl module for simulating genomic restriction digests

Neel — Tue, 12 Jun 2018 13:17:13 -0500

RestrictionDigest can simulate the reference genome digestion and generate comprehensive information of the simulation. It can simulate single-enzyme digestion, double-enzyme digestion and size selection process. It can also analyze multiple genomes at one run and generates concise comparison of enzyme(s) performance across the genomes. For more information, please see the academic paper published online (http://www.sciencedirect.com/science/article/pii/S071734581630001X).

Address of the bookmark: https://github.com/JINPENG-WANG/RestrictionDigest

Benchmarking Perl Module !

Rahul Nayak — Sat, 25 Aug 2018 11:40:42 -0500

The benchmark module is a great tool to know the time the code takes to run. The output is usually in terms of CPU time. This module provides us with a way to optimize our code. With the advent of petascale computing and other multicore processor it is becoming a neccesity to know about the CPU time taken by our perl program.

This is the simple way to use the module

Example1:
use Benchmark;
$first_time = Benchmark->new;
our code……
$second_time = Benchmark->new;
$final_difference = timediff($first_time,$second_time);
print “the code took, timestr($final_difference),”\n”;

that was a very simple way to know the time diff , we can use it to know the time taken by some part of the code in the program.

More sophisticated way:
use Benchmark;
sub first {
my(arguments) = @_;
}
timethese(100, { first => ‘first_sub(arguments)’});
The first argument to timethese is 100 (evaluate 100 times).

Hope this very small tutorial with Benchmark will help people get started.

Machine learning in Perl

Jit — Sun, 16 Feb 2020 15:32:03 -0600

this is a fourth blog post in the Machine learning in Perl series, focusing on the AI::MXNet, a Perl interface to Apache MXNet, a modern and powerful machine learning library.

If you're interested in refreshing your memory or just new to the series, please check previous entries over here: 1 2 3

https://metacpan.org/pod/AI::MXNet

Address of the bookmark: http://blogs.perl.org/users/sergey_kolychev/2018/07/machine-learning-in-perl-kyuubi-goes-to-a-modelzoo-during-the-starry-night.html

Download mutliple fasta file from NCBI in one GO!!

Rahul Agarwal — Wed, 21 Aug 2013 08:13:30 -0500

if you have less time, then use three ways mentioned in bookmark link to extract/download all fasta sequences in single click given that you already have a list of GIs or accession IDs .

Alternatively, use one liner perl script:

perl -ne 'if(/^>(\S+)/){$c=$i{$1}}$c?print:chomp;$i{$_}=1 if @ARGV' GIs.txt >sequence.fasta

where GIs.txt contains a list of GIs or accession IDs.

(from :http://edwards.sdsu.edu/labsite/index.php/robert?start=5)

Address of the bookmark: http://edwards.sdsu.edu/labsite/index.php/robert/380-ncbi-sequence-or-fasta-batch-download-using-entrez