Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • CREST outputs: what's wrong?

    Spent many days to get all the CREST dependencies to run, and anyone knows what this means?

    Code:
    CREST.pl -f tumor.bam.cover.somatic.cover -d tumor.bam -g normal.bam --ref_genome human_g1k_v37_decoy.fasta -t human_g1k_v37_decoy.2bit --blatserver localhost --blatport 2014 -r 1
    I got a bunch of stuff like:

    Code:
    SV filter starting....
    low compexity filter
    -seq value = query.fa
    -repsiz range value = 2,8
    -repsiz values = 2 to 8
    -minrep value = 4
     INFO: Program read 167 bases 
     total bases read=167 
    Sorting 2 repeat records
    End 1st sort
    Filtering repeats
    nreps=2
    AFTER FILTERS nreps=2
     End of Sequence file!
     total bases read=0 
     End of nucleotides. Exiting!
    Type distance filter
    1       223254531       +       hs37d5  33068960        -       14      hs37d5  33068960        +       1       223254531       -       19      CTX     0       0
    Germline sclip filter
    Loaded 201 letters in 1 sequences
    Searched 1260 bases in 47 sequences
    FAILED
    SV filter starting....
    low compexity filter
    -seq value = query.fa
    -repsiz range value = 2,8
    -repsiz values = 2 to 8
    -minrep value = 4
     INFO: Program read 106 bases 
     total bases read=106 
     End of Sequence file!
     total bases read=0 
     End of nucleotides. Exiting!
    Type distance filter
    Germline sclip filter
    Loaded 201 letters in 1 sequences
    Searched 398 bases in 17 sequences
    FAILED
    tumor.bam.predSV.txt file is empty.

    What's wrong here? Any idea?

  • #2
    hello,I have met the same question. The .predSV.txt is empty. Do you have solved it? Can you help me?

    Comment


    • #3
      Originally posted by gaoyf View Post
      hello,I have met the same question. The .predSV.txt is empty. Do you have solved it? Can you help me?
      I forgot what my original problem was.
      The .predSV.txt file is pretty small, sometimes only ~ 50 lines for the whole genome.
      Have you tried the whole genome and saw an empty file?

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Essential Discoveries and Tools in Epitranscriptomics
        by seqadmin




        The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
        04-22-2024, 07:01 AM
      • 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

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 04-25-2024, 11:49 AM
      0 responses
      20 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-24-2024, 08:47 AM
      0 responses
      20 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-11-2024, 12:08 PM
      0 responses
      62 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-10-2024, 10:19 PM
      0 responses
      61 views
      0 likes
      Last Post seqadmin  
      Working...
      X