Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Sow
    Member
    • Feb 2016
    • 16

    Bowtie/1.1.2 Encountered reference seq with gaps

    Hello,

    Can anyone decipher this code for me..

    } else if(cat >= 2) {
    if(lc != -1 && off == 0) off++;
    lc = -1;
    off++; // skip over gap character and increment
    } else if(c == '>') {
    if(off > 0 && lastc == '>') {
    cerr << "Warning: Encountered reference sequence with only gaps" << endl;
    } else if(lastc == '>') {
    cerr << "Warning: Encountered empty reference sequence" << endl;
    }
    lastc = '>';
    //return RefRecord(off, 0, false);
    return RefRecord((uint32_t)off, 0, first);
    }
    c = in.get();
    if(c == -1) {
    // End-of-file
    if(off > 0 && lastc == '>') {
    cerr << "Warning: Encountered reference sequence with only gaps" << endl;
    } else if(lastc == '>') {
    cerr << "Warning: Encountered empty reference sequence" << endl;
    }
    lastc = -1;
    //return RefRecord(off, 0, false);
    return RefRecord((uint32_t)off, 0, first);
    }
    }
    assert(!rparms.color || (lc != -1));
    assert_eq(1, asc2dnacat[c]); // C must be unambiguous base
    if(off > 0 && rparms.color && first) {
    // Handle the case where the first record has ambiguous
    // characters but we're in color space; one of those counts is
    // spurious
    off--;
    }

    I'm trying to build a bowtie index on a genome file but it returns the warning
    "Encountered reference sequence with gaps".

    I'm not sure what that means. My genome is not Rep Masked to have N (at least that I know of). And when I do downstream aligning only 20% of my reads aligns with the genome. When I created an index for an older version of the genome I did not encounter these warnings and I had more reads aligning to the genome (more than 50%). I am limited to using the new version of the genome since I am using a pipeline and the gff file I use later in the pipeline have the same headers as the new version of the genome.

    My question is
    1. What does the warning "encountered....." mean and how do I solve the issue?
    2. Can the low percentage of reads aligning to the genome be because of these warnings?

    Thanks in advance for your help!

Latest Articles

Collapse

  • mylaser
    Reply to Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
    by mylaser
    Kheloyar – Everything You Need to Know About Kheloyaar Login and Kheoyar Id
    If you are looking for an online gaming platform that offers a user-friendly experience, Kheloyar has become a name that many users search for. Whether you're interested in creating a new account, accessing your dashboard through Kheloyaar Login, or learning how to obtain a Kheoyar Id, understanding the platform's features and account process is essential.
    This guide explains everything you need to know about...
    Yesterday, 01:13 AM
  • 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.
    ...
    07-09-2026, 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...
    07-08-2026, 05:17 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 07-09-2026, 10:04 AM
0 responses
17 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-08-2026, 10:08 AM
0 responses
10 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-07-2026, 11:05 AM
0 responses
24 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...