Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Crash when no fusions are present #84

Open
nrhorner opened this issue Jul 26, 2022 · 0 comments
Open

Crash when no fusions are present #84

nrhorner opened this issue Jul 26, 2022 · 0 comments

Comments

@nrhorner
Copy link

Hi,

Using JAFFA v2.2.

I run JAFFAL and get a several fusion transcripts reported as expected.
In order to test what happens if no fusion transcripts are identified, I removed the previously found reads and rerun with the following command.

`bpipe run -n 4 -p jaffa_output="$JAFFAOUT/" -p refBase=chr20 -p genome=hg38_chr20 -p annotation=genCode22 -p fastqInputFormat="*.fastq" JAFFAL.groovy nofusion_genes_full_length_reads.fastq

`╒══════════════════════════════════════════════════════════════════════════════════════════════════╕
| Starting Pipeline at 2022-07-26 12:04 |
╘══════════════════════════════════════════════════════════════════════════════════════════════════╛

========================================= Stage run_check ==========================================

====================================== Stage get_fasta (all) =======================================

================================ Stage minimap2_transcriptome (all) ================================

================================== Stage filter_transcripts (all) ==================================

=============================== Stage extract_fusion_sequences (all) ===============================

=================================== Stage minimap2_genome (all) ====================================

================================ Stage make_fasta_reads_table (all) ================================

==================================== Stage get_final_list (all) ====================================

================================ Stage report_3_gene_fusions (all) =================================

==================================== Stage compile_all_results =====================================

options(echo=F)
[1] "No fusions were found for the following samples:"
[1] "all"
Compiling the results from:
all
Error in file(file, "rt") : cannot open the connection
Calls: read.delim -> read.table -> file
In addition: Warning message:
In file(file, "rt") : cannot open file 'NA': No such file or directory
Execution halted
ERROR: stage compile_all_results failed: Command in stage compile_all_results failed with exit status = 1 :


Is there a way to prevent this happening?
Thanks, Neil

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant