Download the genome from NCBI using bash script/command
...mkdir "G$filename" cd "G$filename" # take action on each file. $f store current file name wget --input $f #cat $f cd .. done #Reference #ftp://ftp.ncbi.nlm.nih.gov/...2525 days ago
Genetic Algorithms demonstration with word DNA in Perl
...my $pop_ref = []; # a reference to a population array...$index = sample(\@weights); # we pass a reference to the weights array here...PAN Sample module) sub sample { # get the reference to the weights array my $we...2365 days ago
2203 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 my $seqfile=""...2133 days ago
2119 days ago
2116 days ago
2068 days ago
Perl script to create a consensus of nucleotide sequences !
...@cons, [ sort { ($H{$b}->[$B]||0) ($H{$a}->[$B]||0) || $prefOrder{$b} $prefOrder{$a} } qw/A T G C/ ]->[0]; } return @cons; } #reference https://www.perlmonks.org/bar...2025 days ago
1718 days ago
Samtools commands for bioinformatician !
...command) # -g : output is bcf (binary call format) file # -f : use reference genome given samtools mpileu...sal_vars.bcf > sal_vars.vcf ### normalize (realign) indels # -f : reference fasta, needed to left align a...1596 days ago