Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • lethalfang
    Member
    • Aug 2011
    • 95

    #1

    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?
  • gaoyf
    Junior Member
    • Jun 2015
    • 1

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

    Comment

    • lethalfang
      Member
      • Aug 2011
      • 95

      #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

      • kinpict.com
        Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by kinpict.com
        The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
        Today, 06:37 PM
      • SEQadmin2
        16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by SEQadmin2





        NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
        10-30-2026, 06:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 11-12-2026, 07:59 AM
      0 responses
      27 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 11-05-2026, 09:25 AM
      0 responses
      35 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 11-03-2026, 11:54 AM
      2 responses
      56 views
      0 reactions
      Last Post kinpict.com  
      Started by SEQadmin2, 10-28-2026, 09:57 AM
      0 responses
      175 views
      0 reactions
      Last Post SEQadmin2  
      Working...