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 sequ...:RepeatMasker->new(); my $rm_report = $rm->run($genome_...474 days ago
Download lumpy skin disease data !
Location https://www.ncbi.nlm.nih.gov/sra?linkname=bioproject_sra_all&from_uid=880745 The raw genome sequence data from the 2022 outbreak in India is available in the SRA Project PRJNA880745459 days ago
143 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'; my $genome_...143 days ago
Raku script to find repeats in sequences !
sub find-repeats($sequence, $min-repeat-length = 3) { my @repeats; for ^($s...return @repeats; } # Example usage my $genome-sequence = "ATCGATCGATCGATCG"; my @result = find-repeats($genome-...143 days ago
Raku script to find microsatellites in DNA fragments !
sub find-microsatellites($sequence, $min-repeat-length = 2, $max-repeat-length = 6, $mi...microsatellites.unique; } # Example usage my $genome-sequence = "ATCGATCGATCGATCGATCG"; my @result = find-microsatellites($genome-...143 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'; # Create Bio::SeqIO object to...nt the computed statistics and information print "Genome...143 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 statistics...# Print the computed statistics and information print("Genome...143 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...= 0 g_count = 0 # Read the genome sequence from the FASTA file...content gc_content = (gc_count / genome_...13 days ago