Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.org/ has...dna = Entrez.efetch(db='nucleotide',id=nm, rettype = 'fasta', retmode= 'text').read().spl...1274 days ago
R script to visualize fastANI core-genome comparison
####### # Purpose: Visualize fastANI core-genome comparison # Usage: Rscript # Output: .pdf # Uses genoP...s(TRUE)[2] fastANI_visual_file=commandArgs(TRUE)[3] library(genoPlotR) #Read fastA...1240 days ago
3054 days ago
Read lines from input file – print lines that match a regular expression
#!/usr/bin/perl use strict; use warnings; my $line; #read the line-by-line for each line ask if it matches the regex print it if it matches while($line = ){ chomp $l...3055 days ago
Needleman-Wunsch Algorithm in Perl
#!/usr/bin/perl # USAGE: perl nw.pl HEAGAWGHEE PAWHEAE BLOSUM50.txt -8 # See:...fixed MATCH and MISMATCH scores we will use values read from BLOSUM50) my $MATCH...any gap my %BLOSUM50 = (); my @aalist = (); # read subst...3044 days ago
Perl script to extract lines with matching ids !!
#!/usr/bin/perl use strict; use warnings; my %patterns; #USAGE: perl extactByIds.pl Idsfile1 file2 > Result # Open file and get patterns to search for open(my $fh2,"3013 days ago
3010 days ago
Perl subroutine to read and write files
# Input output (InOut) the file # usage: # @array = InOut('read',$file) # $string = InOut('read',$file) # InOut('write',$file,\$string) # InOut('write',$fi...g); sub InOut { my($bit,$file,$data) = @_; if($bit eq 'read'){...3004 days ago
2831 days ago
Transpose the file coordinates and plot dendrogram in R
#Save this as tr.awk { for (i=1; i2755 days ago