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" %i...2286 days ago
2266 days ago
Plot custom gene density with R
library(karyoploteR) pp2251 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_read_2.fastq #-t -treads=unit32 Number of tread...2248 days ago
2211 days ago
2211 days ago
Perl script to find coding regions in DNA sequences
#!/usr/bin/perl -w use strict; # if the number of input arguments is lower th...# extract the two values or columns with a regular expression # A group of letters and a decimal number ($c...2157 days ago
Perl script to find the distance beetween all the contigs and scaffolds
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; $| = 1; #Script to see the distance beetween all the contigs and scaffolds #Usage: perl clustalReads.pl genome....2142 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...reference multifasta file name here my $seqfile=""; my $queryfile = "genome.fasta"; # Ur query genome...2142 days ago
Perl script to convert GFF 2 FASTA !
#!/usr/bin/perl use strict; use warnings; use Bio::Seq; use Bio::SeqIO; use Bio::D...ly including the promoter) ### First, index the genome my $file_fasta = $ARGV[0];...y $db = Bio::DB::Fasta->new($file_fasta); print ("Genome...2141 days ago