Perl script to break the contigs by 'N'
#!/usr/bin/perl -w use Bio::SeqIO; use strict; my $fasta = Bio::SeqIO->new( -file => "$ARGV[0].parts.fasta",-format=>'fasta'); open(SCAFF,">$ARGV[0].parts.scaff"); while ( my $se...2090 days ago
Perl script to run in parellel !
#!/usr/bin/perl use strict; use warnings; use Parallel::ForkManager; use Bio::SeqIO; my ($sequence_data_ref) = parse_genome_files($ARGV[0]);...$file=shift; my (%sequence_data); my $file_content = new Bio::Seq...1714 days ago
1541 days ago
Installing docker for Bioinformatics on Ubuntu !
...E sh -c echo "deb [arch=amd64] https://download.docker.com/linux/ubuntu bionic stable" > /etc/apt/sources...ge the docker engine image Manage images network Manage networks node Manage Swarm n...1541 days ago
1536 days ago
Python script to check sequence length in multifasta file
#!/usr/bin/python from Bio import SeqIO import sys cmdargs = str(sys.argv) for seq_record in SeqIO.parse(str(sys.argv[1]), "fasta"): output_line = '%s\t%i' % \ (seq_record.id, len(seq_record)) print(output_line)1243 days ago
1243 days ago
Parse the NCBI taxonomy database with Perl !
use Bio::DB::Taxonomy; use warnings; my $sps="human"; # Get one from a NCBI taxonomy database my $dbh = Bio::DB::Taxonomy->new(-source...$nam); #print $ancestor_taxon; use Bio::Tree::Tree; my $tree_functi...{ my ($lineageId,$classNode)=@_; use Bio::Tre...1200 days ago
Perl script to get all the descendent in a tree !
use warnings; use Bio::DB::Taxonomy; #Rotifera 10190 # Get one from a NCBI taxonomy database my $dbh = Bio::DB:...1192 days ago
Perl script to check perl modules and download NCBI, BUSCO, Taonomy, Silva databases !
use strict; use warnings; use ExtUtils::Installed; use LWP::Simple; use Archive::..., 'Data::Dumper', 'File::Basename', 'Bio::SeqIO', 'FindBin',...s::R', 'Math::Round', 'File::Find', 'Bio::DB:...1188 days ago