Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • zslee
    Member
    • May 2009
    • 29

    #1

    Defuse Calculating spanning p-values Failure

    Hey, could anyone help me with one problem:
    when I run defuse to find gene fusions based on fastq data (illumina HISEQ, paired-end, trimed by trimglore), the software goes well until Calculating spanning p-values:

    Importing fastq files
    Splitting fastq files
    Discordant alignments
    Read Stats
    Fragment mean 182.171678030117 stddev 43.6861609169054
    Read length min 20 max 100
    Generating discordant alignment clusters
    Remove mitochondrial-genomic clusters
    Generating maximum parsimony solution
    Selecting fusion clusters
    Preparing sequences for local realignment
    Performing local realignment
    Filtering concordant clusters
    Generating spanning alignment regions file
    Initializing split read alignments
    Calculating split read alignments
    Evaluating split reads
    Calculating spanning stats
    Calculating spanning p-values
    Failure for defuse command:
    /home/local/bin/Rscript ./defuse-0.6.1//scripts/evaluate_fraglength_mean.R ./temp/defuse/concordant.read.stats ./temp/defuse/spanlength.cov 50 ./temp/defuse/splitreads.span.stats ./temp/defuse/splitreads.span.pval.tmp
    Reason:
    Job command with nonzero return code
    Job output:
    Running on node4
    Rscript execution error: No such file or directory

    real 0m0.099s
    user 0m0.000s
    sys 0m0.003s
    Return codes: 255
    Commands failed after 1 seconds


    when I check I found there is no ./temp/defuse/splitreads.span.pval.tmp produced, what is the reason?
    thanks in advance for anyone's reply.

    code as bellow:
    ~/programs/defuse-0.6.1/scripts/defuse.pl -c ~/programs/defuse-0.6.1/scripts/config.txt -1 $inputdir/$file1 -2 $inputdir/$file2 -o $outputdir -p 24
  • NKAkers
    Member
    • Sep 2011
    • 26

    #2
    /home/local/bin/Rscript ./defuse-0.6.1//scripts/evaluate_fraglength_mean.R

    This is your problem, I believe. The double forward slash betweeen defuse-0.6.1 and scripts.

    Check your config.txt file to make sure you named your directories without a slash at the end. ie:

    /home/usr/defuse-0.6.1
    NOT:
    /home/usr/defuse-0.6.1/

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
      by SEQadmin2



      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
      ...
      10-14-2026, 09:55 AM
    • SEQadmin2
      New Genomics Technologies Take Aim at Long-Standing Limits
      by SEQadmin2


      Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

      We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
      ...
      09-28-2026, 10:25 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    13 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    20 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-13-2026, 09:56 AM
    0 responses
    19 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-09-2026, 11:14 AM
    0 responses
    27 views
    0 reactions
    Last Post SEQadmin2  
    Working...