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() else:...ord in records: print(f"Sequence ID: {record.id}") fo...s(frame)} Translation:\n{protein_sequence}...122 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.substr($i, $min-repeat-length...peats; } # Example usage my $genome-sequence...122 days ago
Python script to find repeats in the DNA sequence !
def find_repeats(sequence, min_repeat_length=3): repeats = [] for i in range(len(sequence) - min_repeat_length + 1): substring = sequence[i:i+min_repeat_length] if s...return repeats # Example usage genome_sequence...122 days ago
Raku script to find microsatellites in DNA fragments !
sub find-microsatellites($sequence, $min-repeat-length = 2, $max-repeat-len...max-repeat-length) { for ^($sequence.chars - $repeat-length * $min...> $i { my $substring = $sequence.substr($i, $repeat-length);...esult = find-microsatellites($genome-sequence)...122 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...gth = 0; my $num_contigs = 0; my @contig_lengths; # Iterate through each sequence...122 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...al_length = 0 num_contigs = 0 contig_lengths = [] # Iterate through each sequence...122 days ago
Python script to convert fastq to fasta
def fastq_to_fasta(fastq_file, fasta_file): with open(fastq_file, 'r') as fq:...header = fq.readline().strip() sequence = fq.readline().strip()...>' + header[1:] + '\n') fa.write(sequence...66 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_bl...'\t') qseq = fields[0] # Extract the aligned query sequence (qseq) #print("Aligned Query Sequence:...31 days ago