Get Blast Database Size
I have piped unknown length sequences into makeblastdb. Now I want to know total length of those sequences (BLAST database size).Example: # "cat" is used as an example # My "real" sequences are piped...
View ArticleError Using Makeblastdb
Hi there,I am trying to run makebastdb to create a new BLAST database. The program starts but cuts out abruptly with an Error:Error: ncbi::s_FixBioseqDeltas() - Bioseq must have Seq-data or Delta...
View ArticleBlast Output With Species Name
Dear all,my question is very simple. I have reads from a human RNASeq and I would like to check for contamination on a specific virus family. So my pipeline looks like1) Align on the human genome...
View ArticleGenerating Protein Consensus Sequence For Blast Input
I have a set of 3 protein sequences for which I need to find the "consensus sequence" which can be used as an input to blast against a local database.how should I go about it?
View ArticleFinding Protein Homology?
I'm looking for something that's more rigorous than BLAST and less rigorous than a global alignment-esque algorithm. Is there a paper out there that compare methods? What's the most popular method?I'm...
View ArticleReciprocal blast and protein identification
Hi Guys,I am doing reciprocal BLAST to identify members of a protein family in different species. Is it necessary for me to set a strict e-value threshold while doing each BLAST? What I am doing is I...
View ArticleWindows Cmd Not Responding With Python Blastall 2.2.17
I'm using the standalone BLAST 2.2.17 with python. However, I'm having a problem with the cmd. When running the script the cmd stops responding at the 15th blast record. It doesn't matter if I change...
View ArticleHow To Blast Against A Local Database
I want to BLAST an mRNA reference sequence against a local database in my directory (a collection of RNA reference sequences) using a BioPython program. I know how to perform BLAST against online...
View ArticleBlastn And Blastcmd Issues - Blast
Hey all, first of all I'd just to say that I am fairly new to the field of bioinformatics and I'm, unfortunately, learning by my mistakes more than my successes atm and I'm having some trouble...
View ArticleError in run blastpgp
Hi, I am trying to have the command line blast work, however, I did not manage to execute it properly. I would appreciate if you can tell me what I am doing wrong. So here are the things I did so far:...
View ArticleHow To Use To Use Standalone Blast In Biopython
I've been trying to get Biopython to use BLAST (the standalone version) and it doesn't seem to work. I'm using blast 2.2.28, python 2.6 and biopython 1.61. I was wondering if anyone could give a short...
View ArticleAligned Ranges In Blast
I have several sequences which I'm blasting against the nt database, and I wanted to know how to retrieve for each hit the intervals matched. As far as I'm aware, the standard output only gives this...
View ArticleDifference Between Blastall And Blastn?
What is the difference, other than subtle different in the output format, between "blastall -p blastn" and running NCBI-Blast+ blastn?Are the results supposed to be the same in terms of...
View ArticleHow To Create A Pssm From Fasta Homologues With Ncbi Blast+ 2.2.23
I have a FASTA sequence file with about 10 homologous proteins. What I would like to do is create a PSSM from them and use it to search a transcriptome database.But how to create it? There is a makemat...
View ArticleWarning Unable To Open .Nin
I am using the latest version of blast from ftp://ftp.ncbi.nlm.nih.gov/blast/executables/release/LATEST/ for linux. I created the database for Bos_taurus using Formatdb commamd of...
View ArticleMaking Specialised A Local Blast Db
Hi,I need to make a local Blast DNA (nt) database that only includes sequences from plants. I'd also like to make one with only chloroplast sequences. Put I can't find any tools to do this. I don't...
View ArticleImport Xml Blastx Results In Geneious
Hi,I have the xml results file of my local blastx and I would like to import the annotations directly on my sequence in geneious. Do you know if there is a way ?Thanks,Aurelien
View ArticleHow to create a Blast database of viruses ?
Hi all,For a metagenomic project a want to make a blast database of viruses. I dont want to blast my reads on the entire nt database. But I dont know how to make it. Downloading the result of a query...
View ArticleIs Blast+ Running As Fast As It Could ?
Hi all, I am trying to run a tblastn search of nt locally using blast+2.2.25, and was wondering if anyone new if it's going as fast as it should - I have 1000's of queries to get through! I should...
View ArticleGet Snps, Indels From Biopython Blast Parser
Hi all,Is there a way through the bioPython API to get the list of variants from BLAST alignments?I have done my alignments using standalone BLAST by providing -query and -subject, and the output is...
View Article