Results for "PacBio read correction"
Bio-Scripts
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...1181 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...1147 days ago
- #!/usr/local/bin/perl -w # Can be easily modified to run any command on every seque...t"; # Go to sequence directory and open it (i.e, read contents) chdir($myDir) || d...: $!"; # Open $myDir foreach $seqFile (sort readdir(DI...
2961 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...2962 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...2951 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,"2920 days ago
- P
2917 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'){...2911 days ago
- #!/usr/local/bin/perl -w # NOTE: YOU MUST CHANGE THE LINE ABOVE TO POINT TO # THE FU...r sequence. This database is not human-readable, and can get...j, $k, $key) = (0, 0, 0, ""); $seq = ""; #Read the s...
2738 days ago
Transpose the file coordinates and plot dendrogram in R
#Save this as tr.awk { for (i=1; i2662 days ago
