Mercurial > repos > dereeper > pangenome_explorer

view COG/bac-genomics-scripts/rod_finder/blast_rod_finder_legacy.sh @ 3:e42d30da7a74 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author dereeper
date Thu, 30 May 2024 11:52:25 +0000
parents
children
line wrap: on
line source

#!/bin/bash
echo "###Running legacy BLASTN with subject '$1' and query '$2'"
formatdb -p F -i $1 -n ROD
blastall -p blastn -d ROD -i $2 -o blastn.out -e 2e-11 -F F
echo "###Running blast_rod_finder.pl with query '$3' and minimum ROD size '$4'"
perl blast_rod_finder.pl -q $3 -r blastn.out -m $4