Troubleshooting

Problems with setup or installation

Please leave us an issue and we will try to help out.

No reads found

• Check your reads contain _R1.fastq and _R2.fastq in the read name.
• Reads cannot be gzipped at present (for the read quality filtering tool prinseq)

If you get missing or 0 sized empty bams in output of wochenende

• Check most recent work folder
• cd work && ls -lrt
• Go into most recent folder and subfolder
• ls -l
• cat .command.err
• It could be that trimming adapter files or the program perldup have not been located correctly
• Check your config.yml that the locations to these tools are correct.
• Also check the aligner bwa is accessible within the wochenende conda environment. Type conda activate wochenende && bwa

Reference problems

If you want to use your own reference sequence or add to our offerings, follow the guide on this Wiki at https://github.com/MHH-RCUG/nf_wochenende/wiki/Building-a-reference-sequence for best results. Contact us via issue if problems persist. Our references have been optimised for human lung samples from diseased and healthy patients, but have been useful for freshwater sampling, mouse guts, contaminated mouse filter paper samples, and more.

No output for growth rate

• Considerable numbers of reads are required for good growth rate analysis.
• In practice, this typically means you need a short read metagenome with >500 reads per sample, with more leading to better results
• This limits growth rate analysis to highly covered metagenomes, with only the more common species getting a growth rate estimate.

No output for plots

• Plots might have considered the data too sparse to produce a results.
• Likely the number of reads was too low.
• Plots has mostly been tested on short read metagenomes, but should work fine for long reads.
• Leave us an issue if trouble persists.

No output for Haybaler integration tool

• Haybaler only works with 2+ samples
• We can't do data integration and prepare data for plots from just one sample :-)

Are you having difficulty producing heatmaps or lack an R server in your cluster ?

• Take the outputs from Haybaler, load into your favourite tool ( R, Python, others), and create some visualizations manually
• R code examples are available in our Haybaler repository on github.

Finally - reliability of metagenomics tools

We always recommend using multiple independent tools and methods (and coverage profiles!) of putatively located bacteria before considering the results to be true. Mock communities are useful (also with public data) to test a tool before applying it.