135 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.gf...135 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()...135 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...135 days ago
Python script to create fastq file with random sequences
import random def generate_random_sequence(length): bases = ['A', 'C', 'G', 'T'] return ''.join(random.choice(bases) for _ in range(length)) def generate...5 days ago
Python script to parse a FASTQ file !
...parse_fastq(fastq_file): # Initialize a list to store parsed sequences sequences = [] # Read the sequen...": record.letter_annotations["phred_quality"] } sequences.append(sequence_info)...5 days ago