657 days ago
Perl script to read the next line of a file !
my $line = ; while(1) { # keep looping until I say so my $nextLine = ; if ($line =~ m/>/ || !defined $nextLine) { ### Do the stuff } ### Do any other stuff; last unless defined $nextLine; $line = $nextLine; }558 days ago
Perl script to find inverted repeats !
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; use Bio::Tools::Run::RepeatMasker; my $genome_file = "genome.fasta"; # read genome sequence my $seqio = Bio::...417 days ago
Perl script for chi-squared test !
...############################################################ # Read sequence file 1 ############...E); ############################################################ # Read sequence file 2 ############...403 days ago
291 days ago
86 days ago
Perl script for six frame translation !
...ings; use Bio::SeqIO; # Path to your input nucleotide sequence file in FASTA format my $input_fasta = 'path/to/your/input.fasta'; # Step 1: Read the input FASTA file my $seq...86 days ago
Raku script to find overlaps between two bed files !
...rrect number of arguments are provided if @*ARGS.elems != 2 { say "Usage: ./compare_bed_files.raku file1.bed file2.bed"; exit 1; } # Read the contents of the two BED f...86 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, '86 days ago
Perl script to calculate the basic stats of the assembled genome !
...::SeqIO; # Input file containing the genome assembly in FASTA format my $input_file = 'genome_assembly.fasta'; # Create Bio::SeqIO object to read the FASTA file my $seqio = B...86 days ago