...
 
Commits (2)
......@@ -228,7 +228,7 @@ sub RunBlast
my $cmd =
"( cd $WORKDIR; cat "
. join('.fasta ', keys %{$rh_analysis_params->{'proteomes'}}) . '.fasta > all.fasta;'
. join('.fasta ', keys %{$rh_analysis_params->{'proteomes'}}) . '.fasta | sed \'s/>/\n>/g;\' > all.fasta; dos2unix all.fasta;'
. "$FindBin::RealBin/../ext/ncbi-blast/makeblastdb -in $WORKDIR/all.fasta -parse_seqids -dbtype prot 2>$WORKDIR/all.fasta.makeudb.log)";
if ($ALIGNER eq 'diamond')
......@@ -236,7 +236,7 @@ sub RunBlast
$cmd =
"( cd $WORKDIR; cat "
. join('.fasta ', keys %{$rh_analysis_params->{'proteomes'}}) . '.fasta > all.fasta;'
. join('.fasta ', keys %{$rh_analysis_params->{'proteomes'}}) . '.fasta | sed \'s/>/\n>/g;\' > all.fasta; dos2unix all.fasta;'
. "$FindBin::RealBin/../ext/diamond makedb --in $WORKDIR/all.fasta --db $WORKDIR/all.fasta.diamond 2>$WORKDIR/all.fasta.diamond.log)";
}
......