Hello As novice in bioinformatics,I have a problem with blast algorithm in order to implement in matlab. I have 32 genomes and I want to find ortholog genes using blastp. However, as I see matlab as well as NCBI blast algorithm I cannot find any solution. In fact, I do not how to run blastp, create local NCBI database, and ... . My main point is to make pairwise aligment between each two genes from different 32 genomes and find the orthologs. I would be grateful if someone gives a short example of how to run this module. Considerinf two genomes from ftp://ftp.ncbi.nlm.nih.gov/genomes/archive/old_refseq/Bacteria/Staphylococcus_aureus_MRSA252_uid57839/ and ftp://ftp.ncbi.nlm.nih.gov/genomes/archive/old_refseq/Bacteria/Staphylococcus_aureus_MSSA476_uid57841/ can anyone give me a sample of code in matlab? Thank you in advance
0 answers
No answers yet.
Log in to answer this question.
-
jannovar download problem
written by curious
-
blastp search shows Blast database error: no alias or index file found for protein database
written by garfield320 •
-
How to download all the bacteria genomes from NCBI?
written by anran04100 •
-
Get genomes in genbank that are NOT in RefSeq
written by rororo •
-
download refseq of thousand of assembly file from NCBI
written by Shelle •
-
How to add ptt file content to fasta header
written by Promi •
-
How to choose NCBI viral database?
written by Tao
-
Which database to select for bacterial ortholog and singleton detection with ncbi-blast+?
written by Naren
-
How To Generate Blast Database
written by Ongkrab •
-
Download All The Bacterial Genomes From Ncbi
written by rehma.ar •
If you need to work in Matlab then you are best off contacting their technical support.
If you are able to work on command line then you can follow NCBI's command line handbook for blast here.
Just curious, but why do you need to reimplement this in Matlab? There are lots of tools and databases available to do this.