Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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

  • #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

    • seqadmin
      Advancing Precision Medicine for Rare Diseases in Children
      by seqadmin




      Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
      12-16-2024, 07:57 AM
    • seqadmin
      Recent Advances in Sequencing Technologies
      by seqadmin



      Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

      Long-Read Sequencing
      Long-read sequencing has seen remarkable advancements,...
      12-02-2024, 01:49 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 12-17-2024, 10:28 AM
    0 responses
    33 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-13-2024, 08:24 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-12-2024, 07:41 AM
    0 responses
    34 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-11-2024, 07:45 AM
    0 responses
    46 views
    0 likes
    Last Post seqadmin  
    Working...
    X