Raku script to find SSRs in fastq file !
sub find-ssrs(Str $sequence) { my @ssrs; for 2..$sequence.chars -> $min-repeats { for $sequence.chars...$min-repeats -> $max-repeat { my...156 days ago
138 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...138 days ago
Perl script for six frame translation !
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; # Path to your input nucleotide sequence file in FASTA format my $input_fasta = 'path/to/your/input.fasta...138 days ago
Python script for six frame translation of sequences !
from Bio import SeqIO from Bio.Seq import Seq def translate_frame(sequence, frame): if frame > 0: translated_seq = sequence[frame-1:].translate() e...138 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_co...8 days ago
Python script to parse a FASTQ file !
#Python script to parse a FASTQ file and extract basic information such as the sequence identifier, sequence, and quality scores #pip install biopython from Bio impor...8 days ago