Plot the density of genes in R
...# check if ggplot2 is installed, if so, load it, # if not, install and load it if("ggplot2" %in% rownames(installed.packages())){ library(ggplot2) } else { install.packages("ggplot2") library...2328 days ago
2308 days ago
2296 days ago
Perl script to remove fasta sequences in multifasta file with certain length threshold
#!/usr/bin/perl use strict; use warnings; my $minlen = shift or die "Error: `minlen` parameter not provided\n"; { local $/=">"; while() { chomp; n...2295 days ago
Plot custom gene density with R
library(karyoploteR) pp2293 days ago
Perl script to find palindromic regions in DNA sequences
use strict; use warnings; my $pp = qr/(?: (\w) (?1) \g{-1} | \w? )/ix; my $filename = $ARGV[0]; open(my $fh, '2198 days ago
Biological Sequence handling with Perl !
package Sequence::Generic; # File: Sequen...e object" unless $new_seq->isa(__PACKAGE__); $to_append =...1; Back to Article Listing Two package Sequence::Nucleotide; #...my ($start,$end) = @_; return (__PACKAGE__ . '::Subsequence')->ne...2225 days ago
Mapping with BWA-mem or BWA-sampe in one go with python script !
...ny read aligner + samtools can be used to obtain such files. Read pairs needs to be aligned in paired read mode. BESST provides...ll for mapping with this script is python reads_to_ctg_map.py /path/to/lib1_A.fq /path/to/lib1_...2168 days ago
2195 days ago
2188 days ago