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
      Genetic Variation in Immunogenetics and Antibody Diversity
      by seqadmin



      The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
      11-06-2024, 07:24 PM
    • seqadmin
      Choosing Between NGS and qPCR
      by seqadmin



      Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
      10-18-2024, 07:11 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Today, 11:09 AM
    0 responses
    22 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Today, 06:13 AM
    0 responses
    19 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 11-01-2024, 06:09 AM
    0 responses
    30 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-30-2024, 05:31 AM
    0 responses
    21 views
    0 likes
    Last Post seqadmin  
    Working...
    X