Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • fjackling
    Junior Member
    • Feb 2011
    • 2

    problem with indexing using bfast

    I'm very new to next-gen sequencing and bioinformatics in general. I'm trying to use bfast to align solid4 colourspace reads to a cDNA sequence. I have completed the following steps:

    1) fasta2brg -A 0 -f cDNA.fasta
    2) fasta2brg -A 1 -f cDNA.fasta

    It is when I attemp step 3 that I encounter a problem:

    3) bfast index -A 1 -f cDNA.fa -m 1111111111111111111111 -w 14 -i 1
    bfast index -A 1 -f cDNA.fa -m 111110100111110011111111111 -w 14 -i 2
    bfast index -A 1 -f cDNA.fa -m 10111111011001100011111000111111 -w 14 -i 3
    bfast index -A 1 -f cDNA.fa -m 1111111100101111000001100011111011 -w 14 -i 4
    bfast index -A 1 -f cDNA.fa -m 111111110001111110011111111 -w 14 -i 5
    bfast index -A 1 -f cDNA.fa -m 11111011010011000011000110011111111 -w 14 -i 6
    bfast index -A 1 -f cDNA.fa -m 1111111111110011101111111 -w 14 -i 7
    bfast index -A 1 -f cDNA.fa -m 111011000011111111001111011111 -w 14 -i 8
    bfast index -A 1 -f cDNA.fa -m 1110110001011010011100101111101111 -w 14 -i 9
    bfast index -A 1 -f cDNA.fa -m 111111001000110001011100110001100011111 -w 14 -i 10

    The error message I get is:
    ************************************************************
    Checking input parameters supplied by the user ...
    Validating fastaFileName cDNA.fa.
    Validating tmpDir path ./.
    Input arguments look good!
    ************************************************************
    ************************************************************
    Printing Program Parameters:
    programMode: [ExecuteProgram]
    fastaFileName: cDNA.fa
    space: [Color Space]
    mask: 1111111111111111111111
    depth: 0
    hashWidth: 14
    indexNumber: 1
    repeatMasker: [Not Using]
    startContig: 0
    startPos: 0
    endContig: 2147483647
    endPos: 2147483647
    exonsFileName: [Not Using]
    numThreads: 1
    tmpDir: ./
    timing: [Not Using]
    ************************************************************
    ************************************************************
    Reading in reference genome from cDNA.fa.cs.brg.
    bfast: RGBinary.c:304: RGBinaryReadBinaryHeader: Assertion `rg->contigs[i].contigNameLength > 0' failed.
    /var/spool/torque/mom_priv/jobs/804690.bruce-m.vlsci.unimelb.edu.au.SC: line 10: 4145 Aborted bfast index -A 1 -f cDNA.fa -m 1111111111111111111111 -w 14 -i 1
    ************************************************************

    Is anyone able to suggest what I might be doing wrong? Any suggestions or assistance would be greatly appreciated.
  • nilshomer
    Nils Homer
    • Nov 2008
    • 1283

    #2
    It looks like there is a problem with your input FASTA. Do you have empty lines or empty contigs? Could you post it somewhere for us to take a look?

    Comment

    • fjackling
      Junior Member
      • Feb 2011
      • 2

      #3
      Ah, yes I had a gap in my fasta file between the > and the sample name but I have fixed this up and repeated the first 2 steps and now the third step worked. Thanks

      Comment

      Latest Articles

      Collapse

      • mylaser
        Reply to Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
        by mylaser
        Kheloyaar: The Complete Guide to Kheloyaar Loginand Kheloyaar ID
        The online gaming industry has transformed the way people enjoy digital entertainment. As technology continues to improve, players are looking for platforms that offer convenience, security, and a seamless user experience. Kheloyaarhas gained attention among users who value an easy-to-use platform, quick account access, and a simple registration process.
        Whether you're exploring Kheloyaar for the first time or want to understand...
        Today, 09:27 PM
      • SEQadmin2
        Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
        by SEQadmin2



        Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
        ...
        Today, 11:10 AM
      • SEQadmin2
        Cancer Drug Resistance: The Lingering Barrier to Rising Survival
        by SEQadmin2



        Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

        There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
        Yesterday, 05:17 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, Today, 10:04 AM
      0 responses
      8 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, Yesterday, 10:08 AM
      0 responses
      7 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 07-07-2026, 11:05 AM
      0 responses
      9 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 07-02-2026, 11:08 AM
      0 responses
      31 views
      0 reactions
      Last Post SEQadmin2  
      Working...