RFAM problem. Died at rfam_scan.pl line 281, line 1128517
Hi everyone. I am working with RFAM. It is blocked in some line with three different commands.perl rfam_scan.pl -blastdb Rfam.fasta Rfam.cm 200pb_Pamicro.fasta --nobig -o rfam_200_nobig.out Died at...
View ArticleWhat does sframe signify in BLAST output?
Hi,I am relatively new to the field of Bioinformatics and have just been introduced to BLAST. I'm trying to understand the output available from the BLAST program. As listed in the NCBI manual, I am...
View ArticleHow Could I Run A Huge Number Of Blast Calls Faster?
For the purpose of my project I need to break down a genome and run blast for each part. Because of the amount of genome file it would be at least 4000 call of blast which would take a lot of time. I'm...
View ArticleProblems In Creating Pssm Using Psiblast
Hi, Everyone.Recently I am using PSIBlast to create PSSM.Command seems like this: psiblast -subject oneseq.fasta -inmsa alignment.fasta -outascii_pssm pssm.txtThe problems is Error: NCBI C++ Exception:...
View ArticleHow To Get A Blast Record From Stdout Of Ncbiblastxcommandline
i am blasting sequences from a fasta file individually since the output will be handled differently depending on the quality of the hit.currently i have this, which works fine:blast_handle =...
View ArticleComparable E-Value Among Blast Programs
If I specify the search space via -dbsize flag among blasts of the same query to multiple databases, the e-values will be comparable among all the results because the search space will be...
View ArticleParsing Blast Output Biopython Error
Hi, I have the following code def runBLAST(self): print "Running BLAST .........." cmd=subprocess.Popen("blastp -db nr -query repeat.txt -out out.faa -evalue 0.001 -gapopen 11 -gapextend 1 -matrix...
View ArticleBlastdbcmd - Too Many Positional Arguments , The Offending Value: %F - Windows
Hey guys, I'm trying to use blastdbcmd - when I type the following on cmdblastdbcmd -db databaseBLAST -entry_batch -outfmt "%f" -out test_query.txtthe following error pops up:Error: Too many positional...
View ArticleIdentifying Viral Sequences In Next Generation Sequencing Data
I am attempting to find viral sequences in the RNA seq data of sequenced tumors. Currently I have changed bam files into fasta files and have used the refseq viral dataset from NCBI, formatted it as a...
View ArticleBlast For Cygwin?
Which version of BLAST is prepared for Cygwin?I tried to compile ncbi-blast-2.2.26+-src.tar.gz under Cygwin 1.7.13. The compilation stopped at "checking python" at first. After python was installed,...
View ArticleGATK BaseRecalibration report
hi everyone I use GATK to do snp call. But I don't understand the BaseRecalibration report as follows . And i really don't known the word 'Observations' 'Errors' stand for. Thanks in advance RaedGroup...
View ArticleCreate a custom .gtf file with a list of genes
Hi GuysI am new to the RNA Seq world and just starting out with linux. I need to create a custom gene annotation file with a list of genes I am interested in analyzing. How do I do that.
View ArticleFormatdb For Chromosome Database?
Hi,I recently downloaded the new Fugu v5 assembly database from IMCB (http://www.fugu-sg.org/downloads/downloads3.htm). This is my first time using a chromosome level database (I am used to contig and...
View ArticleBLAST: Querying a single sequence without input file in C
This is with reference to the query posted in https://www.biostars.org/p/17265/. If I try to run the following using command line it runs perfectly fine: blastn -query <(echo -e...
View ArticleWhat software to use to map and observe substitutional patterns across...
I have protein sequences from about 50 species which I want to align and observe/retrieve key (functionally interesting) sites. Apart from using Blast (blastp) to align sequences and view query...
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 ArticleLooking For A Command Line Version Of Ncbi Taxonomy 'Stratification' Tool
Hi allI am building a pipeline for taxonomically identifying a blast result at each taxa level, with a bespoke reference dataset. I want it to be easily reproduceable, which it is, other than an ugly...
View ArticleNcbi Blast+ Taxid And Taxid_Map
2 related questions. How does the -taxid argument for the NCBI's blast-2.2.24+ makeblastdb command work? When I execute:makeblastdb -title "BO" -dbtype nucl -out BO -in BO.fna -parse_seqids -taxid 3712...
View ArticleBlast Results Wont Return
This isnt my usual way if doing things but because im using Plone it has to be done this way. The function below takes a variable "x" which is the query sequence. Its is then formatted to a blast-able...
View Article