Perl script to find coding regions in DNA sequences
#!/usr/bin/perl -w use strict; # if the number of input arguments is lower than 2 # return a message showing the error if (scalar(@ARGV) < 2) { print "dnalogl...2209 days ago
Biological Sequence handling with Perl !
package Sequence::Generic; # File: Sequence/Generic.pm use strict; use Carp; use overload '""' => 'asString', 'neg' => 'reverse', '.' =>...2235 days ago
Perl subroutine for reading multifasta file !
sub readfasta { (my $file)=@_; my %sequence; my $header; my $temp_seq; #suppose fasta files contains multiple sequences; open (IN, "2198 days ago
Perl script to run SATSUMA in loop !
#!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage perl 1by1.pl for SATSUMA analysis # User need to set the reference multifasta file name here my...2194 days ago
Perl script to convert GFF 2 FASTA !
#!/usr/bin/perl use strict; use warnings; use Bio::Seq; use Bio::SeqIO; use Bio::DB::Fasta; $| = 1; # Flush output my $outfile_cds = Bio::SeqIO->new( -format =>...2193 days ago
2180 days ago
2177 days ago
Perl script to extract sequence by Ids from multifasta file !
#!/usr/bin/perl -w use strict; my $idsfile = "$ARGV[0]"; my $seqfile = "$ARGV[1]"; my %ids = (); open FILE, $idsfile; while() { chomp; $ids{$_} += 1; }...2150 days ago
2129 days ago
Perl script to extract a sequence from multifasta with range !
# filterfastarange.pl #!/usr/bin/perl use strict; use warnings; #perl filterfastarange.pl 301 600 contigs.fasta > contigs-gt300-lte600.fasta my $minlen = shift or d...2120 days ago