Tutorial: Troubleshooting issues with SURPI pipeline for pathogen identification from complex metagenomic NGS data.
8
gravatar for geek_y
5.1 years ago by
geek_y10k
Barcelona
geek_y10k wrote:

SURPI is a pipeline to find out the pathogens from clinical metagenomics samples. It is tested only on Ubuntu and it assumes many things about your installation. I recently installed in on CentOS. These are few key points you need to take care before running the pipeline.

The github page of SURPI https://github.com/chiulab/surpi and published in Genome Research.

  1. The  create_snap_to_nt.sh program uses -Ofactor as 1000, on line 29, which may not work for your machine. You need to figure out the correct value and make necessary changes. Read snap aligner document for details.
  2. The abyss instalation requires mmap. Make sure you have installed it before compiling abyss.  http://hackage.haskell.org/package/mmap-0.5.9/mmap-0.5.9.tar.gz
  3. Make sure formatdb is there in your path. It can be downloaded from ftp://ftp.ncbi.nlm.nih.gov/blast/executables/release/LATEST/
  4. The taxonomy_lookup.pl program, at line 84 has sort --parallel=$cores, where you may need to remove --parallel=$cores option, if the sort utility on you machine does not support --parallel option.
  5. The abyss_minimus.sh program tries to use mpirun to make it parallel. If the mpirun is not configured properly, you need to remove the option 'np=$cores' in line 86, so that it will not be run parallelly. 
  6. The ribo_snap_bac_euk.sh program is hardcoded to use the 10,75 as arguments to crop_reads.csh, which you may need to change in line 43.
  7. The coveragePlot.py program uses mlab.load() at line 47, which is deprecated in latest version of matplotlib. Hence, you may need to change it to np.loadtxt()

I will update this post as and when I find more issues with the pipeline.

I forgot to mention that the configuration file that gets created, will keep the wrong path for the reference sequences. For example, in your case the path to db is dbname=/reference/taxonomy/gi_taxid_nucl.db but I am sure it's not /reference/ instead either it should be reference/ or <full path>/reference/...

Change all the paths carefully in the configuration file.

ADD COMMENTlink modified 2.8 years ago • written 5.1 years ago by geek_y10k
0
gravatar for stu111392
5.0 years ago by
stu11139230
Germany/Kiel
stu11139230 wrote:

Wow, that is a long list. Apart from that I installed SURPI and tryed to start the pipeline. But when I start it there is instantly a file created which contains the following error message:

nohup: ignoring input
DBI connect('dbname=/reference/taxonomy/gi_taxid_nucl.db','',...) failed: unable to open database file at /home/jsusat/surpi-1.0.18/taxonomy_lookup_embedded.pl line 106.
/home/jsusat/surpi-1.0.18/SURPI.sh: line 507: [: =: unary operator expected

Do I have to set up a reference before is start? Or is the taxonomy lookup skript itself fetching a db? I just followed the steps from the website:

http://chiulab.ucsf.edu/surpi/documentation/How_to_run_SURPI.pdf

With kind regards,

Julian

ADD COMMENTlink written 5.0 years ago by stu11139230
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1547 users visited in the last hour