User Tools

Site Tools


blast_and_plast

This is an old revision of the document!


Here is a shell to run plast:

#!/bin/bash
#$ -S /bin/bash
. /etc/profile
#$ -o logo
#$ -cwd
#$ -pe threaded 10
cd $PWD
CPUs=10
DB=database
QF=yourquery.fasta
plast -e 1e-10 -max-hit-per-query 1 -outfmt 1 -a $CPUs -p plastn -i $QF -d $DB -o $QF.plout

to parse the output see http://129.173.88.134:81/dokuwiki/doku.php?id=dayana_salas_-_utility_scripts_taxonomy_coloring_trees_phylogenetics_mixture_models_domain_architecture_and_more

Here is a shell example to run blast using MMETSP database:

#!/bin/bash
#$ -S /bin/bash
. /etc/profile
#$ -cwd
source activate blast DB=/db1/extra-data-sets/MMETSP/MMETSP_db/MMETSP_DB_clean.v2018.fa
blastn -db $DB -query yourquery.fasta -out yourqueryresults.blout
source deactivate

blast_and_plast.1530726204.txt.gz · Last modified: by 24.138.71.142