131 days ago
Perl and BioPerl script to extract protein sequences using GFF file !
#!/usr/bin/perl use strict; use warnings; use Bio::DB::Fasta; use Bio::SeqIO; # Paths to your GFF file and genome FASTA file my $gff_file = 'path/to/your/file.gff'...131 days ago
Raku script to find repeats in sequences !
sub find-repeats($sequence, $min-repeat-length = 3) { my @repeats; for ^($sequence.chars - $min-repeat-length + 1) -> $i { my $substring = $sequence.s...131 days ago
Raku script to find microsatellites in DNA fragments !
sub find-microsatellites($sequence, $min-repeat-length = 2, $max-repeat-length = 6, $min-repeat-count = 3) { my @microsatellites; for my $repeat-length ($min-r...131 days ago
Perl script to calculate the basic stats of the assembled genome !
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; # Input file containing the genome assembly in FASTA format my $input_file = 'genome_assembly.fasta'; #...131 days ago
Python script for basic stats of the assembled genome !
from Bio import SeqIO import statistics # Input file containing the genome assembly in FASTA format input_file = 'genome_assembly.fasta' # Variables for computing...131 days ago
Python script to finds extact similar sequence between two multi fasta files !
from Bio.Blast.Applications import NcbiblastnCommandline import os import sys def perform_local_blast(query_file, subject_file, output_file): # Set up the BLAST...40 days ago
Python script to calculate basic genome stats !
from Bio import SeqIO def calculate_genome_stats(fasta_file): # Initialize variables to store genome statistics genome_length = 0 gc_count = 0 a_cou...22 hours ago