This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem With Rfam Search

I am trying to do a Rfam search for a unique sequence tags for small RNA. I am getting the following error message:

Command line: rfam_scan-1.0.2.pl -o rfamunmapped.txt-blastdb Rfam.fasta PolyAreomevedS1.fastq_unmapped_reads_collapsed.txt
WARN: No BLAST database specified.  Your search may be extremely slow.
Also, the model-specific global/local mode is not used, therefore
curated Rfam thresholds may not be meaningful.
Use of uninitialized value $fh in <HANDLE> at rfam_scan-1.0.2.pl line 527.
readline() on unopened filehandle at rfam_scan-1.0.2.pl line 527.
Use of uninitialized value $rfamid in hash element at rfam_scan-1.0.2.pl line 276.

Error: Failed to open covariance model save file Rfam.fasta

Can anyone help me with Rfam search? Thanks,

database small
• 5,052 views • link
updated by Paul Gardner Paul Gardner • written by Ajw • Ajw

3 answers

If what you ran on the command line is exactly as you typed it here, the problem could be a missing space between the -o switch and the -blastdb switch:

You typed:

rfam_scan-1.0.2.pl -o rfamunmapped.txt-blastdb Rfam.fasta...

Should be:

rfam_scan-1.0.2.pl -o rfamunmapped.txt -blastdb Rfam.fasta...

That might explain why the program cannot find a BLAST database, which I assume is Rfam.fasta.

• 0 views • link
updated by Ram Ram • written by Neilfws Neilfws

Have you tried indexing Rfam.fasta for blast searching? eg. formatdb.

Also, you are missing the latest Rfam.cm file.

Then try:

rfam_scan-1.0.2.pl -o rfamunmapped.txt -blastdb Rfam.fasta Rfam.cm PolyAreomevedS1.fastq_unmapped_reads_collapsed.txt
• 0 views • link
updated by Ram Ram • written by Paul Gardner Paul Gardner

Thanks Paul - yes, it looks like it is running now. AJW

• 0 views • link
written by Ajw • Ajw

Thanks neilws. I corrected the error and still gives the same results. The other thing that I am not sure is that whether or not Rfam.fasta is the correct file.

• 0 views • link
written by Ajw • Ajw

Please comment under answers, rather than writing a new answer.

• 0 views • link
written by Neilfws Neilfws

Sure, I will do this next. Thanks again

• 0 views • link
written by Ajw • Ajw

Log in to answer this question.

More posts like this
  • RSEM calculate expression troubleshoot
    written by prabin.dm prabin.dm
    Hi, I am trying to use RSEM following this [\[tutorial\]][1] This is the error code I am getting. readline() on closed filehandle INPUT at /share/pkg/RSEM/1.3.0/bin/rsem-calculate-expression …
  • Error In Genome Music Path-Scan
    written by Andreas Weller • Andreas Weller
    <p>I'm using Genome Music 0.4.1 on an Ubuntu 10.04 server. Most functions work fine, so I guess my installation and my input files are generally …
  • Parsing A Multi-Record Genbank File
    written by robjohn7000 robjohn7000
    <p>Hi,</p> <p>I tried parsing a multi-record genbank file (from this site: <a href='http://biopython.org/DIST/docs/tutorial/examples/ls_orchid.gbk'>http://biopython.org/DIST/docs/tutorial/examples/ls_orchid.gbk</a>) using the code below. </p> <p>The code returned an error: </p> <pre><code>readline() …
  • Format Of Pathway File For Genome Music: Error Upon Running
    written by DoubleD DoubleD
    <p>Hello,</p> <p>I have created a pathway file from the KEGG database and believe that I have it formatted correctly (per the geno music help files). …
  • Error When Running The Music Bmr
    written by yliubu • yliubu
    <p>I am now using MuSiC to calculate BMR, but there is an error. The command I use is </p> <pre><code>&gt; genome music bmr calc-wig-covg --wig-list=/home/share/input/wig_list.txt …
  • Error When Running The Varscan Mergesegments.Pl Tool
    written by sousuffer • sousuffer
    <p>I generated the correct columns and the mergeSegments.pl still fails:</p> <pre><code>perl ../../apps/mergeSegments.pl outp.file Use of uninitialized value $input in at ../../apps/mergeSegments.pl line 446. readline() on …
  • Merging Segments From Varscan2 Output...
    written by sousuffer • sousuffer
    <p>I am a little new to this so any help would be appreciated - I ran a variation of the standard example workflow for VarScan2 …
  • Questions About Mergesegment.Pl
    written by Yu Fu • Yu Fu
    <p>I am using VarScan to analyze somatic copy number alteration. In the recommended workflow found here:</p> <p><a href="http://varscan.sourceforge.net/copy-number-calling.html">http://varscan.sourceforge.net/copy-number-calling.html</a></p> <p>, I am stuck at the fifth …
  • Velvet Postprocessing
    written by khikho khikho
    <p>when I run Postprocessor script (version 1.6 ) on my contigs (make them by velvet) I got these errors with tones of repeats , could …
  • Find The Match Entries In The Columns And Print Them Out By Perl
    written by Ar Es • Ar Es
    <p>Hi, I have asked this question before , but as Admin said If i want to change my question , I have to ask new …