589 days ago
Perl script to parse blast results and plot basic stats !
#!/usr/bin/perl use strict; use warnings; use List::Util qw(sum); # Usage: ./parse_blast.pl blast_result.txt die "Usage: ./p...last_file = shift @ARGV; my @blast_entries = parse_blast($blast_...dentity = calculate_average(\@blast_entries, 'PERCENT_IDENTITY');...92 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_blast(query_file, subject_file, output_file): # Set up the BLAST command with format 6 (tab-de...dentity=100) # Run BLAST stdout, stderr = blastn_...15 hours ago