wget ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz && gunzip nr.gz &>log.nrdb &
wget ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/swissprot.gz && gunzip swissprot.gz &> log.swissprotdb &
wget ftp://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.fasta.gz &>log.uniprotdb &
wget ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_trembl.fasta.gz &>log.trembldb &
curl -o plant.#1.protein.faa.gz ftp://ftp.ncbi.nlm.nih.gov/refseq/release/plant/plant.\[1-198\].protein.faa.gz
#Create uniprot_database makeblastdb -in uniprot_sprot.fasta -dbtype prot &>log &
#Run database
blastx -num_threads 64
-query alatum.transcripts.fasta
-db /NABIC/HOME/senthil/002_DB/uniprot_sprot.fasta
-out uniprot.xml -outfmt 5 &>log &
#Create nr_database
makeblastdb -in nr -dbtype prot &>log &
#Run database
blastx -num_threads 64
-query alatum.transcripts.fasta
-db /NABIC/HOME/senthil/002_DB/nr
-out nr.xml -outfmt 5 &>log &