Raku script to find SSRs in fastq file !
...return @ssrs; } sub process-fastq-file(Str $filename) { my $fh = open $filename, :r; my $line-number = 0; while $fh.readline -> $header { $line-number++; my $sequence = $f...126 days ago
Raku script to find overlaps between two bed files !
#!/usr/bin/env raku # Check if the correct number of arguments are provided if @*ARGS.elems != 2 { say "Usage: ./compare_bed_files.raku file1.bed file2.bed"; exit 1; }...108 days ago
Perl script to find overlaps between two bed files !
#!/usr/bin/perl use strict; use warnings; # Check if the correct number of arguments are provided if (@ARGV != 2) { die "Usage: $0 file1.bed file2.bed\n"; } # Read the contents of the two BED files my $file1 = shift @ARGV; my $file2 = shift @ARGV; open my $fh1, '108 days ago
Perl script to calculate the basic stats of the assembled genome !
...cs and information print "Genome Assembly Statistics:\n"; print "---------------------------\n"; print "Total Length: $total_length\n"; print "Number of Contigs: $num_contigs\n";...108 days ago
Python script for basic stats of the assembled genome !
...stics and information print("Genome Assembly Statistics:") print("---------------------------") print(f"Total Length: {total_length}") print(f"Number of Contigs: {num_contigs}")...108 days ago