2880 days ago
2880 days ago
Reformat the file names with Perl
...use strict; use warnings; use File::Copy qw(copy);; $| = 1; my %hash; my @files = glob "*.scf"; if (!$ARGV[0]){ print "Table file needed\n USAGE: $0 table\n"; exit;} my $ifh = read_...2344 days ago
Plot the density of genes in R
#column1 = chromosome name and column2 = start position of the gene # check if ggplot2 is installed, if so, load it, # if not, install and load it if("ggplot2" %...2293 days ago
2273 days ago
Plot custom gene density with R
library(karyoploteR) pp2258 days ago
Estimate Genome Size with Jellyfish and R
jellyfish count -t 8 -C -m 19 -s 5G -o 19mer_out --min-qual-char=? /common/Tutorial/Genome_estimation/sample_read_1.fastq /common/Tutorial/Genome_estimation/sample_rea...2254 days ago
Perl script to find coding regions in DNA sequences
...the filenames my $filecodontable = $ARGV[0]; my $filesequence...loglkh.pl: the file $filecodontable can not be opened\n"; exit...\d+)/); # load the hash table $pcodons{$codon}=$freq; }...probability computed from the table of codon frequencies divided...2163 days ago
2135 days ago
2084 days ago