2924 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...',$file,\@array) #$string = "YO!"; #InOut('write','file.txt',\$string...: $!\n"; my @file = ; close InOut; return wantarray ? @f...2918 days ago
Perl script introduces control structures, arrays and hashes.
...ctual size of an array my $size_of_array = scalar(@first_array); my $another_way_of_getting_size_of_ar...ay\n\n"; print "Perl's index size of array: $#first_array\n\n"; print "Another way of getting size: $ano...2914 days ago
Blast script to index and extract sequence !!
...coli O157:H7 str. EC4115 chromosome, complete genome. AGCTTTTCATTCTGACTGCAACGGG...coli O157:H7 str. EC4115 chromosome, complete genome AGCTTTTCATTCTGACTGCAACGGGC...coli O157:H7 str. EC4115 chromosome, complete genome TTAAAA...2904 days ago
Perl script to insert sequence in contig !!
...'.($s = insertSEQintoCONTIGatLOC( "CCCC ", $s , 26 ))."\n"; print " (now, s = '$s' ) \n"; # OUTPU...TATGATGATAGATGATAGTAGATAGATAGATAGATAGATAG # insert = 'CCCC ' # (now, s = ATATGATGATAGATGATAGTAG...2760 days ago
2756 days ago
2745 days ago
2745 days ago
Perl script to insert the DNA string in genome
#!/usr/bin/perl use warnings; use strict; use Bio::SeqIO; use Bio::Seq; my $file = $ARGV[0]; # input fasta file (genome file) my $out = $ARGV[1]; # output fasta file my...2711 days ago
BASH script for SelfBLAST a genome
#!/bin/bash #self BLAST a genome -- Expecting you have blast and samtools installed in your system #Au...[ -f $MYDB.nhr ] then echo "BLAST database for MergedContigs.fasta genome exists" else echo "Than...2704 days ago