Blast result parser with Perl and Bioperl
...sion number, length, description, E value, bit score, query frame, query start, qu...sion_number\tlength\tdescription\tE value\tbit score\tframe\tquery_start\t"; prin.... "\t"; #get the bit score of the hit print OUT $hit...2924 days ago
Needleman-Wunsch Algorithm in Perl
...matrix and gapcost from command line my ($seq1, $seq2, $smfile, $gapcost) = @ARGV; # scoring scheme (instead of using fixed MATCH and MISMATCH scores we will use values read from...2915 days ago
Create genome scaffolding with Perl
...y $alSeqL = ""; for (my $i = 0; $i {"prefix"}, $nextScaffoldID++); if(!exists($replacementSeqs{$sName}{score}) || ($trimTotal < $replacementSeqs{$sName}{score}) || (($trimTotal =...2318 days ago
Plot the clock using Lastz -gerenal outfile
...h) { #next if $pHash{$val} > 10; my @sChr = split '\_', $val; my $score=$pHash{$val}/$sChr[1]; push @allKeys, $sChr[0]; push @allFreq, $score; } my $allKeys=join ',',...2309 days ago
Estimate Genome Size with Jellyfish and R
...ation #-o -output=string Output file name #--min-quality-char Base quality value. Version 2.2.3 of Jellyfish uses the “Phred” score, where "?" = 30 jellyfish...2256 days ago
2137 days ago
2109 days ago
2085 days ago
Perl script for Smith-Waterman Algorithm
...that match my $MISMATCH = -1; # -1 for letters that mismatch my $GAP = -1; # -1 for any gap # initialization my @matrix; $matrix[0][0]{score} = 0; $matrix[0][0]{pointe...960 days ago
Perl script to parse blast results and plot basic stats !
...lue: ", $entry->{EVALUE}, "\n"; print "Bit Score: ", $entry->{BITSCORE}, "\n"; print "----...e\n"; print "Bit Score Range: $min_bitscore - $max_bitscore\n"; print "Average Bit Score...104 days ago