cgeb2001's DokuWiki!

You've loaded an old revision of the document! If you save it, you will create a new version with this data.

Here is a shell to run **PLAST**:\\

<code>
#!/bin/bash
#$ -S /bin/bash
. /etc/profile
#$ -o logo
#$ -cwd
#$ -pe threaded 10
cd $PWD
CPUs=10
DB=/db1/blast-may-2024/nt.nal
QF=yourquery.fasta
plast -e 1e-10 -max-hit-per-query 1 -outfmt 1 -a $CPUs -p plastn -max-database-size 10000000000 -i $QF -d $DB -o $QF.plout -force-query-order 1000
</code>

to parse the output see https://perun.biochem.dal.ca/user-wiki/doku.php?id=taxonomy_recovery

Here is a shell example to run **BLAST**:\\

<code>
#!/bin/bash
#$ -S /bin/bash
. /etc/profile
#$ -pe threaded 1
#$ -cwd
source activate blast
export BLASTDB=/db1/blast-may-2024/
DB=nt
query=your_query.fasta
blastn -db $DB -query $query -out yourqueryresults.blout -num_threads 1 -outfmt "6 qseqid sseqid evalue pident qcovs length slen qlen qstart qend sstart send stitle"
source deactivate
</code>

Both shells using NCBI nt database (/db1/blast-may-2024/nt.nal), but the formats for specifying DB are different for BLAST and PLAST.

Please fill all the letters into the box to prove you're human. E᠎ V P Y H Please keep this field empty:

Edit summary

Note: By editing this page you agree to license your content under the following license: CC Attribution-Noncommercial-Share Alike 4.0 International