Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • valei
    Member
    • Sep 2010
    • 14

    #1

    TopHat2 with genes lacking introns

    Dear all,
    I'm working with RNA-seq paired-end data and I'm using TopHat-BowTie in order to align my reads against the Human reference genome.
    My problem is the following: I' have some reads coming from genes that are duplicated in the genome and they are very similar in terms of sequence, the only difference between them is that the second gene doesn't have the introns (for example this is the case of GLUD and GLUD2, or NACA and NACA2, etc..). When TopHat works with these reads it prefer align with mismatches in the contiguous sequence (GLUD2 or NACA2) instead of splitting the read in two parts and align it in two separated exon without any mismatch (GLUD or NACA).
    Obviously the correct alignment is the second one, I'm sure of this because I got the same mismatches in all the samples that I analyzed.
    How can I overcome this problem?
    Can I use some option or particular parameter? I would like to constrain TopHat to prefer no-mismatch alignments even if they are splitted in two of more segments.
    I hope my question is clear....
    Many many thanks for your help!
    Last edited by valei; 06-08-2012, 02:16 AM.
  • Dario1984
    Senior Member
    • Jun 2011
    • 166

    #2
    You have to upgrade to TopHat version 2. You supply your transcriptome annotation with the --GTF option. It maps the reads first to the spliced transcript models, which will take away your problem.

    The name for the genes that don't have introns and are the same sequence as ones that do is pseudogenes.

    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...