Perl script to calculate the basic stats of the assembled genome !
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; # Input file contain...0; my @contig_lengths; # Iterate through each sequence in the assembly while (my $s...t($input_file); # Print the computed statistics and infor...108 days ago
Python script for basic stats of the assembled genome !
from Bio import SeqIO import statistics # Input file containing the genome assem...= 0 contig_lengths = [] # Iterate through each sequence in the assembly for record i..._percentage, 2) # Print the computed statistics and infor...108 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...52 days ago
Python script to finds extact similar sequence between two multi fasta files !
from Bio.Blast.Applications import NcbiblastnCommandline import os import sys...(output_file, "r") as result_handle: for line in resul...# Extract the aligned query sequence (qseq) #print("Aligned Query Sequence:", qseq) # Print...17 days ago