2128 days ago
R script to visualize fastANI core-genome comparison
...ualize fastANI core-genome comparison # Usage: Rscript # Output: .pdf # Uses genoPlotR package: http://genoplotr.r-forge.r-project.org #Parse command line arguments query...1086 days ago
Parse a genbank file using regular expressions
#! /usr/local/bin/perl -w $genbank = "genbank_file.txt"; open (GENBANK, $genbank) || die "cannot open $gb_report for reading: $!"; # Flag for multiline transl...2901 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 output the to...2899 days ago
2898 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 her...2124 days ago
Perl script to convert GFF 2 FASTA !
...moter) ### First, index the genome my $file_fasta = $ARGV[0]; my $db = Bio::DB::Fasta->new($file_fasta); print ("Genome fasta parsed\n"); ### Second, parse the GFF3 my %CDS; my %CDNA;...2123 days ago
Perl script to extract sequence by Ids from multifasta file !
...cal $/ = "\n>"; # read by FASTA record open FASTA, $seqfile; while () { chomp; my $seq = $_; my ($id) = $seq =~ /^>*(\S+)/; # parse ID as first word in FASTA hea...2080 days ago
1988 days ago
1537 days ago