Python script to convert Multi-line Fasta to Single-line Fasta
...ilename} in single-line FASTA format.") except FileNotFoundError: print(f"Error: File '{input_filename}' not found.") # Example usage: # multi_to_single_line_fast...139 days ago
Perl script to convert Multi-line Fasta to Single-line Fasta !
...close $output_file; print "Successfully converted $input_filename to $output_filename in single-line FASTA format.\n"; } # Example usage: # multi_to_single_line_fast...139 days ago
138 days ago
Raku script to find repeats in sequences !
...substring) > 1 && $substring !~~ any(@repeats) { @repeats.push($substring); } } return @repeats; } # Example usage my $genome-sequence = "ATCGA...138 days ago
Python script to find repeats in the DNA sequence !
...if sequence.count(substring) > 1 and substring not in repeats: repeats.append(substring) return repeats # Example usage genome_sequence = "ATCGATCGA...138 days ago
Raku script to find microsatellites in DNA fragments !
...unt) { @microsatellites.push($substring); } } } return @microsatellites.unique; } # Example usage my $genome-sequence = "ATCGA...138 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; }...138 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, '138 days ago
Perl script to parse VCF file !
#!/usr/bin/perl use strict; use warnings; # Usage: ./parse_vcf.pl input.vcf die "Usage: ./parse_vcf.pl input.vcf\n" unless @ARGV; my $vcf_file = shift @ARGV; my @vcf_entries = parse_vcf($v...138 days ago
Perl script to parse blast results and plot basic stats !
#!/usr/bin/perl use strict; use warnings; use List::Util qw(sum); # Usage: ./parse_blast.pl blast_result.txt die "Usage: ./parse_blast.pl blast_result.txt\n" unless @ARGV; my $blast_file = shif...138 days ago