51 days ago
Perl script for chi-squared test !
...file1_ob{$di}\t$file1_ex{$di}\t$chi1\t$file2_ob{$di}\t$file2_ex{$di}\t$chi2\n"; $chi_total += ($chi1+$chi2); } printf "Chi squared value = %6.2f\n", $chi_total; pri...470 days ago
358 days ago
Raku script to find microsatellites in DNA fragments !
sub find-microsatellites($sequence, $min-repeat-length = 2, $max-repeat-length = 6, $min-repeat-count = 3) { my @microsatellites; for my $repeat-length ($min-repeat-len...153 days ago
Python script to finds extact similar sequence between two multi fasta files !
..., subject_file, output_file): # Set up the BLAST command with format 6 (tab-delimited) blastn_c...mmandline(query=query_file, subject=subject_file, out=output_file, outfmt=6,...62 days ago
Python script to parse GFF file
...'start': int(fields[3]), 'end': int(fields[4]), 'score': fields[5], 'strand': fields[6], 'phase...97 days ago