Perl script to split fasta sequence / overlaps
#!/usr/bin/perl use strict; use warnings; my $len = 5000; my $over = 200; my $seq_id=$ARGV[0]; my $seqFile = $ARGV[1]; my $seq; open(my $fh, "2041 days ago
Perl script to split fasta sequence and create overlaps
#!/usr/bin/perl use strict; use warnings; my $len = 5000; my $over = 200; my $seq_id=$ARGV[0]; my $seqFile = $ARGV[1]; my $seq; open(my $fh, "1970 days ago
1723 days ago
Perl script to run in parellel !
...allel::ForkManager; use Bio::SeqIO; my ($sequence_data_ref) = parse_genome...{name}, $genome{$chr_set}{nuc_seq}); } }...nce { my ($chr_set, $chr_seq) = @_; print "starting X...length; my $GCcount = $sequence =~ tr/GC|gc//; my...1684 days ago
1555 days ago
Perl onliner to print fasta headers !
#Save all your fasta in seq.fa and run the following ... perl -ne 'print if /^>/' seq.fa #Print header with line number perl -ne 'print "$. $_" if /^>/ ' seq.fa994 days ago
Bash script to simulate a genome !
# Reference https://github.com/chhylp123/hif...activate bcftools1.10.2 # Use Seqtk to convert soft-masked base...s with bgzip # https://github.com/lh3/seqtk /genetics/elbers/bin/seqtk...ploid.fasta.gz|\ /genetics/elbers/bin/seqtk/seqtk seq -L0|paste - - |gr...916 days ago
Perl script to find inverted repeats !
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; use Bio::Tools::Run::Repeat...genome.fasta"; # read genome sequence my $seqio = Bio::SeqIO-...file, -format => "fasta"); my $seqobj = $seqio->next_seq(); my...q "Inverted") { my $rm_seq = substr($seq, $rm_start-1, $...423 days ago
R script to covert and export html page to png
# Library library(streamgraph) # Create data: data416 days ago
Perl script for chi-squared test !
#!/usr/bin/perl # # chidi.pl # # A scri...############################## # Read sequence file 1 ################...FAlite(\*FILE); # loop through each sequence in file 1 while(my $ent...y $entry = $fasta->nextEntry) { my $seq = uc($entry->seq); # to cou...409 days ago