Download genome using NCBI esearch command line !
#Download using command line "NC_005956.1" / "CP003784.1" esearch -db nucleotide -query "NC_005956.1" | efetch -format fasta > seq1.fa esearch -db nucleotide -query "...1116 days ago
Blast result parser with Perl and Bioperl
#!/usr/local/bin/perl # # Dr. Xiaodong Bai # It may be freely distributed under GNU General Public License. # This script will parse a NCBI blastx output file and ou...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-Wun...2920 days ago
Calculate some statistics for a DNA alignment with Perl
use Bio::AlignIO; use Bio::Align::DNAStatistics; my $stats = Bio::Align::DNAStatistics->new(); my $alignin = Bio::AlignIO->new(-format => 'emboss',...2664 days ago
Mapping with BWA-mem or BWA-sampe in one go with python script !
BAM files and mapping BESST requires sorted and indexed BAM files as input. Any read aligner + samtools can be used to obtain such files. Read pairs needs to be aligned i...2139 days ago
Biological Sequence handling with Perl !
package Sequence::Generic; # File: Sequence/Generic.pm use strict; use Carp; use overload '""' => 'asString', 'neg' => 'reverse', '.'...2196 days ago
Perl script to run SATSUMA in loop !
#!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage perl 1by1.pl for SATSUMA analysis # User need to set the reference multifasta file name here m...2155 days ago
2141 days ago
Install Parrot Virtual Machine !
#Parrot is a virtual machine designed to efficiently compile and execute bytecode for dynamic languages jit@jit-HP-Pro-3335-MT:~/Downloads/parrot-8.1.0$ wget ftp://ftp....1552 days ago
Perl script to merge LastZ overlaps
#!/usr/bin/perl use strict; use warnings; use 5.010; # Filter out the exact/direct overlaps from tab seperated alignment file. (lastz format=general- ready) # Do not inclide header in lastz outfile # USAGE: perl filterOverlaps.pl infile > outfile open my $fh, '2140 days ago