Insert the sequence at desire location in multi-fasta file with Perl
...omechrName locationStart AlienGene AlienLength # The coordinate should not overlaps --- next postition shold be bigger than firstpos+alienLen open(my $fh, '2310 days ago
Create genome scaffolding with Perl
...$projOpts->{"query"}), -verbose => 0}); } print(STDERR "Loading query sequences into memory..."); open(my $queryFile, "=21 tab-separ...2304 days ago
Plot the clock using Lastz -gerenal outfile
...clockPlot.pl Palindrome.palfc 1500 my $R = Statistics::R->new() ; $R->startR ; my $fileN=$ARGV[0]; my $mSize=$ARGV[1]; open (my $fh2, $fileN) or die "Could not open file $fileN $!"; my (@allKey...2295 days ago
Remove duplicate lines with perl
#! perl -sw use strict; my %lines; #open DATA, $ARGV[0] or die "Couldn't open $ARGV[0]: $!\n"; while () { print if not $lines{$_}++; } __DATA__ apple apple plum vinegar apple banana banana banana apple2288 days ago
846 days ago
Plot custom gene density with R
library(karyoploteR) pp2245 days ago
Perl script to read multi fasta sequence one by one
...le='tmp.fa'; my $fastaSeq_ref = readfasta ("$ARGV[0]"); my %fastaSeq = %$fastaSeq_ref; foreach my $key ( keys %fastaSeq) { open (OUT, ">$outfile") or die "couldn't open the file $outfile $!"; print...2217 days ago
Perl script to extract the uniq Ids
#!/usr/bin/perl -w use strict; use warnings; use List::Uniq ':all'; open(my $val, "2235 days ago
Perl script to find missing and move to desire folder
#!/usr/bin/perl -w use strict; use warnings; open(my $ids, "2235 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, '2150 days ago