Blast result parser with Perl and Bioperl
...(-format => 'blast', -file => $infile); open (OUT,">$outfile") or die "Cannot open $outfile: $!"; # print the header info for tab-deliminated columns print OUT "query_name\tquery...2922 days ago
Needleman-Wunsch Algorithm in Perl
#!/usr/bin/perl # USAGE: perl nw.pl HEAGAWGHEE PAWHEAE BLOSUM50.txt -8 # See: "Biological sequence anaysis" Durbin et al. ed. CUP 1998, Pg. 19 # Needleman-...2914 days ago
2880 days ago
Plot the density of genes in R
...ges())){ library(ggplot2) } else { install.packages("ggplot2") library(ggplot2) } # import a text file with gene positions # columns should be: chr, position (no...2292 days ago
Perl script to find coding regions in DNA sequences
...%pcodons; my ($codon,$freq); my $line; # for each line in the file (tableofcodons) do while ($line = ) { # extract the two values or columns with a regular expression #...2163 days ago
2083 days ago
Samtools commands for bioinformatician !
...breadth samtools depth -a sorted_dupremoved.bam | awk '{c++; if($3>0) total+=1}END{print (total/c)*100}' ### flagstat # output file has two columns: QC-passed reads and QC-faile...1611 days ago