1205 days ago
Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.org/ has...d().split("\n")[1:] allseq[nm] = ''.join(dna) import json with open("data/allseq.j...1124 days ago
R script to visualize fastANI core-genome comparison
####### # Purpose: Visualize fastANI core-genome comparison # Usage: Rscript # Output: .pdf # Uses genoPlotR package: http://genoplotr.r-forge.r-project.org #P...1089 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 match...?A{3,10}$/) { print "$line\n"; } } exit(); __DATA__ ATGCC...2905 days ago
Parse a genbank file using regular expressions
#! /usr/local/bin/perl -w $genbank = "genbank_file.txt"; open (GENBANK, $genbank) || die "cannot open $gb_...nslation product not found.\n"; } } else { # Skip this data } }2905 days ago
A multilayer perceptron (MLP) neural network in Perl
#!/usr/local/bin/perl -w ############################################...eate an #extra input in the training data #and set to 1 #########...my $LR_HO = 0.07; ################ Data dependent settings #########...s initWeights(); # load in the data initD...2904 days ago
Perl program to implement sliding window !
#!/usr/bin/perl -w my $filename = 'data.txt'; open(my TR, '2902 days ago
2902 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,"2863 days ago
Perl subroutine to read and write files
# Input output (InOut) the file # usage: # @array = InOut('read',$file) # $string...ile.txt',\$string); sub InOut { my($bit,$file,$data) = @_; if($bit eq 'read'){...en $file for output: $!\n"; print InOut ref $data eq 'ARRA...2853 days ago