Download genome using NCBI esearch command line !
#Download using command line "NC_005956.1" / "CP003784.1" esearch -db nucleotide -query "NC_005956.1" | efetch -format fasta > seq1.fa esearch -db nucleotide -quer...1114 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...1114 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 = ){...2930 days ago
Extract a random sequence from a file
#!/usr/local/bin/perl -w use strict; use warnings; use autodie; use List::Util qw/ shuffle /; my $outputfile = 'randomoutput.txt'; open my $in_fh, '', $o...2928 days ago
2927 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-...2919 days ago
Count GC Content in nucleotide sequence with Perl
#!/usr/bin/perl -w ### Usage: get_gc_content.pl ### #-----------------------------------------------------...2919 days ago
Perl script to extract fasta sequence by matching name/ids !!
#!/usr/bin/perl use strict; use warnings; use Text::Trim qw(trim); #Usage perl extractSeqbyID.pl ids.txt seq.fasta Result.fasta $ARGV[2] or die "use extractS...2888 days ago
Converting from Windows-style to UNIX-style line endings with dos2unix
Lenovo-ideapad-320-15ISK:~/Downloads/abc/bin$ sudo apt install dos2unix [sudo] password for bioqueen: Reading package lists... Done Building dependency tree...1190 days ago
2706 days ago