- or download this
sub getSeqfromfasta2lmers
{
...
} #end while
return @seqs;
}
- or download this
my $nofseq = scalar( getSeqfromfasta2lmers( $file ) );
- or download this
my @input_seqs = getSeqfromfasta2lmers($file);
- or download this
my $ip = @input_seqs;
- or download this
sub getSeqfromfasta2lmers {
my $file = shift;
...
my $in = Bio::SeqIO->newFH( -format => 'fasta', -fh => $fh );
return <$in>; ### MUST BE CALLED IN A LIST CONTEXT!!!!
}
- or download this
## Get the sequences.
my @input_seqs = getSeqfromfasta2lmers($file);
## And a count of them
my $nofseq = @input_seqs;
- or download this
sub main_process {
....
...
## replace references to $ip with $nofseq
## DELETE my $ip = @input_seqs;
- or download this
sub getlmersfromseq {
my ($seqsarr,$l)= @_;
...
my @uniq_lmers = uniq @lmers;
return @uniq_lmers;
}
- or download this
sub getlmersfromseq {
my ($seqsarr,$l)= @_;
return uniq map {substrings $_, $l} @{$seqsarr};
}
- or download this
sub gen_param {
my ( $file, $file_neg, $nofseq ) = @_;
...
}
return @param_groups;
}
- or download this
$hash{ 'some complicated key' } = [@matches];
- or download this
$hash{ 'some complicated key' } = \@matches;
- or download this
my $output = main_process( @{$_}{
qw/ file file_neg submt_len submt_d e W_size lp support_thresh
+old min_inst_lower min_inst_upper polyTA_lim poly_lim /
});
$result{ 'ParamGroup' . $count++ } = $output;
- or download this
} # ----- end foreach $mcands -----
return;
}