Perl script to find inverted repeats !
#!/usr/bin/perl use strict; use warning...eqIO; use Bio::Tools::Run::RepeatMasker; my $genome_file =...= $seqobj->seq(); # run RepeatMasker my $rm = Bio::Tools:...n($genome_file); # parse RepeatMasker output while (my $rm...} print "Inverted repeat found at positions $rm_star...487 days ago
Perl script to find edit distance between two sequences !
#!/usr/bin/perl use strict; use warnings; sub edit_distance { my ($s1, $s2) = @_; my $len1 = length($s1); my $len2 = length($s2); my @dp; for (my $i = 0; $i487 days ago
Raku script to find palindrome in genomes !
sub is-palindrome(Str $str) returns Bool { $str.=uc; # convert to uppercase $str.=subst:g/\s+//; # remove any spaces return $str eq $str.flip; } sub find-palindro...487 days ago
R script to covert and export html page to png
# Library library(streamgraph) # Create data: data480 days ago
54 days ago
Perl script for chi-squared test !
#!/usr/bin/perl # # chidi.pl # # A script to perform a chi-squared test of the di...; # hashes for obersered and expected dinucleotide frequencies...############################# open(FILE,"$ARGV[0]") || die "Can...############################# open(FILE,"$ARGV[1]") || die "Can...473 days ago
Raku script to find SSRs in fastq file !
...for 2..$sequence.chars -> $min-repeats { for $sequence.chars...$min-repeats -> $max-repeat {...sequence.substr($min-repeats - 1, $max-repeat - $min-repeats + 1);...{ start => $min-repeats, end...174 days ago
Perl script to calculate GC content !
#!/usr/bin/perl sub calculate_gc_content { my ($sequence) = @_; $sequence = uc($sequence); # Convert the sequence to uppercase my $gc_count = () =...158 days ago
Python script to convert Multi-line Fasta to Single-line Fasta
def multi_to_single_line_fasta(input_filename, output_filename): try: with open(input_filename, 'r') as input_file: with open(output_filename, 'w') as output_file: cu...157 days ago
Perl script to convert Multi-line Fasta to Single-line Fasta !
#!/usr/bin/perl use strict; use warnings; sub multi_to_sin...y ($input_filename, $output_filename) = @_; open my $input_file, '', $output_filename or die "Error: Could not open file '$output_filename': $!"...157 days ago