2870 days ago
2870 days ago
2870 days ago
Perl subroutine to read and write files
...ile.txt',\$string); sub InOut { my($bit,$file,$data) = @_; if($bit eq 'read'){ open InOut,"< $file" or die "Cannot open $file for input: $!\n"; my @file = ; close InOut; return wantar...2864 days ago
Perl script introduces control structures, arrays and hashes.
...!/usr/bin/env perl use strict; use warnings; my @first_array = ('DNA', 'ATGCGTGC', 5, 'RNA', 'AUGC'); pri...first_array[0], "\n\n"; # Scalar gives actual size of an a..._array = @first_array; # implicit way print "Scalar of array...'RNA' => 'AUGC', 'Number of...2861 days ago
Blast script to index and extract sequence !!
# look at the file $ head EC4115.fa >NC_011353.1 Escherichia coli O157:H7 str. EC4115 chromosome, complete geno...ATATGTCTCTGTGTGGATTAAAAAAAGAGTCTC TGACAGCAGCTTCTGAACTGGTTACC...base $ makeblastdb -dbtype nucl -out EC -in EC4115.fa -parse...TAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGAC...2851 days ago
Perl script to insert sequence in contig !!
# sub signature: #insertSEQintoCONTIGatLOC( SEQ , CONTIG , LOC ) ; sub insertSEQintoCONTIGatLOC{ my ( $SEQ , $CONTIG , $...@_; substr( $CONTIG , $LOC , -length($CONTIG) ) = $SEQ ; return $CO...rt = '.($s = insertSEQintoCONTIGatLOC( "CCCC ", $s , 26 ))."\n"; p...2707 days ago
2703 days ago
2692 days ago
2692 days ago