Python script for six frame translation of sequences !
from Bio import SeqIO from Bio.Seq import Seq de...translated_seq def six_frame_translation(fasta_file): records = list(SeqIO.parse(fasta_file, "fasta")) for re...uence}\n") # Replace 'path/to/your/input.fasta' with the actual path to your...138 days ago
Perl script to calculate the basic stats of the assembled genome !
#!/usr/bin/perl use strict; use warnings; use Bio...nput file containing the genome assembly in FASTA format my $input_file = 'genome_assembly.fasta'; # Create Bio::SeqIO obje...eqIO->new(-file => $input_file, -format => 'fasta'); # Variables for computi...138 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...bly for record in SeqIO.parse(input_file, 'fasta'): length = len(record.s...138 days ago
92 days ago
Python script to convert fastq to fasta
def fastq_to_fasta(fastq_file, fasta_file): with open(fastq_file, 'r') as fq: with open(fasta_file, 'w') as fa:...# Write to the FASTA file fa.writ...82 days ago
Python script to finds extact similar sequence between two multi fasta files !
from Bio.Blast.Applications import NcbiblastnCommandline import os import sys...if len(sys.argv) != 4: print("Usage: python script.py query.fasta subject.fasta output.txt") sys.exi...47 days ago
Python script to calculate basic genome stats !
from Bio import SeqIO def calculate_genome_stats(fasta_file): # Initialize variables to stor...# Read the genome sequence from the FASTA file for record in SeqIO...t: {gc_content:.2f}%") # Example usage fasta_file = "path/to/your/genome.f...8 days ago