Results for "line"
Bio-Scripts
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 -query "CP003784...1129 days ago
R script to visualize fastANI core-genome comparison
...I core-genome comparison # Usage: Rscript # Output: .pdf # Uses genoPlotR package: http://genoplotr.r-forge.r-project.org #Parse command line arguments query_fasta=comman...1129 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 $line;...2944 days ago
Extract a random sequence from a file
...my $outputfile = 'randomoutput.txt'; open my $in_fh, '', $outputfile; my $size = 21; my $count = 10; while (my $line = ) { next unless $line =~ /^([ATGCN]+)/; my $g...2942 days ago
- #!/usr/bin/env perl # Usage: fastaRead.pl data.fa use strict; use warnings; my $filename = $ARGV[0]; my $sequence; open my $fileH, "
2941 days ago
Needleman-Wunsch Algorithm in Perl
...alignment algo (GOTHO 1982 mod) # usage statement die "usage: $0 \n" unless @ARGV == 4; # get sequences, matrix and gapcost from command line my ($seq1, $seq2, $smfile, $...2933 days ago
Count GC Content in nucleotide sequence with Perl
...seq = ""; while () { chomp; if (/^>/) { #finish up previous line. if (length($seq) > 0) { &process_it; } #start new line. $id = $_; $id =~ s/^>(.+...2933 days ago
Perl script to extract fasta sequence by matching name/ids !!
...; my $fasta = shift @ARGV; my $out = shift @ARGV; my %select; open LINE, "$list" or die; while () {...0]} = 1; } my $size = keys %select; print "Total Ids $size\n"; close LINE; $/ = "\n>"; open OUT, ">...2902 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...1204 days ago
- .../local/bin/perl -w # NOTE: YOU MUST CHANGE THE LINE ABOVE TO POINT TO # THE FULL...Databases size scales linearly with sequen...%hash = (); } dbmclose %TEMPFILE; $img->line($x0, $y1, $x1, $y1, $gray);...
2720 days ago
