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
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM
    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    32 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    37 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    31 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    53 views
    0 likes
    Last Post seqadmin  
    Working...
    X